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INTRODUCTION 

The promise of filamentary composite materials, whose de- 
velopment may be considered as entering its second generation, 
continues to generate intense interest and applications activ- 
ity. Such interest and activity are well founded, since they 
are based on the possibility of using relatively brittle mate- 
rials with high modulus, high strength, but low density in 
composites with good durability and high tolerance to damage 
and which, when they do fail, do so in a non-catastrophic 
manner. Fiber reinforced composite materials of this kind 
offer substantially improved performance and potentially 
lower costs for aerospace hardware. 

Much progress has been achieved since the initial devel- 
opments in the mid 1960's. Rather limited applications to 
primary aircraft structure have been made, however, mainly in 
a material -substitution mode on military aircraft, except for 
a few experiments now being conducted on large passenger air- 
planes. 

To fulfill the promise of composite materials completely 
requires a strong technology base. NASA and AFOSR recognize 
the present state of the art to be such that to fully exploit 
composites in sophisticated aerospace structures, the tech- 
nology base must be improved. This, in turn, calls for ex- 
panding fundamental knowledge and the means by which it can 
be sucessfully applied in design and manufacture. 
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As technology of composite materials and structures moves 
toward fuller adoption into aerospace structures, some of the 
problems of an earlier era are being solved, others which 
seemed important are being put into perspective as relatively 
minor, and still others unanticipated or put aside are emerg- 
ing as of high priority. The purpose of the RPI program as 
funded by NASA and AFOSR has been to develop critical advanced 
technology in the areas of physical properties, structural 
concepts and analysis, manufacturing, reliability and life 
prediction. 

Our approach to accomplishing these goals is through an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constituents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their non-destructive and proof testing, to the logical 
conclusion of such activities; namely research into the com- 
posite structure's long term integrity under conditions perti- 
nent to service use. Inherent m the RPI program is the moti- 
vation which basic research into the structural aspects pro- 
vides for research at the materials level, and vice versa. 

Second, interactions among faculty contributing to pro- 
gram objectives - which is a group wider than that supported 
under the project - is on a day to day basis, regardless of 
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organizational lines. Program management is largely at the 
working level, and administrative, scientific and technical 
decisions are made, for the most part, independent of con- 
siderations normally associated with academic departments. 
Involvement of this kind includes - depending on the flow of 
the research - faculty, staff and students from chemistry, 
chemical engineering, civil engineering, materials engineer- 
ing, and the department of mechanical engineering, aeronau- 
tical engineering and mechanics. 

Both of these characteristics of the NASA/AFOSR program 
of research in composite materials and structures foster the 
kinds of fundamental advances which are triggered by insights 
into aspects beyond the narrow confines of an individual dis- 
cipline. This is a program characteristic often sought in 
many fields at a university, but seldom achieved. 

Overall program emphasis, is on basic, long-term re- 
search in the following categories: (a) constituent mater- 

ials, (b) composite materials, (c) generic structural ele- 
ments, (d) processing science technology and (e) maintaining 
long-term structural integrity. Progress in the program will 
be reported in the following pages under these headings. 

Those computer software developments are also undertaken which 
both support Rensselaer projects in composite materials and 
structures research in the areas listed above and which also 
represent research with the potential of widely useful re- 
sults m their own right. 



6 


In short, the NASA/AFOSR Composites Aircraft Program is 
a multi-faceted program planned and managed so that scien- 
tists and engineers in a number of pertinent disciplines will 
interact to 'achieve its goals. Research in the basic compo- 
sition, characteristics and processing science of composite 
- materials and their constituents is balanced against the 

mechanics, conceptual design, fabrication and testing of gen- 
eric structural elements typical of aerospace vehicles so as 
to encourage the discovery of unusual solutions to present 
and future problems. In the following sections, more de- 
tailed descriptions of the progress achieved in the various 
component parts of this comprehensive program are presented. 
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transverse properties of anisotropic fibers) 



8 




9 


iv/i. 


n ur 


Ji I 


:J Preceding page blank 


I I -A TRANSVERSE PROPERTIES OF ANISOTROPIC FIBER REINFORCED 
EPOXY MATRIX COMPOSITES (A step in determining the 
transverse properties of anisotropic fibers) 

Senior Investigator: R. J. Diefendorf 

1. Introduction 

Unlike the longitudinal properties, the transverse 
properties of composites are a function of packing and homo- 
geneity. If, as is often the case, the transverse proper- 
ties dictate or limit composite utilization, then homogene- 
ity and packing take on added importance. Inhomogeneity 
produces critical resin areas which may induce crack propa- 
gation (see Figure II-A-1) . Therefore, a measure of packing 
uniformity is essential for describing transverse properties. 

Several researchers have attempted to make provisions 

f 11 * 

for homogeneity and packing. Jones and Charms and Sen- 

r 2 1 

deckyj use a semi-empirical method and introduce a term 
called a contiguity factor, C, where C = 0 for areas where 
one fiber is isolated from others, and C = 1 for areas where 
matrix material is surrounded by fibers (see Figure II-A-2) . 
Values between 0 and 1 are determined by comparison of the 
composite cross-section to experimental results (i.e., prop- 
erties) and interpolating. For either the case of C = 1 or 

C = 0, there is no provision for distinguishing between the 

r 3 1 

effects of a square or hexagonal array. Tsai 1- J includes 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 64 . 
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Figure II-A-1. Schematic Diagram Depicting Various Resin Rich Geometries 
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Figure II-A-2. Schematic Diagram Depicting Various Packing Geometries 
for Continuity Factor, C = 0 and C = 1 
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contiguity in an elasticity approach in order to obtain a 
value for the composite modulus transverse to the fiber axis. 
As can be seen from Equations (1) through (5) , the formula- 
tion is somewhat cumbersome. 


E„ 


2[i~ v + (v -v lv 3 

f r m n 


( 1-0 


K (2K + G ) 
f m m 


G (K - K )V 
m f mm 


(2K + G ) + 2 (K - K ) V 
mm f m m 


K (2K 4- G ) + G„(K -K_)V 
f m f f m f m 


(2K + GJ - 2 (K — K ) V 
m r m t m 


( 1 ) 


where 


Kj- = 


E f 

2(1 - v f ) 


E 


K = 
m 


m 


2(1 -v 


E, 


Gj- - 


2(1 + v f ) 


E 


G m " 


m 


2 (! + v ) 
m 


(2) 

(3) 

(4) 

(5) 


For C = 1 the first term in brackets drops out, while 
the same is true for the second term when C - 0 . Intermedi- 
ate values were determined by comparing curves of E 2 versus 
volume fraction. A shortcoming of the formulation is that 
it makes the assumption of isotropic constituents. 

Another formulation that considers geometric configura- 
tion, to a somewhat limited degree, is the Halpin-Tsai equa- 

T41 ’ r 5 1 

tions . These equations are based on Herman " s 
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r 6 1 

generalized solutions of Hill's model 1 1 and can be reduced 
to the approximate forms given below. 


e 2 1 + enV f 



(E f /E) - 1 
f m 



( 6 ) 

(7) 


where e is related to fiber geometry, packing geometry and 
load conditions. For circular fibers in a square array, 
e = 2 and r] = 3, and a value of ri = 0 indicates a homo- 
geneous material. Difficulties arise in determining appro- 
priate values for e, since e can vary from zero to infinity. 
Values are usually determined via curve fitting and by com- 
paring Equations (6) and (7) to exact elasticity solutions. 
However, Halpin and Tsai were able to predict with high ac- 
curacy the properties of a composite with 55% reinforcing 
(see Figure II-A-3) . This theory also assumes isotropic con- 
stituents and does not account for fiber anisotropies. An 
interesting concept is the theory of fiber spacing, 6/r. 
However, the model only considers regularly spaced (i.e. , 
non-random) arrays, since 6 in the X direction is always e- 
qual to 5 in the Y direction. A modification which would 
allow analysis of random geometries - may be possible and would 
be highly beneficial in including geometry in composite prop- 
erties calculations. 
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Figure II-A-3. Comparison of Halpm-Tsai Calculations (circles) 

with Adams and Doner's Calculations for E of Cir- 
cular Fibers m a Square Array 

(after Halpin and Tsai, Ref. 4) 
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r 7 1 

An example of packing effect is given by Piggot . In 
this work, packing is related to shear modulus with the 
result being: 

r 

U R - U f = —~Hn (R/r) (8) 

m 

where U_. is the nearest neighbor fiber displacement under 

K. 

load, is the fiber displacement under load and the quan- 
tity R/r is the argument dependent upon packing. 

The respective expressions for a square and hexagonal 
array are: 


Jin (R/r) = hln(ir/V € ) 

s n 

(9) 

Jin (R/r) ^ - *5 Jin (2'rr//3V^) 

(10) 


For a volume fraction of 50%, Equation (9) yields a value of 
Jin (R/r )‘ equal to .92 for the square array and Equation (10) 
a value of .99 for the hexagonal array. These equations can 
be generalized to yield the following results: 

In (R/r) = k&n (P f /V f ) (11)* 

where 

P_g = geometrical packing factor . 

The difference in packing factor between square and hexagonal 
arrays is .07. While this may not be significant for the 


* 

NOTE: For a more detailed account, the reader is referred 

to the material in Reference [7] . 
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shear modulus, two points must be mentioned. First, packing 
may have an important effect on failure. Second, the effect 
of random geometry has not been fully described and quanti- 
fied . 


2. Status 

In order to measure homogeneity in a laminate, a speci- 
men was cut from a plate which was fabricated from 27 layers 
of graphite prepreg using a standard 250°F cure at a pressure 
of 200 psi with 8 bleeders. The specimen was prepared as a 
metallographic sample, and photomicrographs were taken. 

Thepe pictures were analyzed with the ASTM's grid technique 
and the Bausch & Lomb image analysis system, to determine 
the sample's fiber volume fraction. 

As a first step towards describing homogeneity, the im- 
age analysis system was used to perform subfield analysis on 
the sample. That is, the area of scan was subsequently sub- 
divided from whole field to halves, quarters., eighths, six- 
teenths and finally thirty-seconds. Values of volume frac- 
tion were obtained for each set of subfields so as to gain 
some insight into the variation of volume fraction across 
the sample and, therefore, the "variation of homogeneity". 
From this data a plot of volume fraction versus area was 
obtained. A plot of variance versus reciprocal area was 
also constructed so as to determine the machine's effect in 
subfield analysis. 
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The variance versus volume fraction was determined, us- 
ing data from theoretical arrays of hexagonal (48.0%) packing 
and square (34.9% and 66.1%) packing. A parabolic relation- 
ship should exist between variance and volume fraction, with 
a minima occuring around 50% fibers. The reason for the 50% 
minima was attributable to the fact that even amounts of 
fiber and matrix material are spaced uniformly, thus result- 
ing in a drop in variance. Experimental data, as seen in 
Figure II-A-4, validated the hypothesis. 

3. Progress During Report Period 
The work reported here is largely that of graduate stu- 
dent J. F. Helmer. 

a. Transverse Composite Homogeneity 

Theoretical arrays of hexagonal and square geometries 
were constructed, and volume fraction versus area was meas- 
ured. The purpose was twofold. First, the data from the 
arrays was needed to understand the "best case" limit of 
fiber packing. Second, the information was useful in evalu- 
ating actual composite performance via comparisons. The 
next step was to characterize a real composite by means of 
comparison to the square and hexagonal arrays . 

The volume fractions of the square and hexagonal arrays 
were 66.1% and 48.0%, respectively. As was the case for the 
composite specimen, the two theoretical arrays were analyzed 
for volume fraction at full, half, quarter, eighth, sixteenth 
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VOLUME FRACTION % 

Figure II-A-4. Variance versus Volume Fraction - Theoretical Array 
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and thirty-second divisions. Averaged data for groups of 
ten measurements is presented for the square array in Table 
II— A— 1 . Table II-A-2 is a comparison of the results for a 
square array and the composite sample. It can be seen that 
the dispersion, as expected, is far less for the regular, 
square array than it is for the composite. This is also 
evident in Figures II-A-5 and -6, which are plots of volume 
fraction versus area for the square array and for the com- 
posite and square array, respectively. Even in the case of 
the hexagonal array (.shown m Figure II-A-7) the theoretical 
dispersion is less than that of the composite, so -that the 
composite dispersion may be used as a measure of the degree 
of overall inhomogeneity. 

Plots of standard deviation of volume fraction versus 
average volume fraction were constructed for the composite 
and for the square array (.see Figures II-A-8a and -8b) . 

Since the slopes have different signs and the scattering of 
a few of the data points determines the slope, it was con- 
cluded that the standard deviation was independent of volume 
fraction over the range studied. An estimate of technique 
error can be obtained by constructing a plot of frequency 
versus deviation from mean including all 320 data points for 
the 1/3 2nd subfield analysis. That is, taking the average 
of each individual data set and subtracting each of the ten 
data points in the set to obtain the individualized data 
points. Figure II-A-9a is a plot of this data for the 
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TABLE II-A-1 

VOLUME FRACTION DATA FROM SUBFIELD ANALYSIS 


Subfield 

Divisions 

32 

nds 

lG ths 

g ths 

4 ths 

Halves 

Full 


68.1 

67.7 

68.6 

66.3 

68.6 

67.6 

68.3 


66.7 

68.4 

68.2 

68.3 

68.7 

86.8 

— 


68.6 

68.8 

70.3 

68.5 

67.2 

— 



67.5 

68.1 

69.6 

67.0 

68.7 




67.2 

67.7 

68.3 

67.1 

— 




68.3 

68.1 

71.0 

69.2 




Volume 

67.9 

67.4 

66.3 

67.6 




Fraction 

m 

67.9 

68.8 

70.1 

67.8 




Percent 

66.2 

68.1 

65.6 

— 





65.8 

66.3 

65.5 






69.2 

69.6 

67.4 






66.9 

69.3 

69.3 






68.4 

70.1 

67.0 






67.4 

67.9 

70.0 






67.0 

74.0 

68.6 






66.3 

71.7 

70.2 





Mean 

68. 

2 

68.5 

67.7 

68.3 

66.7 

68.3 

Standard 

Deviation 

1 . 

612 

1.730 

0.932 

0.735 

1.273 

— 

Variance 

2. 

518 

2.806 

0.759 

0.405 

0.81 
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TABLE II-A-2 

COMPARISON OF SUBFIELD DATA FOR COMPOSITE AND SQUARE ARRAY 


Subfield 

Divisions 

32 nds 

16 thS 

8 th5 

. ths 
4 

Halves 

Composite 






Mean 

64.2 

63.9 

64.1 

66.4 

66.8 

Standard 

Deviation 

6.325 

5.003 

1.842 

1.771 

2.616 

Variance 

38.754 

23.465 

2.970 

2.353 

3.423 

Square Array 






Mean 

68.2 

68.5 

67.7 ' 

68.3 

66.7 

Standard 

Deviation 

1.612 

1.730 

0.932 

0.735 

1.273 

Variance 

2.518 

2.896 

0.759 

0.405 

0.81 


Full 


66.0 


68.3 






VOLUME FRACTION % 


K> 



32/N 

N- NUMBER OF FIELDS 

Figure XV-A-7. Volume Fraction of Fibers versus Area for Theoretical Hexagonal Array 
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Figure II-A-8b. Standard Deviation versus Volume Fractions 32nd Subfield Analysis - Square Array 
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Figure II-A-9a. Frequency versus Normalized Deviation, 320 Points: 32nd Subfield Analysis - Square Array 
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composite, while Figures II-A-9b and -9c are similar plots 
for the square array. Again, both data sets appear to be 
normally distributed, but in this case, the breadths should 
be the same. The breadth is a measure of the standard devi- 
ation in the technique. 

Further insight can be gained by constructing a plot of 
frequency versus volume fraction for the thirty-second sub- 
field analysis for both the composite and the square array. 
Such are shown in Figures II-A-10a and -10b. As can be as- 
certained by inspection, the distribution of data for both 
the square array and composites are quite Gaussian. The 
distribution for the composite, however, while fairly normal, 
is skewed towards low volume fraction or resin rich areas 
(see Figure II-A-10c) . 

It is possible to plot volume fraction dispersion ver- 
sus area for a periodic array, in terms of volume fraction 
versus length in unit cell dimensions. Such a curve, in 
theory, takes the shape of Figure II-A-11, assuming no in- 
strumental error. The dispersion of the measured volume 
fractions is shown to decrease with increasing number of 
unit cell areas, since the relative importance of that part 
of the subfield area which is not an integer multiple of a 
unit cell decreases with increasing viewing area. In fact, 
as the number of unit cells gets very large, the dispersion 
bands will asymptotically approach zero, since the frac- 
tional unit cells are negligible in comparison to the area 
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Figure II-A-9b. Frequency versus Normalized Deviation, 320 Points: 32nd Subfield Analysis - Square Array 
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Figure II~A~9c. Frequency versus Normalized Deviation, 320 Points: 32nd Subfield Analysis - Square Array 
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Figure II-A-lQa. Frequency versus Volume Fraction: 32nd Subfield Analysis - Composite Data 
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Figure II-A-10b. Frequency vs. Volume Fraction: 32nd Subfield Analysis - Square Array Data 
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Figure II-A-lOc. Normal Curves {based on data for frequency m terms of probability, 32 points 
each curve) versus Volume Fraction for the 32nd Subfield Analysis Comparison 
of Composite versus Square Array 
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Figure II-A-11. Volume Fraction Dispersion versus Area 
for a Periodic Array 
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Figure II-A-12. Comparison of Volume Fraction versus 

Number of Unit Cells for Various Arrays 
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of view. When the field of view is exactly an integer mul’- 
tiple of a unit cell, the dispersion bands should also be 
zero. 

Figure II-A-12 is a schematic close-up of the composite 
C [A] ) , periodic square array ([B]) and theoretical square 
array ( [C] ) curves. The experimental square array curve is 
observed to be broader than the corresponding theoretical 
curve due to machine broadening. This is similar m nature 
to the effect encountered in crystallography and x-ray 
studies. In particular, errors are attributable to the 
pixel count of the machine. In the case of the composite 
curve, assuming no machine error, inhomogeneity would in- 
crease the dispersion at the nodal points and smooth out the 
modulation. If the inhomogeneity were at the level of a 
fraction of a unit cell ( ID] ) , the modulation of the dis- 
persion in would smooth, but not drastically increase. 

For large inhomogeneities, the disperson increases, as well 
as smoothing to mono tonic behavior. 

b. Resin Characterization 

Since the overall major goal of the project is the pre- 
diction of transverse fiber properties via derivation from 
the transverse composite properties, the properties of the 
matrix are a major consideration. To this end, and for the 
purpose of determining the aging effects in resins used in 
the ultralight glider, the RP-1 (see Part V-A of this 
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report) , several resin systems were tested for glass tran- 
sition as well as strength data. 

The resin systems were chosen for their properties and 
their widespread usage in the composites laboratory and 
the composites structural program. The systems picked were 
Ciba Giegy epoxies 508 and 509 mixed with Ciba Giegy hard- 
eners XU-224 and XU-225 and a Shell system consisting of 
epon 828 resin mixed with epon V-40 curing agent. For the 
case of the Ciba Giegy resin, one of the epoxies was mixed 
with varying amounts of hardeners. The hardener itself was 
a mixture of XU-224 and XU-225 which varied in the range 
from 3-to-l to 12-to-l, while the ratio of epoxy to harden- 
ers was varied from 3-to-l to 4-to-l. In the Shell system, 
the percentage of the V-40 curing agent ranged from 40% to 
60%, in steps of 5%. The rationale behind varying the ratio 
of mixtures is to determine the best matrix system for trans- 
verse properties research as well as supplying data for the 
processing science experiments described in Part V-B. The 
cure time for the V-40 was approximately 24 hours, but sam- 
ples were not tested until three days had passed. The 50 9- 
based system required approximately three weeks to cure, 
while the 508 system did not totally cure. Even after three 
months, the 508 system was still tacky and could easily be 
plastically deformed. All cures were made at room tempera- 
ture in an open environment without vacuum. Samples, in disc 
shape, were not stored in a dessicator so as to reflect with 
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greater accuracy the actual operating conditions encountered 
by the components of the RP-1- sailplane. Some of the par- 
tially cured 508 system samples were postcured at 212 °F for 
between one to five days in a furnace. The result in each 
case was still a rubber-like material with the only obvious 
effect of differing postcure times being a deepening of 
color with increasing time. Due to these characteristics, 
the 508 system was dropped from further study. 

While it was apparent that the glass transition temper- 
atures for the 508 systems were at or below room temperature, 

it was felt that the T information for the other systems 

<3 

would be useful in characterizing them. The results of work 
performed using several samples on a differential scanning 
calorimeter (DSC) indicated that the T for the 508 system 
also occurred at or below room temperature. A verification 
run of one postcured sample validated the "room temperature 
or below" T^ hypothesis for the 508 system. 

In addition to thermal testing, the discs were measured 
for diameter and thickness at the center and then mechani- 
cally tested to further determine the optimal system. In 
these tests, edges of the specimens rested around their com- 
plete circumferences on a ledge 1/8 inch wide and the load 
was applied to their centers by means of a blunt rod. The 
various specimen data and test results are listed in Table 
II-A-3. In an attempt to estimate the long-term behavior of 
the matrix material used in the structure of the RP-1 



TABLE II-A-3 


MODIFIED SHORTBEAM SHEAR TEST RESULTS FOR SEVERAL COMPOSITE RESIN SYSTEMS 


co 

CO 


Sample 

Number 

Resin Systems 
509 224 225 

Dia. 

(in.) 

Thick, 
(in. ) 

Scale 

(load) 

P yld. 

(lb.) 

P 

ult . 
(lb.) 

6 and 
P 

y 

(in.) 

S = K. w/t 2 

max 1 

(psi) 

Kwa 2 

Y t 3 

max 

(psi) 

1 

24 

9 

1 

1.508 

.103 

.25 

122 

126 

.0394 

11,822 

474,643 

2 

36 

9 

1 

1.487 

.168 

.25 

161 

NA 

.0787 

5,864 

72,267 

3 

36 

12 

1 

0,1.516 

.141 

.25 

98 

371 

.0197 

5,067 

297,247 

4 

36 

9 

1 

1.504 

.189 

.25 

217 

NA 

.0344 

6,245 

156,506 

5 

32 

8 

1 

1.500 

— 

.25 

153 

472 

.0394 



6 

32 

8 

1 

1.498 

.181 

.25 

294 

634 

.0492 

9,225 

168,796 

7 

30 

10 

1 

1.517 

'll. 242 

.25 

444 

0.5 

.0443 

7,794 

118,454 

8 

27 

9 

1 

1.511 

0,. 185 

.5 

273 

616 

.0443 

8,200 

163,027 


9 

33 

11 

1 

1.515 

'v. 188 

.5 

224 

560 

.0443 

6,525 

127,463 


10* 

24 

8 

1 

1.506 

.198 

.5 

No 

Curve 

No Curve 




11* 

30 

10 

1 

1.511 

.213 

.5 

308 

336 

.0344 

6,979 

155,192 

, o q 

12* 

32 

8 

1 

1.499 

.203 

.5 

175 

259 

.0197 

4,366 

177,868 

i C 1 

13* 

33 

11 

1 

1.502 

.169 

.5 

147 

— 

.0148 

5,291 

344,675 

o 

14* 

36 

12 

1 

1.499 

.195 

.25 

266 

— 

.0344 

7,191 

174,676 

; o 

XT' f*** 

15* 

24 

8 

1 

1.510 

.260 

.25 

294 

500 

.0148 

4,471 

189,313 

1 

16* 

24 

9 

1 

1.508 

.185 

.25 

140 

262 

.0148 

4,205 

250,245 

c r.< 

17* 

27 

9 

1 

1,508 

.205 

.25 

318 

343 

.0295 

7,779 

209,584 

, > 

* r- Uj 

18* 

30 

10 

1 

1.502 

.152 

.25 

No 

Curve 

No Curve 




19* 

32 

8 

1 

1.503 

.198 

.25 

224 

— 

.0197 

5,874 

245,358 


20* 

33 

11 

i 

1.506 

.185 

.25 

140 

— 

. 0148 

4,205 

250,245 

1 

21* 

36 

12 

l 

1.506 

.173 

.25 

147 

168 

.0246 

4,415 

193,312 

L 

22 

24 

8 

l 

1.516 

.266 

.25 

595 

1,071 

.0443 

8,645 

119,532 


23* 

24 

8 

l 

1.508 

.246 

. .5 

385 

770 

.0246 

6,540 

176,090 


24 

24 

? 

l 

J. • S1.X 

.166 

.5 

217 

602 

.0394 

8,095 

201,676 


25 

24 

9 

l 

1.512 

.212 

.5 

378 

770 

.0442 

8,646 

150,001 



•k 

post-cured or treated 


[Continued] 



TABLE II-A-3 (Continued) 


Sample 

Number 

Resin Systems 
509 224 225 

Dia. 

(in.) 

Thick, 
(in. ) 

Scale 

(load) 

P yld. 

(lb.) 

P ult. 

(lb.) 

8 and 
P 

y 

(in.) 

S = K w/t 2 

max 1 

(psi) 

r Kwa 2 

Y t 3 

max 

(psi) 

26 

24 

10 

1 

1.513 

.202 

.5 

322 

791 

.0394 

8,112 

166,081 

27 

24 

10 

1 

1.515 

.194 

.5 

329 

714 

.0492 

8,986 

153,405 

28 

27 

9 

1 

1 . 511 

.207 

.5 

266 

595 

.0443 

6,382 

133,392 

29 

27 

9 

1 

1.506 

.200 

.5 

245 

574 

.0394 

6,296 

130.195 

30 

30 

10 

1 

1.506 

.182 

.5 

182 

441 

.0443 

5,648 

114,148 

31 

30 

10 

1 

1.511 

.210 

.5 

259 

588 

..0443 

6,037 

105,743 

32 

32 

8 

1 

1.506 

.187 

.5 

210 

504 

.0541 

6,173 

99,429 

33 

32 

8 

1 

1.504 

.173 

.5 

189 

497 

.0541 

6,492 

113,016 

34 

33 

11 

1 

1.510 

.193 

.5 

217 

511 

.0443 

5,989 

114,130 

35 

33 

11 

1 

1.507 

.196 

.5 

161 

427 

.0492 

4,308 

72,796 

36 

36 

9 

1 

1.507 

.096 

.5 

— 

301 


33,575 


37 

36 

9 

1 

1.498 

.158 

.5 

■ 

329 


13,548 


38 

36 

12 

1 

1.512 

.224 

.5 

196 

497 

.0394 

4,016 

74,136 

39 

36 

12 

1 

1.507 

.136 

.5 

126 

392 

.0492 

7,003 

170,531 

40* 

36 

12 

1 

1.512 

.183 

.5 

328 

679 

.0492 

10,068 

182,209 

41* 

30 

10 

1 

1.513 

.202 

.5 

No Curve No Curve 


— 

42* 

24 

9 

1 

1.513 

.223 

.5 

252 

- — 

.0246 

5,209 

154,727 


* 

post-cured or treated 


u> 


* v -- ^ - 

; OF POOR QU&UI * 



TABLE II-A-3 



Resin Systems 

b ampxe 
Number 

V-40 

828 

43* 

40 

60 

44 

40 

60 

45 

40 

60 

46*() 

40 

60 

47* 

40 

60 

48 

45 

55 

49 

45 

55 

50*- 

45 

55 

51*0 

45 

55 

52 

50 

50 

53 

50 

50 

54*() 

50 

50 

55*() 

50 

50 

56*~ 

50 

50 

57 

55 

45 

58 

55 

45 

59* A 

55 

45 

60*- 

55 

45 

61 

60 

40 

62 

60 

40 

63*() 

60 

40 

64*- 

60 

40 


* 

post-cured 

^no change in thickness 

A 

1 mil change in thickness 

2 mil change in thickness 


Dia. 

Thick. 

Scale 

(in.) 

(m. ) 

(load) 

1.515 

.140 

.1 

1.515 

.137 

.5 

1.512 

.120 

.25 

1.513 

.109 

.1 

1.515 

.148 

.1 

1.510 

.098 

.25 

1.513 

.096 

.25 

1.507 

.078 

.1 

1.508 

.105 

.1 

1.509 

.132 

.1 

1.514 

.181 

.1 

1.512 

.141 

.1 

1.514 

.163 

.1 

1.514 

.162 

.1 

1.512 

.104 

.1 

1.510 

.122 

.1 

1.511 

.109 

.1 

1.511 

. 115 

.1 

1.511 

.095 

.1 

1.511 

.160 

.1 

1.513 

.103 

.1 

1.505 

.076 

.1 



P yld . 
(lb.) 

P , 
ult . 

(lb.) 

6 and 
P 

y 

(in. ) 

S = K,w/t 2 

max 1 

(psi) 

„ Kwa 2 

E = 

Y t 

max 

(psi) 

126 



.0344 

6,609 

223,585 

105 

217 

.0246 

5,751 

278,040 

112 

— 

.0345 

7,996 

314,681 

193 

— 

.0984 

16,699 

253,687 

141 

— 

.0246 

6,617 

296,151 

56 

63 

.0344 

5,994 

289,712 

112 

— 

.0738 

12,493 

287,318 

71 

— 

.0738 

11,997 

339,573 

132 

370 

.0738 

12,308 

258,800 

161 

— 

.0443 

9,499 

264,676 

396 

«-• — 

.0640 

12,426 

174,781 

217 

430 

.0640 

11,220 

202,599 

286 

441 

.0590 

11,066 

187,485 

209 

288 

.0492 

8,187 

167,360 

104 

— 

.0886 

9,885 

174,789 

428 

451 

.3297 

29,561 

119,746 

374 

444 

.2756 

32,360 

175,520 

272 

378 

.3199 

21,143 

93,643 

81 

- — 

.0443 

9,226 

357,211 

153 

— 

.0246 

6,144 

254,338 

168 


.0935 

16,279 

275,423 

94 

— 

.0886 

16,730 

404,826 


o 


ou’g^al :>■ v ; 
of- poo:^ quality 
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sailplane, several of the samples were postcured and/or sub- 
merged in water prior to testing. The purpose was to gain 
added insight into just how brittle the matrix would become 
after prolonged exposure (postcuring) . 

A schematic of the apparatus used in testing the sam- 
ples is shown in Fxgure II-A-13. Since the discs are sup- 
ported at their edges over a radial distance of only an 
eighth of an inch and centrally loaded, the test is actually 
a modified three-point bend test. Further, because of the 
dimensions of the specimens, failure will occur m shear. 
Maximum shear stress and modulus were calculated for a cir- 
cular plate with this kind of loading L 1 . It was assumed 
that the blunt rod applied a uniform load over a circular 
area. The radius of the blunt piston, r Q , was measured and 
used as the radius of the circular load area. The fact that 
the rod end was not perfectly flat was ignored, since the 
disc need only deform a little before the entire radius of 
the rod was in contact with the specimen. 

The values of K and Kj , geometric shape factors , for 
this specific case were determined. Using interpolation for 
an a/r Q value of 2.478 ('a* is radius from center of sample 
to support) yielded values of 1.440 and 1.028 for K and K lt 
respectively. The equations used in the calculations are 
given as follows : 


S 


max 


= KiW/t 2 


( 12 ) 
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Y = Kwa 2 /Et 3 (13) 

max 

where 

S = maximum shear stress 

max 

Y = maximum deflection 

max 

E = isotropic Young's modulus 
w = load 

a = one-half support distance (constant) 
t = specimen thickness. 

Rearranging Equation (13) to solve for modulus yields: 

E - Kwa‘A max t 3 (14) 

Table II-A-3 lists the test results and the calculated values 

for S and E. Table II-A-4 lists special conditions, if 
max 

any, for the discs listed in Table II-A-3. Tables II-A-5 to 
-7 illustrate some of the effects on dimensional stability 
of heat-treatment and moisture absorbance encountered with 
some of the resin systems. 

c. Short Beam Shear Strength/Homogeneity 

The effect of homogeneity on the properties of different 
samples is to be determined using shear strength as a meas- 
ure. This choice was based on the fact that it is easier to 
pick up packing differences in the short beam shear (SBS) 
test than m a tensile test. Two different fabrication 
techniques for the test specimens were employed. The first 
was a typical lay-up (Specimen #20) , consisting of 27 layers 
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TABLE II-A-4 

ENVIRONMENTAL CONDITIONS ON TEST DISCS 


Specimen 

Number 


Heat-Treatment 
Time (hrs.) Temp. C°C) 


10 

48.5 

100 

11 

48.5 

100 

12 

48.5 

100 

13 

48.5 

100 

14 

48.5 

100 

15 

73 

100 

16 

73 

100 

17 

73 

100 

18 

73 

100 

19 

73 

100 

20 

73 

100 

21 

73 

100 

23 

18 

100 

40 

6 

100 

41 

15 

100 

42 

12 

100 

43 

9 

60 

46 

12 

60 

50 

• 9 

70 

51 

12 

60 

54 

55 

56 

12 

60 

21 

60-75 

59 

60 
63 

12 

60 

12 

60 

64 

36 

60 


Other 


60 hrs. submerged in ^0. 
Submerged in H 2 O bath 
for duration. 


60 hrs. submerged in ^0 
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TABLE I I -A- 5 

WEIGHT GAIN CHARACTERISTICS OF RESIN EMERSED IN WATER 


Sample 

Time 

Temp. 

Weight 

(g) 

A wt. 

A 

Number 

(hrs.) 

(°C) 

Before 

After 

(g) 

(%) 

55 

60 

Room 

5.570 

5.602 

.032 

.57 

56 

60 

Room 

5.5851 

5.7179 

.1328 

2.38 

60 

36 

60-75 

3.9132 

4.0512 

.1380 

3.53 



TABLE II-A-6 


SHRINKAGE DATA FOR TYPICALLY POST-CURED RESIN MATERIALS 


Sample 

Number 

Initial 
Diameter 
(in. ) 

Diameter After 
Post-Curing at 
60°C for 12 Hrs. 

A 

(%) 

Initial 

Thickness 

(in.) 

Thickness After 
Post-Curing at 
60°C for 12 Hrs. 

A 

(%) 

46 

1.513 

1.510 

-.20 

.109 

.109 

0 

51 

1.508 

1.509 

.07 

.105 

.105 

0 

54 

1.512 

1.513 

.07 

.141 

.141 

0 

59 

1.511 

1.510 

-.07 

.109 

.110 

0.92 

63 

1.513 

1.512 

-.07 

.103 

.103 

0 

64 

1.505 

1 . 510 

.33* 

.075 

.076* 

4 ; 

1.33 


* 

Data is for post-curing at 100°C for 36 hours . 
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TABLE I I -A- 7 



PERIODIC SHRINKAGE 

DATA FOR RESIN 

SAMPLE 

NUMBER 

t 

(hrs. at 100°C) 

Diameter 

(in.) 

Thickness 
in. (in) 

0 

1.508 

.228 

(5.800) 

2 

1.516 

.231 

(5.874) 

4 

1.512 

.244 

(6.206) 

6 

1.514 

.246 

(6.264) 

18 

1.513 

.246 

(6.256) 

1 hour after 

1.512 

.246 

(6.250) 



Initial 

Final' 

1 Hr. Cool Off 

A 

A 


(in. ) 

(in.) 

(in. ) 

(in. ) 

(%) 

* 

Diameter 

1.508 

1.513 

1.512 

.004 

.26 

Thickness 

. .228 

.246 

.246 

.018 

7.89 


* 

in inches 
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of graphite prepreg cured in a Dake 75-ton press. The tem- 
perature was increased to 174 °F in approximately 40 minutes 
at which time 300 psi pressure was applied. The temperature 
was then further increased to 250°F in approximately twenty 
minutes. Temperature and pressure were held constant for 
two hours. The part was then cooled to 75°F, the pressure 
removed and the cure process considered complete. Twelve 
layers of bleeders were used to absorb the excess resin and 
to increase the volume fraction of the fibers. The second 
technique was to make composites by using fiber bundles and 
a room temperature-curing resin. The idea was to attempt to 
characterize the two ends of the fabrication spectrum; i.e., 
well-made and poorly-made composites.'' 

Specimen #20 was chosen for investigating the prepreg 
fabrication technique. Three test samples were cut for each 
fiber orientation, longitudinal and transverse. The samples 
were then tested in an Instron 1333 mechanical tester. Spe- 
cifications, test results and calculated values for maximum 
shear stress are presented for the specimens m Table II-A-8 . 

The equation utilized for the calculation is S = 3P/4bd 

max 

F91 

and is given m the ASTM Book of Standards, Part 36 , 

where P = applied load, and b and d are the width and depth 
of the cross section, respectively. 

Results to date are available only for the 250°F, 300 
psi cured prepreg. Sideviews of the failures in the longi- 
tudinal specimens are given in Figure II-A-14. Parts (A) and 
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TABLE II -A- 8 
SHORTBEAM STRESS DATA 


Sample 

Length 
(in. ) 

d 

Thickness 
at Midspan 
(in. ) 

b 

Width 
at Midspan 
(in. ) 

Load 

(lbs.) 

Shear 

Stress 

(psi) 

20 -XT A 

.921 

.129 + 

.003 

-.004 

* 267 +.001 

116 

2526 

20-XTB 

.982 

.129 + 

.001 

-.008 

* +.002 

84 

1789 

20-XTC 

.911 

.131 + 

. 002 

.284 “* 005 

143 

2883 


Sample 

Load 

(lbs.) 

bd 

(in. 2 ) 

Shear 

Stress 

(psi) 

20-XLa 

468 

.0826 

4249 

20 -XL 3 

431 

.0782 

4134 

20-XLy 

456 

.0772 

4430 


U1 

o 



Figure II-A-14. 


Sideviews of Failed Short Beam Shear Specimens 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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and (C) show planar failure surfaces, while part (B) exhibits 
a beveled type of surface failure. A head-on look at the 
failure surface is provided in Figure II-A-15. As can be 
seen, the specimens failed in a brittle manner. High vari- 
ability of results often is associated with brittle fracture, 
and, in these cases, failure stress cannot be described with 
a single number but as a probability of failure at a given 
stress . 

One problem encountered in using a three-point bending 
test is in choosing the appropriate span-to-depth ratio 
which insures that the specimen will fail in shear. If the 
span-to-depth ratio is too large, then the specimen fails 
in bending on the tension side, while a bearing -compressive 
type failure, under the loading point, results if the span- 
to-depth ratio is not high enough. The longitudinal speci- 
men #20 samples are suspected of failing in a compressive 
mode, although some shear cracks are visible. This is sub- 
stantiated by the pictures in Figures II-A-16a and -16b and 
-17 for y, a and B , respectively. This type of microstruc- 
ture is indicative of a compressive failure in a bend test. 
Figure II-A-18 is an exploded view of the B specimen. A 
side experiment that would add insight into the phenomenon 
of failure via compression would be a photographic examina- 
tion of a specimen with a single fiber failed in a buckling 
mode. Future experiments will be conducted with further 
steps taken to insure failure in the desired mode. 


U1 
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Figure II-A-15. 


Photographs of Failed Cross-Sections Short Beam Shear Specimens 





(A) Specimen 20-XL y (lOOX) 



(B) Specimen 20-XL a (lOOX) 


Figure II-A-16. Evidence of Compressive Failure in 
Short Beam Shear Specimen 20-XL 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 
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ORIGIN AC PAGE 

BLACK AND white photograp 



Specimen 20-XL g (lOOX) 


Figure II-A-17 


Evidence of Compressive Failure in 
Short Beam Shear Specimen 20-XL 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 

** 



Reproduced from 
best available copy. 



Figure II-A-18. Exploded Views of Specimen 20-XL 8 (500X) 
(showing compressive failure and laminar 
cracking and misalignment; note the 
moisture spot) 
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The materials used in the specimen molding were Hercu- 
les' HMS fiber strands and the Miller-Stephenson resin sys- 
tem. The V-40 and 828 compositions were varied and are 
given along with specimen and test dimensions in Table II-A-9. 
In many cases post-test examination showed that the interior 
fibers were not completely wet, thus explaining the low 
shear failure stress results. Future work will emphasize 
refining fabrication techniques to prevent this phenomenon 
from reoccurring. 

d. Thermal Expansion Coefficients 

The thermal coefficient of expansion, CTE, is an impor- 
tant quantity in the process of determining composite and 
fiber constituent data. Further, transverse CTE may also 
provide another way to evaluate property variations with 
homogeneity. 

Specimens for these kinds of tests have already been 
fabricated . To date, a number of samples have been 
tested for CTE. The first, a practice run, was tested at a 
rate of 5°C/min. This rate was found to be too fast to a- 
chieve thermal equilibrium in the sample. One means avail- 
able to achieve equilibrium is long thermal stabilization 
times. However, there are advantages to finding a heating 
rate that facilitates a constant heating mode, thereby alle- 
viating the need to reach complete thermal equilibrium. The 
problem then translates into one of compensating for the poor 
thermal conductivity of the material. 


TABLE II-A-9 


SHORTBEAM SHEAR TEST RESULTS OF VARIOUS COMPOSITIONS OF V-40/828* 


b 

Length Width 


Specimen 

Composition 

(in.) 

(in.) 

A 

50% 

V-40 

1.751 

.245 

B 

50% 

V-40 

1.751 

.243 

C 

40% 

V-40 

1.654 

.185 

D 

50% 

V-40 

1.546 

.245 

E 

50% 

V-40 

1.739 

.248 

F 

50% 

V-40 

1.732 

.245 


d 

Thickness 

(in.) 

Span 

(m.) 

Span/Depth 

S ** 

max 

(psi) 

.112 

.750 

6.70 

2400 

.129 

.750 

5.81 

5275 

.095 

.750 

7.89 

3875 

.102 

.750 

7.35 

4200 

.049 

.750 

15.62 

2400 

.082 

.250 

3.05 

5500 


* 

V-40/828 - reinforced with Hercules As-4 Fiber 

** 

S 

max 


= 3P/4bd 
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In order to estimate the thermal transmittance of the 
samples, a sample was placed between two blocks of copper, 
and one of the blocks was placed on a hot plate heating 
source . One would expect that a uniform temperature would 
be obtained in the copper block. Even with the large heat 
sink provided by the second copper block, the temperature 
difference never exceeded 6 3 s°C. The temperature difference 
in a specimen undergoing CTE testing at a heating rate of 
l°C/min. is, therefore, estimated to be less than one-half 
this value for the composites that are to be studied. The 
effect of not reaching equilibrium with a constant heating 
rate seems likely to be just a lag effect on the sample. It 
may, therefore be a reasonable approximation to consider that 
the effective terminus of the testing temperature range is 
shifted 2 to 4°C lower. To test this hypothesis, a sample 
will be run at a rate of l°C/min., and then the same sample 
will be allowed to stabilize and the two results compared to 
determine how close the CTE measured with the l°C/min. con- 
stant rate is to that measured using stepped equilibrium con- 
ditions. 

Specimens for which results have been obtained thus far 
were tested on a Dilatronic brand dilatometer manufactured 
by Theta Industries. The furnace used, a 160 OC, is capable 
of temperatures over 1500°C. The process is regulated and 
controlled by a transducer signal conditioner, and the ap- 
paratus is monitored by built-in programming capabilities. 
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Data is recorded automatically at desired intervals by a 
Fluke data logger, model 2200B. Three critical parameters 
are monitored every two minutes: 1) furnace temperature, 

2) sample temperature and 3) dilatometer millivolt output. 

The change in millivolts can be easily calibrated in terms 
of expansion by a micrometer included in the test apparatus. 
The machine can accommodate specimens of any length up to 
two inches and any width up to one-half inch. Any geometry 
can be handled, but the most accurate results are obtained 
with a flat plate of length approximately one inch and width 
one -ha 1 f inc h . 

Several such samples have already been tested. The 
specimens and their dimensions, weights (before and after 
testing) and drying times at 150°C before testing were re- 
corded. Neglible change in weight indicates that no moisture 
gain or loss occurred during testing. Since the sample 
holder is made of alumina, a fused quartz cylinder was run 
as a standard in order to determine and subtract out the ex- 
pansion of the alumina. The quartz also serves as a check 
on the accuracy of the system and as a means of calibrating 
it or enabling a correction factor to be calculated. The 
expansion of the quartz, the temperature range and the heat- 
ing rate are given in Table II-A-10. This table also pre- 
sents values for some of the samples already tested. The 
curve of the expansion of the quartz must be subtracted from 
(or, as the case may be, added to) the expansion curve of 
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TABLE II -A- 10 

COEFFICIENT OF THERMAL EXPANSION FOR COMPOSITE MATERIALS AND QUARTZ 


CTE (xlO G ) Range 


Specimen 

Lower Temp °C 

Upper Temp °C 

in./xn./°< 

1-3TA 

30.3 

107.9 

65.0 

1-3TB 

28.0 

107.0 

6l.l 

2-3TA 

25.3 

107.0 

64.0 

2 -3TB 

27.8 

107.0 

64 . 6 

3-XTA 

23.4 

107.5 

64.2 

3-XTB 

31.3 

106.9 

62.1 

4-lTA^ 

29-7 

107.7 

62.3 

4-1TA 

28.0 

107.2 

60.3 

4-1TB 

31.3 

106.6 

60.7 

5-XTA* 

27-9 

107.7 

69.9 

5-XTA 

26.3 

107.0 

61.6 

5-XTB 

30.3 

106.8 

62.7 

5-51a 

26.4 

106.4 

-8.4 

6-3TA 

2 6.6 

107.1 

63.1 

6- 3TB 

30.6 

106.9 

58.7 

6-31a 

28.9 

107.9 

-9.4 

8-4TA 

30.8 

107.1 

33.3 

8-5TA 

30.7 

106.5 

40.2 

t8-XT 

22.1 

151.8 

47.3 

8-xla 

2 9-9 

106.8 

-4.9 

9-XTA 

21.7 

107.4 

34.9 

9-XTB 

29.4 

107.4 

34.5 

10-5TA 

29-0 

107.1 

33.1 

10-5TB 

31.8 

107.6 

29.7 

10-51a 

26.6 

107.4 

-6.3 

11-3TA 

30.3 

107.2 

40.5 

11-3TB 

30.1 

106.6 

37.0 

ll-51a 

27.1 

106.5 

-6.7 

12-1TA 

28.7 

107.0 

31.8 

12-1TB 

29.6 

107.7 

31.7 

12-4la 

30.2 

108.0 

-5-7 

1 3-XTA 

23.8 

107.8 

36.9 

13-3TA 

30.5 

106.9 

38.1 

13-51a 

29.1 

106.9 

-5.8 

14-4TA 

29.5 

106.4 

28.1 

i4-4tc 

31.4 

106.8 

26.5 

15-TA 

25-7 

107.0 

28.7 

15-TB 

20.3 

106.9 

22.4 

15-LA 

24.8 

108.1 

-6.2 

16-1TA 

30.8 

107.2 

29.3 

16-1TB 

29.6 

106.6 

26.9 

l6-31a 

27.0 

107.4 

-5.9 

1T-2TC 

21.6 

107.9 

26.2 

17-5TA 

26.7 

107.7 

29-4 

17-31a 

21.3 

106.4 

-5.2 


[Cont * &. ] 



TABLE II-A-1Q (Continued) 


-6 

CTE (xlO ) Range 


Specimen 

Lower Temp °C 

Upper Temp °C 

in./in./°C 

18-4TA 

29.7 

' 10t. 8 

23.3 

18-XTA 

30.5 

107.0 

31.3 

18-lla 

28.1 

107.4 

-6.6 

19-XTA 

28.4 

107.7 

24.8 

19-XTB 

30.5 

106.5 

23.5 

20-4TC 

30.3 

106.9 

24.9 

20-5TA 

17*5 

107.3 

25.4 

20-51a 

28.9 

107.3 

-5.9 

21-1TA 

19.8 

106.6 

30.3 

21-1TC 

29.4 

107.6 

28.1 

21-Ula 

29.7 

107.6 

-5.2 

22-4TA 

23.4 

107.2 

26.3 

22-4TB 

30.4 

107.1 

25.6 

23-5TA 

19.6 

107.4 

24.1 

2 3- 5TB 

31.4 

106.8 

21.1 

24-3TA 

28.0 

108.0 

29.7 

24-3TB 

27.4 

107.8 

27.3 

24-51a 

29.1 

107.7 

-3.0 

25-XTB 

33.4 

106.0 

27.6 

25-XTC 

31.1 

107.9 

27.4 

26-1TA 

29.1 

106.4 

26.3 

26-1TB 

29.6 

107.1 

24.7 

26-4la 

2 9.6 

106.4 

-5.6 

27-XTA 

29.6 

106.8 

24.0 

27-XTB 

30.1 

107.4 

21.5 

27-xla 

27.5 

106.8 

-6.1 

28-XTA 

24.9 

106.9 

40.6 

28-XTB 

29-3 

107.3 

44.9 

29-5TA 

19.6 

106.4 

29.9 

29-5TC 

31.3 

108.3 

24.5 

30-4TA 

28.1 

106.5 

38.0 

30-XT 

. 44.7 

106.0 

50.5 

30-xla 

30.5 

108.0 

-1.5 

31-XTA 

23.7 

106.7 

26.5 

31-XTB 

29.6 

106.5 

25.5 

32-1TA 

28.4 

107.6 

43.9 

32-1TC 

30.6 

106.8 

43.3 

32-lla 

28.9 

107.3 

76.3 

33-4TA 

28.6 

107.6 

30.2 

33-4TC 

31.1 

107.3 

25.1 

33-lla 

23.1 

107.4 

-3.8 

34-itc 

18.5 

106.5 

34.9 

34-1TA 

28.7 

106.3 

36.8 

35-2TA 

26.2 

106.7 

42.7 

35~ 2TB 

30.9 

106.9 

39.0 

36-xtb 

17.4 

106.5 

28.0 

36-xtc 

24.2 

105.5 

24.6 


[Cont ’ d. ] 
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TABLE II-A-10 (Continued) 


_6 

CTE (xlO ) Range 

Specimen Lower Temp °C Upper Temp °C in./ in./°C 


37-4TA 

26.9 

106.4 

49.7 

37-4TC 

27.3 

107.0 

50.0 

38-1TA 

26.7 

IO 6.5 

30.4 

38-lTC 

29.6 

107.2 

26.1 

39-XTA 

20.5 

106.5 

23.6 

39-1TB 

27.4 

106.8 

24.7 

QUARTZ 

35-9 

106.8 

3.6 

quartz 

22.8 

142.1 

3.6 

1908 2* 

28.1 

107.0 

69.9 

1908 2 

25.8 

107.3 

67.7 

1908 3 

25.7 

107.0 

63.8 

979-A1* 

24.6 

106.9 

63.2 

979-A1 

25.9 

107.8 

56.2 

979-B1 

24.6 

106.7 

56.6 


t All samples were run at a programmed heating rate of 1 °C/min. 
Sample 8-XT was run at a rate of 5 °C/min. 

* All samples were postcured unless denoted "by a In this 

ease the sample was tested prior to post curing. [Subsequently 
postcured and tested. ] 

NOTE: Items 1-6 Kevlar 

Items 8-39 — Graphite 

Items 9T9 & 1908 Neat Resin 
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the samples to eliminate the expansion of the alumina, since 
the recorded data reflects the differential expansion between 
the alumina boat and the sample. This can be facilitated by 
using a Bascom-Turner 800 smart recorder. After running a 
significant number of CTE experiments, the next phase will 
be to determine the volume fraction for each sample after 
they have all been mechanically tested for elastic constants 
and - because of its importance to homogeneity aspects - 
short-beam shear. Two things are then possible. First, a 
plot of CTE versus volume fraction can be constructed, and 
the transverse CTE value of the fiber determined by extrapo- 
lating to a volume fraction of 100 percent. Second, the 
test data can be correlated and the concept of homogeneity 
included in property assessment. 

4 . Plans for Upcoming Period 
Future work will continue in several major areas, one 
of which is the effect of unit cell periodicity on homogene- 
ity of the composite. As reported here, the preliminary 
work performed xn this area tends to support the hypothesis 
that if the view space is n (i.e., an integer multxple) unit 
cells, the only variance is machine error and that dispersion 
is at a minimum. Future work wxll attempt to model this 
problem via theoretical calculations and compare it to actual 
experimental results. Once this phenomenon is more fully 
understood, the shape of the curve in Figure II-A-12 can be 
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described and perhaps the previously mentioned "machine 
broadening" can be accurately predicted. 

Work has also started on the concept of fiber dimensions 
and agglomerated fiber dimension distributions and their re- 
lationship to packing. As can be seen from Figure II-A-19, 
such items as distribution of fiber diameter, maxim-urn hori- 
zontal length to the longest dimension and maximum width to 
area are being investigated as potential shape factors. The 
concept may also be applicable to resin rich areas and war- 
rants further experimentation. 

Further work will also proceed in the area of resin 
characterization. Mechanical testing will also continue. 

As mentioned in the report, appropriate steps will be taken 
to make sure all fibers are wet, that the composite is not 
too densely packed and that the necessary span-to-depth ratio 
will be incorporated in SBS specimens so that they fail in 
the desired shear mode. Thermal expansion runs and analysis 
will continue. Eventually, the work will be assembled in 
such a way as to enable the accurate prediction of transverse 
properties. 


5. References 
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6 . Current Publications or Presentations by 
Professor Diefendorf on this Subject 

"Principles and Potential of Composites m High Temperature 
and Ceramic Matrix Composites for Oxidizing Atmosphere 
Applications" 

Published in NMAB - 376, October 1981, 47-96. 

"Ceramic Composites" 

Presented at the Norton Company, Wor Chester, MA, 

October 22, 1981. 

"Composite Materials" 

Presented at the General Electric Company, Saratoga, 

NY, October 28, 1981. 



67 


"The Structure of Carbon Fibers and the Consequences on Per- 
formance" 

Presented at the DuPont Pioneering Research Lab. , 
Wilmington, DE, November 3, 1981. 

"The Strength of Carbon Fibers" 

Presented at Exxon Enterprises, Fountain Inn, South 
Carolina, November 4-5, 1981. 

"Aerospace Applications of Composite Materials" 

Presented at the ACS and Society of Plastics Engi- 
meers Meeting, Albany, NY, January 6, 1982. 

"Structure of Carbon Fibers and Consequences on Performance" 

Presented at the Aerospace Corporation, El Segundo, 

CA, January 14, 1982. 

"Residual Stresses in Carbon Fibers" 

Presented at the Gordon Conference on Composite 
Materials, Ventura, CA, January 18-22, 1982. 
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III— A ADVANCED ANALYSIS METHODS 

Senior Investigator: E. J. Brunelle 

1. Introduction 

The ultimate purpose of this research effort is to 
understand exactly how each parameter changes the solution 
of any composite plate problem. Standard methods provide 
answers to specific analysis problems, but they do not ad- 
dress the above statement of intent. 

Recognizing that a primary difficulty is that of a mul- 
tiplicity of variables, solutions for orthotropic plates 
have been sought xn terms of generalized variables. These 
efforts have produced many solutions for individual lamina 
and they are expected to form a basis for solving general 
composite plate problems. The first such solutions (by the 
present investigator) for a composite plate are one of the 
topics reported on in this progress report. 

2 . Status 

The principal results of the last period are summarized 
as follows: 

a) The generally orthotropic plate equation (which 
has the same form as the symmetric angle-ply 
plate) was found to depend on only two para- 
meters in its affine plane representation. 

Thus, the coefficients of the operator 


PRECEDING PAGE BLANK NOT FILMED 
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ORIGINAL p/Ufe Vj } 
f OF POOR QUALITY 

+ ;* 2p jy + r .. + 

9x J 3x o 3 *o 3x o 9y o 3x i 3y o K 

depend only on the values of D 11 / D 2 2 an ^ D * -*- n 
principal directions. Furthermore, L(D*,D i;l / 
d 22 /S) = R(D* / d 11 / d 22 ' 90 ” 0 ^ and vxce versa. 

b) A very efficient numerical procedure was found 
for factoring the quartic associated with the 
above operator. 

c) Exact solution bounds were found for the spe- 
cially orthotropic plate. 

d) The incremental behavior of D* with respect 
to the plate coefficients was delineated. 

e) The similarity rules for free vibration of a 
specially orthotropic plate were presented. 

3. Progress During Report Period 

Progress in the last reporting period has been made 
principally m the three areas described below [(a) - (c) 3 . 


a) The existence of double frequencies for other than 

square (a/b=l) isotropic plates has been questioned and 

hotly debated in the literature. Assuming a simply-sup- 

ported plate, an explicit relation is now presented for 

a /b values at which ft Eft with a parametric depend- 
o o mn m*n* ^ * 

ence on D*. 



D* [ (m*n*) 2 - (mn) 2 ] t / [D* 2 - 1) [ (m*n*) 2 - (mn) 2 ] 2 + [ (m*n) 2 - (mn*) 2 ]’ 2 


n 4 - n* 4 


2 
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Oj’ ; i. < L ; u ; 

where the positive square root is chosen. The values of 

{a /b } 2 for the limiting values of D* = 1 and 0 are: 
o o 


m 


*2 _ 


m 


D* = 1 


n 2 - n* 2 



0 



and, finally, if (m*n*) 2 - (mn) 2 = 0, one determines values 

ct 

of _o that are independent of D* , 


a 

2 

f "I 

o 


m 

b 


n* 

°j 

INDEP. OF D* 



From the expression 



fb ] 

4 

fb 

m 4 

o 

+ n 4 + 2D* (mn) 2 

o 


cl 

l oj 

a 

o 

/ 


it is seen that the frequencies corresponding to the above 
special values of a Q /b o are not independent of the value of 
D*. 


b) During the last dozen years many solutions have been 
presented for a wide variety of laminated rectangular compos- 
ite plate problems. Unfortunately, all of these solutions 
are specific (a specific material and a specific aspect 
ratio must be given in order to present numerical results) 
rather than generic. Thus, the general behavior of these 
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solutions as various parameters change is not well under- 
stood. Employing techniques used by the present investigator 
on orthotropic plate problems, both the static deformation 
problem and the eigenvalue problems (vibration and buckling) 
are solved for S2 antisymmetric cross-ply plates with all 
sides simply supported. These solutions are m terms of the 
similarity variable T = which presents the results in a 
highly understandable and intuitively appealing form. 

For example, the frequency spectrum is represented by 

0 2 (B -, ) 2 

= f. (T, D*) - F (T,A,A*)-=— 

n 4 1 1 A ll D ll 

The function f^(T/D*) is the orthotropic plate solution so 
that the function F^(T,A,A*) determines the coefficient of 
the bending-extension parameter (B^) 2 /A D . In expanded 
form. 


F (T,A,A*) 


2A*T~ 4 + (AT~ 2 + 1) T~ 3 + A + T~ 2 
(A + T" 2 ) (1 + AT" 2 ) - (A*?" 1 ) 2 


A A 66 //A ll 


A i ? + A fifi A i 9 

A* = -A!- ii = A + v— — • 

A 11 A 11 


Note that the bending extension matrix element appears 

only in the combination (B.^) 2 / A i 1 d 1 i shown above, i.e., it 
does not appear in the F^ (or f^) function. 
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c) Many advanced composites have large values of Q.^ and 

Q 22 such that Q 11 rQ 2 2 >> ^ 12 '^66* 0ne conse( 3 uence of this 

is that the parameter D* is in the .1 to .2 range. It has 

been demonstrated many times (exclusively at RPI) that the 

"D* = 0 solution" is always in the neighborhood of low D* 

(.1 to .2) solutions. Therefore, the "D* = 0 solution" 

(which also really means that the ratios of A^/VA, ,A„ and 

66 11 2 2 

^ A 12 + A 6 6^ //i/A 11 A 2 2 are neglected) is a very convenient and 

• 4 * 

very useful extremum solution. For example, when consider- 
ing S2 antisymmetric cross-ply plates, the "D* = 0 solution" 
is physically interpreted as neglecting the Poisson coupling 
and the twist-curvature coupling but retaining fully the 
bending-extension composite plate couplings. Thus, the three 
u, v, and w coupled equations of motion take on the simple 
form 


B 


u 


11 


r ~ „ - *W, 
xx A i;l yyy 

B, 


V, 


11 


- * ~ w ' 
yy Ajj_ yyy 


D (w, + w, ) - B (u, - v, ) + Pw, + . . 

11 xxxx yyyy 11 xxx yyy ww 


■ + P w / tt = <i( x ^yrt) 


It is easily verified that the uncoupled w equation is given 
as 


w. 


xxxx 


+ w. 


yyyy (I-j)d w, xx ' (i - j)d . "'tt (l- j)d 


w,. 


„ q(x,y,t) 


11 


11 


'll 


^Upper bounds for deflection problems and lower bounds for 
eigenvalue problems 
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where J = ( b i;l ) 2 / a 11 d 11 * This equation yields extremum solu- 
tions for the cross-ply antisymmetric plate. Notice that xt 
is a separable equation, so that arbitrary boundary conditions 
may be imposed on all four sides. For example, defxning the 
buckling coefficient at k Q = Pb 2 /iT 2 (1 - J)D i;l the k Q versus 
a o /b Q buckling curve for all sides clamped is shown in Fig- 
ure III-A-1. This is an exact solution, which has as one of 
its most important features the fact that J (the effect of 
the bending-extension coupling) appears only in the denomi- 
nator of the k defxnition. 
x 

4 . Plans for Upcoming Period 

Research on this project will be suspended temporarily 

this summer, but will resume with the beginning of the Fall 

*82 Semester. At that time the work on similarity solutions 

will be extended to S3 angle-ply antisymmetric plates, and 

U 2 

earlxer work on the materxal parameter B = will be 
suitably generalxzed and presented in what is hoped to be a 
final form. 


5. Current Publications or Presentations by 
Professor Brunelle on this Subject 


"Thermal Buckling of Initially Stressed Thick Plates" 
To be published in the ASME Journal. 


"Generic Buckling Curves for Specially Orthotropic Rectangular 
Plates", with G. Oyibo. 

Accepted for publication by the AIAA Journal, June 
4, 1982. 
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0 .7 1.0 1.5 2.0 

Figure III-A-1. Critical Buckling Coefficient vs. Plate Aspect Ratio 
(all sides clamped) 
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"Eigenvalue Similarity Rules for a Class of Rectangular Spe- 
cially Orthotropic Laminated Plates" 

To be presented at the Ninth U. S. National Congress 
of Applied Mechanics, Ithaca, NY, June 21-25, 1982. 

"Similarity Solutions for the Eigenvalue Problems of S2 Anti- 
symmetric Cross-Ply Composite Plates" 

To be presented at the SIAM 30th Anniversary Meeting, 
Stanford, CA, July 19-23, 1982. 
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III-B FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator: E. Krempl 

1. Introduction 

The deformation and failure behavior of graphite-epoxy 
tubes under biaxial (tension, torsion) loading is investi- 
gated with the assistance of Ph.D. candidate T.-M. Niu. The 
aim of this research is to provide basic understanding and 
design information on the biaxial response of graphite-epoxy 
composites. 

2. Status 

Using graphite-epoxy [±45]^ thin-walled tubes, the 
static axial and torsional strength were determined together 
with the elastic moduli. Fatigue tests under completely re- 
versed, load-controlled axial loading showed comparatively 
poor axial fatigue strength; furthermore, fatigue life was 

affected by frequency. A paper on this subject will appear 

I’ll * 

m the Journal of Composite Materials . 

3 . Progress During Report Period 
During the report period, static combined loading tests 
as well as combined tension-torsion, completely reversed, 
load-controlled fatigue tests were performed with the loads 
in-phase and out-of-phase. A total of 58 tubes were tested. 

_ 

Numbers in brackets in this section refer to the references 
which are listed on page 85. 
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Publications are in. preparation and only highlights are 
given here. 

a. Combined Static Tests 

Test results are shown in Figure III-B-1 where cr and x 
represent the axial and shear stress components, respec- 
tively. Previously only pure axial and torsion tests were 
reported ^ . It is seen that the combined tests complement 
the previously reported data very well. 

rn 

It was noted m our earlier publication that the low 
negative torsional strength of the tubes is due to local 
buckling. The combined loading results suggest that local 
buckling reduces the strength in cases where negative torque 
{the fibers of the outer layer are tensed) is applied com- 
pared to the cases where positive torque (the fibers of the 
outer layer are compressed) is involved. The failure enve- 
lope can be represented by two half ellipses with the same 

* 

"axial half axis" as shown. One represents the local buck- 
ling failure mode; the other can be thought of as represent- 
ing fiber breakage. 

b. Combined Fatigue Tests 

All tests were run m the completely reversed, load- 
controlled mode using a 1 Hz sine wave. If the variations 
in applied axial and shear stress are in phase with each 
other (<j> = 0) , the image of this loading in the o - x 
plajie is a straight line which passes through the first 
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CD— STRESSES AT FAILURE 


Figure III B-l. Biaxial Failure Locus for Static Tests 

(Two failure mechanisms - buckling and fiber breakage - 
are operative and result m two distinct surfaces.) 
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and third quadrant (see Figure III-B-2),. For a phase shift 
4> = it (tf) denotes the phase shift) between the axial and the 
shear stress, a straight line in the second and fourth quad- 
rant is obtained. For phase shifts of tt/2 and 3ir/2, an el- 
lipse is obtained in the a - t plane, with clockwise and 
counterclockwise rotation for tt/2 and 3ir/2, respectively. 

Since <j> = 0 and <J> = ir result in straight line images on 
the a - x plane, they will be referred to as in-phase load- 
ings. 

The ratio between the axial and shear stress amplitude 
is characterized by the angle 0 such that 

x 

tan 0 = — — 
a 

a 

where the subscript 'a' designates the amplitude. Values of 
0 > tt/2 designate the tests with images in the second and 
fourth quadrants. 

Figures III-B-3 and -4 show the results of the combined 
load fatigue testing. In this figure relative axial stress 
amplitude is plotted versus cycles to failure. Figure III-B 
-4 shows the results on the basis of the relative shear 
stress amplitude. (The data are normalized by o = 156 MPa 
and by = 120 MPa.) 

As expected, the pure torsional fatigue performance is 
very good. No fatigue failure was obtained at 95% stress 
amplitude in 10 6 cycles (the failure recorded at 10 5 cycles 
was obtained on a specimen with an initial flaw; an arrow 
next to a symbol designates a run-out) . 
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Figure III-B-3. 


Fatigue Test Results Plotting a /a versus Cycles 
to~Failure (a = 156 MPa) max U 




85 


Although the uniaxial case - where deformation is strong- 
ly matrix dominated - was expected to be the worst case in 
fatigue performance, the data in Figures III-B-3 and -4 show 
that the biaxial loading case with 0 = 49.1 degrees exhibits 
the worst fatigue performance of all the combinations tested. 

From the static test results shown in Figure III-B-1, 
two failure mechanisms are apparent. Buckling was shown to 
be limiting for negative twist. It appears that buckling is 
also limiting in the fatigue case. At amplitudes just below 
the torsional buckling strength, no fatigue failure is ob- 
served in pure torsion. This behavior is made evident in 
Figure III-B-5 which shows a cross-plot of the {sometimes 
extrapolated) data in Figures III-B-3 and -4. (Only in- 
phase test results can be plotted in Figure III-B-5.) 

All points are inside the small ellipse which represents the 
buckling mode of failure in static loading. 

4 . Plans for Upcoming Period 
The fatigue test results will be analyzed and manuscripts 
for publication will be prepared. 

5. References 

1. Krempl, E. and T.-M. Niu, Journal of Composite Materials, 
Vol . 16, 1982, pp. 172-187. 
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Figure III-B-4. Same Test Results as in Figure III-B-3 

fr _ A is used as ordinate T = 120 MPa) 
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Figure III-B-5. 


Surfaces of Constant Fatigue Lives 
(obtained using Figures III-B-3 and -4) 
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6. Current Publications or Presentations by 
Professor Krempl on this Subject 


"Graphite-Epoxy [±45] s Tubes. Their Static Axial and Shear 
Properties and Their Fatigue Behavior under Completely Re- 
versed Load Controlled Loading" , with T.-M. Niu. 

Published in the Journal of Composite Materials, 

Vol . 16, 1982, 172-187. 

"Biaxial In-Phase and Out-of -Phase Behavior of Graphite- 
Epoxy Tubes" 

To be presented at the ASTM Conference on Biaxial/ 
Multiaxial Fatigue, San Francisco, CA, December 
15-17, 1982. 
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III-C MOISTURE AND TEMPERATURE EFFECTS ON THE MECHANICAL 
PROPERTIES OF LAMINATES 

Senior Investigator: S. Sternstein 

1 . Introduction 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
physical properties of the matrix resin. To date, the re- 
search has involved the precise viscoelastic characterization 
of epoxy neat resins, interlaminar failure of composites and 
the inhomogeneous swelling of and the effects of moisture on 
composites. 


2 . Status 

It has been shown that the dispersion characteristics 
of neat resin and in situ resin are equivalent and that time- 
temperature superposition is obeyed in the glass transition 
region of the epoxy. Using the centro-symmetric deformation 
(CSD) testing method, it was found that when a postcured 
laminate sample is subjected to moisture its glass transi- 
tion temperature varies directly with moisture content as 
measured by weight gam. Further, that a sample which is 
dried will recover its original glass transition temperature. 
To the contrary, laminate samples which have not been pre- 
viously post-cured, when initially exposed to moisture, ex- 
perience a decrease of glass transition temperature which is 
not recovered after drying. Typically, the decrease in T 

y 

is five to ten degrees, depending on the moisture treatment. 
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Mechanical tests to determine dynamic modulus variations 
versus time, at constant temperature on wet samples, indicate 
that moisture interacts strongly with the postcure kinetics. 
It appears that even after the sample has dried there is 
still a residual effect on postcure kinetics. 

A broad out-of-phase modulus peak in the dynamic be- 
havior of wet neat resin samples of Narmco 5208, extending 
from 160°C to 240°C, contrasted with a 20°C to 30°C wide 
peak for a dry sample. This broadness of the wet glass 
transition peak may be direct rheological evidence for the 
presence of structure in the resin. Previously, other in- 
vestigators have suggested that an amorphous structure ex- 
ists in thermosets, based on electron microscopy of fracture 
surfaces. In the work reported here, the results represent 
bulk behavior, not surface behavior. These results may be 
of considerable utility in probing the structural features of 
epoxy resins. 


3. Progress During Report Period 
The characterization of matrix dominated mechanical 
properties of carbon-epoxy laminates using dynamic mechanical 
techniques is being continued . The sample test geometry in- 
volving centro- symmetric deformation (CSD) of a circular 
disc has been described previously and some of the early re- 
sults with this technique reported. We are now using dynamic 
techniques to probe the damage induced by water absorption 
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and desorption. The CSD test geometry was used in earlier 
project periods to measure the delamination strengths of 
laminates of various stacking sequences and the results re- 
ported. In that study, it was found that the delamination 
strength and load drop at the first failure event were both 
markedly altered by the presence of moisture uptake-produced, 
delamination cracking, even prior to sample loading. 

As a result of this previous study, we felt it necessary 
to examine in some detail the development of damage in lami- 
nates as a result of the moisture uptake process itself and 
without the influence of externally applied loads on the 
sample. Dynamic mechanical techniques appear to be a good 
method for characterizing the damage, per se, since this 
method, utilizes sample loads which are about a factor of 30 
lower than the first delamination event load level, and con- 
sequently, is not likely by itself to promote additional 
damage. 

In order to isolate the damage caused by the uptake of 
water from the possible additional damage caused by subse- 
quent drying of the sample, we felt that it would be best to 
obtain the dynamic data at 30°C and not at elevated tempera- 
tures. Thus, we would be able to control the drying condi- 
tions in a separate experiment. 

The question which remains to be answered is whether or 
not dynamic mechanical measurements are reflective of the 
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damage process and if so, how sensitive are they? The pre- 
liminary results are most encouraging. 

The storage (m-phase) and loss (out-of -phase) modulii 
of a 12-ply carbon-epoxy -laminate which has been kept dry 
are shown versus frequency for 30 °C in Figure III-C-1. After 
boiling the sample for 210 hours in a water bath, the results 
(again at 30°C) are shown in Figure III-C-2. There is a 
small but appreciable decrease in the storage modulus as a 
result of the water treatment. The loss modulus, however, 
is increased significantly (a factor of two at 0.1 Hz) . The 
sample was then partially dried until its residual water 
weight gam reached one percent and the data m Figure III-C-3 
were obtained. Note that the storage modulus is significant- 
ly lower than in Figure III-C-2, even though the sample has 
been partially dried. Finally, the sample was dried to a 
residual weight gain of 0.25% and the results in Figure III- 
C-4 were obtained. There is a slight increase in the storage 
modulus relative to Figure III-C-3 but the loss modulus re- 
mains high (compare Figures III-C-1 and -4) . 

The storage modulii results are summarized in Figure 
III-C-5. Note that the values shown (actually stiffnesses) 
have been normalized to the thickness of the dry sample by 
correcting for the thickness-cubed effect of out-of-plane 
deformation. Clearly, the effect of moisture absorption 
followed by drying results m a permanent reduction of flex- 
ural stiffness. 
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Figure III-C-1. Storage and Loss Moduli of a Dry Carbon/Epoxy Laminate 
(30°c) 
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Figure XII-C-2. Storage and Loss Moduli of a Wgfc Carbon/Epoxy Laminate 

{30°C; after boiling the Figure III-C-1 specimen in water) 
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Figure III-C-3. Storage and Loss Moduli of a Partially Dried Carbon/Epoxy Laminate 
(30°C; after drying the Figure III-C-2 specimen to 1% weight gam) 
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The loss moduli! data are summarized in Figure III-C-6 
where the loss factor (tan 6 , the ratio of M'VM') relative 
to the dry state is shown versus moisture level. Several 
experiments and samples are summarized on thus graph. If 
moisture level is increased only (no drying) , all samples 
indicate a linear increase of tan 6 with moisture level. If 
the sample is dried beginning with any level of moisture, 
then a permanent residual increase in tan <5 is observed. 

What is somewhat suprising is that even the sample which has 
only been subjected to a weight gain of ca. 1% shows a hys- 
teresis effect on drying 

The results of Figure III-C-6 strongly suggest that dy- 
namic mechanical loss factor can be used as a measure of 
moisture induced residual damage in carbon-epoxy laminates. 
This work is being continued. 


4. Current Publications or Presentations by 
Professor Sternstein on this Subject 


"Matrix Dominated Mechanical Properties of Composites" 

Presented at an S.P.E. Short Course, Louisville, 

KY, November 4, 1981. 

" Composites" 

Presented at the CIBA-Geigy Research Center, Tarry- 
town, NY, November 19, 1981. 

"Viscoelastic Characterization of Composites" 

Presented at the Naval Research Labs., December 1, 
1981. 

"Matrix Dominated Mechanical Properties of High Performance 
Composites" 

Presented at the International Conference on Defor- 
mation, Yield and Fracture of Polymers, Churchill 
College, Cambridge, England, March 29, 1982. 
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Preceding page blank 


III-D NUMERICAL INVESTIGATION OF MOISTURE EFFECTS 

Senior Investigator: M. S. Shephard 

1 . Introduction 

This analysis builds on the theoretical developments 
of Professor S. Sternstein (reported in progress reports dated 
December '79 and June ’80) in which the fully-coupled, non- 
linear thermomechanical equations of inhomogeneous swelling 
of composites m the presence of moisture and temperature 
are analyzed. Here the one-dimensional case of a single 
fiber in an infinite matrix is being extended to two-dimen- 
sional, multiple fiber cases. 

2 . Status 

As discussed in earlier progress reports, moisture ef- 
fects were introduced into the problem through a nonlinear 
constitutive relation, the principle of virtual displace- 
ments was employed to develop a nonlinear matrix equation 
using displacement-based finite elements and these nonlinear 
equations were solved using quasi-Newton methods. Subse- 
quently, the constitutive equation was built into the finite 
element analysis. The results of a single fiber comparison 
analysis were also given. 

3. Progress During Report Period 

Graduate student Frida Lumban-Tobing has continued, 
during the last reporting period, testing the program with 
sample problems and establishing additional features in it. 
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The areas of major analysis improvement are in the quasi- 

Newton solution. Improvements have also been made in the 

two-dimensional postprocessor, allowing the user to plot 

additional information. These additions include contours 

of the liquid volume fraction, principle strains, strain 

components and additional stress components. The following 

* 

two sections, extracted from Reference [1] , discuss the 
numerical solution procedure and the results of a two-fiber 
problem. 

a. Numerical Solution 

a.l) Quasi-Newton Method . Newton's method could be used 
to solve the nonlinear algebraic equations resulting from the 
application of the principal of virtual work. For a large 
system of equations, however, this method is expensive, be- 
cause the Jacobian matrix has to be inverted at every itera- 
tion. To avoid this, the quasi-Newton method was applied to 

update the inverse of the stiffness matrix at every iteration. 

F2 31 

Of the several quasi-Newton methods available ' J , the Da- 
vidon update, which is a rank one correction, was selected 
here for two major reasons. First, this is a symmetric up- 
date, i.e., the updated stiffness matrix remains symmetric. 
Second, it is less costly per iteration than rank two cor- 
rection methods and appears to be less expensive overall. 


Numbers in brackets in this section refer to the references 
which are listed on page 118 . 
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In the quasi-Newton method, the Jacobian matrix is ap- 
proximated by evaluating the residual at two successive 
points, that is 

G ( A) (A - A) = r (A) - r ( A) 
or 

G(A)d = y - (1} 

\ 

where 

G(A) is an approximation of the Jacobian matrix 

A and A are the column matrices of displacements 
for two steps 

r(A) and r(A) are the column matrices of residuals 
for the same two steps 

d = A - A 

y. = r (A) - r (A) 

til 

For the Davidon update, the matrix G for the (k + 1) 
iteration is updated in the following form 


G, . -» — G, + 
..k + 1 _k 


( ?k ~ ( yk - 

( ?k - ?k?k )T ?k 


T 


( 2 ) 


The superscript T indicates the transpose of the matrix. To 

avoid inverting the matrix G at every iteration, it must be 

possible to directly update the inverse of the matrix G for 
til 

the (k+1) iteration. This can be done by inverting equa- 
tion (2) using Householder's formula 

(A + gcb ^) -1 - A " 1 - yxz T 
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where 


q = a scalar, 

b and c are column matrices 

and 

x = A 1 c 
z = A" X b 

Y = q(l + qb A c) 

' A> >V 

which for equation (2) gives 


G 


-1 

k + 1 




~k ( ~k ~ ~k ~k^ 

( 5k ' 


(3) 


For a given set of initial values A n and G , the inverse up- 

~ u ~ u 

date m equation (3) can be written as 


?k + 1 



k 

+ E 
i=0 


T 

s.sT/3. 

-i-i 7 l 


(4. a) 


where 


s . 
~1 

-l 

= d - G. y. 

~i ~i ~i 

(4 . b) 

p i 

T 

= s . y . 

~i~i 

(4.c) 

2i 

= r . , _ - r . 

~i + 1 ~i 

(4.d) 

d. 

~i 

= t.d 

i~i 

(4 . e) 

d. 

~i 

= -GT 1 r. 

~i ~i 

(4.f) 

r . 

~i 

= r(A ) 



and t^ is a scalar obtained using a line search which will be 
discussed in the next section. 
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a. 2) Devices to Insure and Speed Convergence . To in- 
sure and speed up convergence, it is important that the ap- 
proximate Jacobian be as accurate as possible. In this 
study, two devices were used to insure that. One is the use of 
a line search to select reasonable step sizes. The other is 
a domain check which will force any unacceptable solution 
step back into the feasible domain. 

a. 2.1) Domain of the problem: All of Rivlin's 

strain invariants have to be positive because 
they are squares of functions of the extension 
ratios. In addition, the Flory-Huggins equation 
requires the volume fraction of polymer in the 
swollen body, v 2 , to be less than 1. This re- 
striction can be stated as 

v 2 = V*/VT; < 1 (5. a) 

or 

l 3 > (v*) 2 (5.b) 

Thus, the following conditions must be met to in- 
sure that the domain of the problem is not vio- 
lated 

X 1 > 0 

1 2 > 0 (6) 

1 3 > (V*) 2 

If at any point during the analysis, these domain 
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requirements are not met, the line search tech- 
nique is used to bring the analysis back inside 
the domain. 

[4 51 

a. 2. 2) Line search : Bathe and Matthies ’ sug- 

gested the use of line search techniques to improve 
convergence of the quasi-Newton method. The sug- 
gested line search requires the use of step length 
t, in the direction of d. such that 

X K 

I Sk£ ( 5k + W I i 11 I 5kE ( 5k> I 

where n is convergence tolerance. They suggested 

[31 

the value of r\ to be 0.5. Geradxn et al. con- 
cluded that this line search did not give substan- 
tial improvement to the convergence rate, and 
[4] 

Bathe indicated that , for some problems , a 
still larger value of n yields good convergence 
rates. 

For the problem reported here, the line 
search is needed in some steps to insure that the 
solution remains in the domain of the problem. 
Comparing the convergence rate for values of 
n = 0.5 and n - 0 . 9 indicated that the convergence 
rate for n = 0.5 is not better than the conver- 
gence rate for n = 0.9. Since the line search is 
an expensive operation it is only conducted while 

I ?kE ( *k + W I > °- 9 I 5k€ ( i k > I 
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Ocasionally, the value of t^. is negative. This 

[6 71 

is not surprising since Bard and Broyden L ' 1 men- 
tioned that negative steps are sometimes neces- 
sary. 


a. 3) Solution Algorithm . For small problems , where the 

number of independent variables is reasonably low, the 
T 

s.s./3- term from equation (4. a) can be stored in a matrix 

~ X~ X X 

form. In this case, the summation process of the updating 
T T 

matrix s.s./S- is simple. The new s.s./$* can be added into 

~i~i l ^ -i-i 7 l 

the matrix that is obtained in the previous iteration. 

T 

For large problems, it is .not convenient to store 

in matrix form. In this case, the p.'s and the vectors s.'s 

i ~i 

must be stored , and the summation process has to be done at 

every iteration. The algorithm shown below indicates how to 

avoid this summation every time G 1 is encountered at each 

T 

iteration. In this way, only one summation of s^s^/3^ per 
iteration needs to be performed. 

Step 1: read A 

u 


Step 2: k = 0 


Ho “ H ( V 


% = ? ( £o> 

Ho “ ?o Ho 


So - -So 
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ORIGINAL, rr-™ uj 
OF POOH QUAUfY 


Step 3 : k = k + 1 


tk-i = 1 


Step 4 


~k ~k - 1 + fc k - l~k - l 

(find t, , using line search if the solu- 
tion is outside the domain or if 

I?k-i5 ( ?k-i + VlSk-i’ I 1 > °- 9) 

£k = ;<5k) 

check convergence: | | r^ | ] a specified 

value. If satisfied, A is solved, stop. 
Otherwise, continue to Step 5. 


Step 5= y k = - £ k -l 


5k So Hk 


k-1 


°k = < i ^ 0 5i?i /B k ) 5k 


Sk = -5k - 5k 


5k (t k - 1 - l, 5k - 1 + Sk 


B k = 5k?k 


store s^ and 8^ 



(s* s. /6, 
~k~ k k 



Step 6 : go back to step 3 . 

In this algorithm, the only matrix to be inverted is G Q . By 
factoring G into upper and lower triangular matrices, the 

•V U 
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forward reduction and backward substitution can be done 
later when the vector b needs to be calculated. 

For the problem under consideration, two methods can be 
used to start the iteration. The first method is to guess a 
starting v 2 value, apply this to the entire region and use 
it to estimate A Q and G Q . This does not work well for a 
problem with a large range of v 2 values, which is usually 
the case when high tensile load is applied. The second 
method is to use the displacements due to zero load or lower 
load level as the initial displacements. 

b. Example Problems 

b.l) A Pair of Fibers in an Epoxy Matrix Without Ap - 
plied Load . In this example, the interaction between the 
two fibers at a given humidity is studied by varying the 
distance between those fibers. The basic model of this prob- 
lem is shown in Figure III-D-1. It is analyzed using six- 
nodal isoparametric elements. Figure III-D-2 summarizes the 
results for the case of 100% humidity. The values used in 
this figure are for the maximum value anywhere on the inter- 
face, for point A on the interface (see Figure III-D-1) and 
at the center of the model (point B in Figure III-D-1) . Fig- 
ure III-D-2 indicates a strong dependence of fiber spacing 
with the peak values occurring in the one radius case. For 
spacings of less than five radii, the interaction between 
the fibers is significant and thus differs appreciably from 
the case of a single fiber in an infinite matrix. The results 
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also indicate that at reasonable fiber distances the magni- 
tude of the principal stress a 2 is larger than c^. As the 
distance becomes smaller, the interaction becomes stronger, 
and at a certain distance the magnitude of the principal 
stress a 2 becomes less than a 1 . The maximum value of the 
liquid volume fraction increases with decreasing fiber spac- 
ing and at a certain distance it becomes larger than the 
saturated value for epoxy under no load. 

The interaction between the fibers is also affected by 
the level of humidity. Figure III-D-3 shows this interac- 
tion as a function of humidity for the case where the dis- 
tance between the fibers is one fiber radius. This figure 
shows that the interaction becomes weaker as the level of 
humidity is reduced. 

b.2) A Pair of Fibers in an Epoxy Matrix Under Applied 
Load . Here, three groups of problems are presented. The 
particular groupings are selected because the results are 
similar for all the cases within a group. The first group 
consists of cases where the applied load is either uniform 
tension in the x-direction, uniform compression in the y- 
direction, or both tension in x and compression in y-direc- 
tions. The second group consists of cases where the applied 
uniform load is either compression in the x-direction, ten- 
sion in the y-direction, or both compression in x and tension 
in y-directions. The cases in the third group include those 
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where the applied load is one of the following: either 

in both x and y-directions , tension in both x and y-direc- 
tions or the no-load condition. 

The shapes of the contours of the principal stresses 
and volume fraction of liquid in the swollen body wi-thiri ' ' . 
each group have similar shapes but with different values as- 
sociated with those contours. Figures III-D-4 to 9 show the 
shapes of the and c? 2 contours for all three groups under 
100% humidity, where the fiber spacing is one radius. More 
specifically. Figures III-D-4 and -5 are and a contours 
when load is applied in the positive x-direction; Figures 
III-D-6 and -7 are a 1 and o 2 when the load is applied in the 
positive y-direction; and Figures III-D-8 and -9 are a 1 and 
a 2 when the load is applied in both x and y-directions. As 
might have been expected, these results show that the dif- 
ference between the cases in each group cannot be eliminated 
by simple addition or subtraction of one case with another, 
as can be done m linear problems where the superposition 
principle holds. 

Figure ITI-D-10 shows the maximum of a (anywhere on the 
interface) versus applied uniform load, for the problem with 
100% humidity and with fiber spacing of one fiber radius. 

The trends in the maximum value of the liquid volume fraction 
are similar to the trends in the maximum value of the princi- 
pal stress a The maximum values of the stress and the li- 
quid volume fraction always occur at a point on the interface 
which depends on the applied load. 
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b.3) Discussion of Results . The first example indi- 
cates that the stress level is affected by the distance be- 
tween the fibers. When this distance is less than five fiber 
radii, the combination cannot be approximated by a single 
fiber in an infinite medium. As the distance becomes smaller, 
the stress magnitude becomes larger and at a certain distance 
cr becomes larger than the magnitude of o 2 . The liquid vol- 
ume fraction also increases with the decrease in distance 
between the fibers. The lower the humidity level the less 
important is fiber interaction. 

When there is applied load, the effects of nonlinearity 
become more pronounced and the superposition principle can- 
not be used. The trends in the maximum values of the liquid 
volume fraction and those of cf^ are found to be similar. 

In all two-fiber analyses, the stresses developed are 
quite high for cases with both high humidity and small fiber 
spacing, even when there is no applied load. It seems from 
these studies that permanent damage at the interface between 
fiber and epoxy could occur for combinations of these para- 
meters as they are found in practice. 

4. Plans for Upcoming Period 
During the next reporting period, effort -will be ex- 
pended on additional analyses to gain more insight into load- 
humidity interactions. A different form of the constitutive 
relation may also be investigated. 
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III-E NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF 
COMPOSITE FRACTURE 

Senior Investigator: M. S. Shephard 

1 . Introduction 

To understand the mechanism's of failure in composites 
it is necessary to develop insight into the micromechanical 
behavior, including interactions between matrix and fibers 
as the load is increased from zero to that corresponding to 
failure. Investigations of these phenomena, either experi- 
mental or numerical, are difficult. The purpose of this 
project, being carried out by graduate student Nabil Yehia, 
is to develop a nonlinear finite element analysis capability 
for performing numerical investigations of the micromechan- 
ical failure of composites. 

2 . Status 

A survey of the applicable literature has been conducted 
the conclusion drawn that crack modeling ability is not avail 
able in existing programs and a decision was made to design 
and develop a special nonlinear code. The basic structure 
of the required FEM program has been designed and many of 
the procedures to be included started. Coding of the first 
set of features was begun. 

3. ‘Progress During Report Period 
Since the last progress report, effort has been concen- 
trated on designing and implementing the program data base 



120 


and initial program features. Work has also been done on 
the process of piecing together a consistent set of numerical 
fracture mechanics procedures to track the entire fracture 
problem from crack initiation through crack propagation to 
identifying the point of unstable crack growth. 

The program data structure has been designed and the 
linear elastic portions of the program implemented and 
tested. The major components of the linear elastic static 
portions are: 

a. Six different 2-D element types and a vari- 
able-node, 3-D brick element. 

b. An equation solver based on the frontal solu- 
tion algorithm for single and multiple load 
cases . 

c. A generalized stress recovery algorithm. 

Geometric nonlinearities in terms of Green's strain 

tensor and nonlinear elastic material properties based on 
the polynominal strain energy function have been implemented 
but not yet checked. Portions of the material yield cri- 
teria have been coded and are currently being tested. 

The major challenges of the cracking/fracture problem 
include the predicting of crack initiation, the direction of 
crack propagation, the rate of crack growth, debonding at 
the matrix-fiber interface and the point at which crack 
growth becomes unstable. Based on an examination of all 
available information, the most promising methods for 
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attacking the cracking/fracture problem have been selected 
and are listed as follows: 

a. Crack initiation and propagation: 

i) For predicting crack initiation, the 
maximum tangential stress criterion 
will be used. 

ii) For crack propagation and for pre- 
cracked problems, the minimum strain 
energy (S) criterion will be em- 
ployed 11 ,2] *. 

b. The use of the minimum strain energy criterion 
will facilitate the process of predicting the 
direction and the rate of crack propagation^ 

c. Stress intensity factors will be evaluated 


13 ] 


using the displacement method, through the use 

r 4 1 

of crack tip singular elements with a layer 
of transition elements to improve the re- 
sults ^ . 

d. At bimaterial interfaces, the problem will be 
solved partially by introducing material dis- 
continuity effects in evaluating the stress 
intensity factors determined from the numer- 
ical results on the crack flanks. 

e. To determine crack growth stability, two 
criteria are proposed and listed as follows: 


* 

Numbers in brackets m this section refer to the references 
which are listed on page 122, 
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i) Compare the size of the crack with 
some typical (global) dimensions of 
the problem. 


ii) Compare the crack growth increment 


(given for a stable crack as 


Aa (n + 1) 
Aa 


n 

where Aa n is the crack growth incre- 
ment in step n) with a stability 
constant. 


A summary of possible crack problem algorithms and the 
reasons for selecting those listed above is given in [6] . 


4 . Plans for Upcoming Period 
During the next reporting period we will concentrate 
our efforts on implementing the various computational mod- 
ules of the crack problem. One such module, not yet solved, 
is the remeshing of the finite elements as cracks move 
through a region. This aspect will be considered. We ex- 
pect the automatic mesh generation capability, discussed in 
Section IV-D of this report, to play a central role in the 
solution to this problem. 
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INTRODUCTION 


The redesign of f light -critical structural components 
currently in service as metal parts on existing aircraft was 
included as part of RPI ' s composite materials and structures 
program to assure relevance and realism to the project's di- 
rection and to contribute to the nation's pool of innovative 
structural design concepts . Two such components were se- 
lected in cooperation with major commercial airframe manu- 
facturers, Boeing Commercial Airplane Co. and Lockheed Cali- 
fornia Company, and with advice from NASA. The first, the 
727 airplane's elevator actuator attachment rib, has com- 
pleted the cycle planned at the project's inception. That 
is, several alternate redesigns were attempted; one was se- 
lected, analyzed and finalized; and two specimens were fab- 
ricated and successfully tested to ultimate load. The sec- 
ond, a redesign of the wing-mounted engine drag strut for 
the L-1011 airplane, has been carried to the point where (a) 
the usefulness of certain basic research performed earlier 
at RPI under the NASA/AFOSR program was made very clear and 
was called on in the design process, and (b) several areas 
of research in which results are needed for such a redesign 
- but are not now available - have been identified. 

As a matter of fact, a number of basic research projects 
motivated by the early work on the 727 elevator rib, have 
been reported in these project reports for some time. Some 
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of these have been completed, others are reported here and 
still others were planned based on challenges encountered 
in the 727 elevator problem. Both these research projects 
and those arising out of the L-1011 engine drag strut rede- 
sign were selected for study because they were seen as hav- 
ing the potential to contribute significantly to classes of 
generic composite aerospace structures. It was decided in 
November 1981 to suspend .any further conceptual design re- 
search on actual aircraft components, in favor of pursuing 
basic studies of such generic structural elements. Progress 
on these projects during the current program period is re- 
ported in the sections that follow. 
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IV- A MECHANICAL JOINTS IN COMPOSITES 

The 727 elevator actuator attachment rib redesign in 
composites at RPI , reported in earlier progress reports, led 
to separate research projects. One deals with heavily load- 
ed holes in composites such as are encountered at the actu- 
ator attachment; the second deals with lightly loaded mechan- 
ical joints, such as occur, for example, in a row of small 
bolts or rivets, typified by the more or less continuous 
joint between rib flange and the elevator skin. Both of 
these problems have required the development of new analysis 
and/or test methods. 

1 

More recently, the importance of optimizing the strength 
of compact lugs has become clear. All three of these projects 
are described in the paragraphs to follow. 

1. Analysis of Heavily Loaded Mechanical Joints 

Senior Investigator: C. Muser 

a . Introduction 

Discontinuxties such as mechanical joints and openings 
are always going to be found in structures independent of the 
materials used. Although fiber reinforced materials have the 
potential to yield efficient designs for such discontinuities, 
these materials are used in only a very few cases for heavily 
loaded mechanical joints. It is believed that this is due to 
a lack of knowledge of how to use compos xtes efficiently for 
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such discontinuities and, therefore, a detailed investigation 
is justified. 

b. Status 

It had been decided to focus on cylindrically orthotro- 
pic materials rather than on materials orthotropic with res- 
pect to cartesian coordinates. A closed form solution for 
the stress distribution in a circular plate with a circular 
hole at its center was found by Dr. N. J. Hoff for the case 
when its material is cylindrically orthotropic . The geome- 
try of that plate can be seen in Figure IV-C-l-a of the pre- 
ill * 

vious report , Later, this analysis was extended to treat 
a similar plate whose elastic properties, however, were a 
function of the distance from the hole. The constitutive 
equations used were: 

E a. = S i j (1 + C<i)( j > rh > 0 i 

Y re - S 66 U + Q 66 rh)T r9 

where i, j = r,0 and () indicates no sum on indices. 

The constants St. and Q. . are a function of the dimen- 

1 J -*-J 

sion of the plate and the parameter h, which describes the 

rate of change of material properties along the radius r. 

The compliances, S. ., can be prescribed at both boundaries 

J 

independent of the choice of the parameter h, so long as the 
convergence criteria are observed. 

_ - . 

Numbers m brackets in this section refer to the references 
which are listed on page 133. 
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It was shown that cylindrical orthotropy can reduce the 
stress concentrations along the edge of the hole substantial- 
ly, especially when the stiffness of the fibers is much higher 
than that of the matrix. The same behavior can be observed 
for the case with varying material properties. In the latter 
case, however, it is possible to link the properties yielding 
low stress concentrations with those normally encountered in 
designing a plate as a structural element. 

c. Progress During Report Period 

Substantial effort was devoted to reporting the results 
obtained, as described in the previous section. One paper 
has been accepted for publication and a second is going to 
be presented at an international composites conference (see 
Section e. Current Publications or Presentations, etc., 
page 133) . 

Progress was also made in applying the same analytical 
methods to the problem of a circular plate pulled at one end 
in one direction, with that load reacted by a bolt inserted 
in the centrally located hole (see Figure IV-A-l-b of the 

rn 

previous report ) . The material properties of that plate 
were allowed to vary along the radius . 

As mentioned in the previous report, this work does not 
attempt to find the exact boundary conditions at the edge of 
the hole. The goal is rather to develop an algorithm which 
is capable of handling, to a satisfactory degree of approxi- 
mation, any conditions at both boundaries of the plate. 
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Therefore, the assumptions commonly used, i.e., there is no 
friction between the bolt and the plate and there is no gap 
between the bolt and the plate, were considered to apply. 
The boundary conditions were expanded in Fourier series. 
This implies that there are n sets of boundary conditions, 
if for every stress acting on the boundaries, n terms are 
considered. For every such set, a solution can be found. 
The total of all n solutions found in such a manner yield 
the solution for the bolt-loaded plate. 

A typical form of the stress function, , used is: 


— oo 4 n oo in + hin 

$ = E E A E K r P cos n0 

— n=l p=l P m=o P 


The coefficients K 1 * 1 can be determined with a recursion formu- 

P 

la which emerges from the process of satisfying the compati- 
bility equations. (There is one compatibility equation for 
every n, where n is the index of terms considered in the ex- 
pansion of the boundary conditions.) The constants A n can 

P 

be determined with the aid of the n— th set of boundary con- 
ditions. The roots, t are those obtained from the case of 

Jr 

homogeneous materials, for which a closed form solution had 


been found, as mentioned earlier, 
a function of n. 


These roots, t , are also 
P 


A computer program, which is in the process of being 
tested, will be used to calculate the stresses m a circular 
bolt-loaded plate whose elastic properties vary radially. 
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d. Plans for Upcoming Period 
The investigation of stress distributions in circular, 
bolt-loaded plates, the materials of which are cylindrically 
orthotropic and vary radially, is planned for completion in 
the next reporting period. 


e. References 


1. 41st Semi-Annual Report, Composite Structural Materials, 
AFOSR and , NASA Grant No. NGL 33-018-003, Rensselaer 
Polytechnic Institute, Troy, NY, December 1981, 117 and 
120 . 


f. Current Publications or Presentations by 
C. Muser on this Subject 


"Stress Concentration Factors for Cylindrically Orthotropic 
Plates, with N. J, Hoff. 

To be published in the Journal of Composite Mate- 
rials, May 1982. 

"Stress Concentrations in Cylindrically Orthotropic Plates 
with Radial Variation of the Compliances", with N. J. Hoff. 

To be presented at the Fourth International Confer- 
ence on Composite Materials (IV-ICCM) , Tokyo, Japan, 
October 25-28, 1982. 
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2. Pin-Loaded Holes in Uniform Composite Plates 

Senior Investigator: R. G. Loewy 

a. Introduction 

Stress analysis is a key step in the design of an op- 
timum mechanical joint in composites as it is in the design 
of any other structural member. To study this problem, an 
idealized model of a single fastener joint has been postu- 
lated (see Figure IV-A-2-1) . Note the existence of pin-hole 
clearance, S. This gives rise to uncertain boundary condi- 
tions , particularly as to the amount of the contact along 
the pin hole boundary under load, as explained in previous 
reports. These unprescribed boundary conditions preclude 
application of ordinary stress solution methods to this prob- 
lem. 


1 Preceding page blank 


b. Status 

The work reported in this section is that of graduate 
student Wonsub Kim. 

A linear constraint equation has been introduced in the 
finite element method (FEM) analysis of this problem to des- 
cribe the motion of pin-hole contact points as the hole 
boundary deflects under load. These constraint equations 
were constructed on the basis that any boundary point within 
the region of pin-hole contact must remain on the surface of 
the pin (assumed rigid and frictionless) for a given applied 
load level. Thus, the only remaining unknown regarding 
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boundary conditions is the extent of the pin-hole contact. 
Next, a "mixed method" was introduced, consisting of the 
finite element analysis (FEM) as an analytic part and photo- 
elastic experiments as an empirical counterpart. In the 
mixed method, the unknown extent of the contact region 0 

c 

is assumed and FEM analyses iterated until a value of 0 c 
can be deduced which produces a distribution of maximum 
shear strain which matches that produced by the photoelastic 
experiments. This contact angle is assumed to be close to 
the actual one, and the numerical method (FEM) is then 
available for all aspects of a solution. 

Two advantages of the mixed method over more conven- 
tional methods are evident. First, it is a self -checking 
method, since experimentally determined results are reflected 
in the analysis and some confidence can be gained as to the 
mathematic modeling and the theory used in the analysis. 
Second, extrapolated analyses can be conducted at failure 
load levels, and reliable failure prediction studies are, 
therefore, possible. 

In the previous reporting periods, analyses were per- 
formed at several load levels for the so-called B-type spe- 
cimen ([±45°/0°] ) to produce the curve relating applied 

ns 

force versus assumed contact angle, which is shown again in 
Figure IV-A-2-2, for convenience. Also, results were obtained 
in photoelastic experiments for laminate specimens type B, 
[±45°/0° 3 0 and type F, [±30°/0°] . 

lib IlO 
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c. Progress During Report Period 
Analyses for the types B and F laminate specimens were 
completed at a number of additional load levels in this re- 
porting period. In addition, two new analysis approaches to 
the pin-hole contact problem were devised, each of which has 
its own distinct advantages. 

An attempt was also made to predict failure loads by 
introducing two distinct existing failure criteria, using 

the FEM stress computation results. These criteria are 0. 

Ql* 

Hoffman's failure criterion and the so-called "maximum 

[21 

stress criterion" . Correlation attempts showed that both 
criteria predicted correct failure modes but were too con- 
servative in predicting joint strength. 

In the following paragraphs, the two new analysis meth- 
ods will be explained, the results of the analyses summarized 
and failure prediction attempts described. 

1) Interactive Method and Energy Criteria 

Figures IV-A-2-3a, b and c show deformed shapes of the 
fastener hole as dotted curved lines for cases in which the 
assumed contact angles are correct, too small and too big, 
respectively. In the same figures, the shape of the pin as 
recognized by the FEM through the linear constraint equations 
(which is the only way the FEM can recognize the existence 


Numbers m brackets in this section refer to the references 
which are listed on page 171. 




J 


1 1 J t J i i l l [ l l l 1 1J.J-4J.W 


Applied Load 


Applied Load 


Applied Load 


Figure IV-A-2-3. Deformation of Hole Due to Assumed Contact Angle (0 ) 
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of the pin) are also depicted. Comparing Figures IV-A-2*-3a 
and -3b, one can see that the Figure 3b case lost some sup- 
port (between points P e , assumed end of contact, and P, the 
actual end of contact) nominally provided by the pin, which 
is fixed in space. In fact, there is an inadmissable inter- 
ference between the pin and the deformed hole edge configur- 
ation. If the pin were not shown, however, there is still 
an identifiable kink in the deformed boundary near point P^. 
From the lost support there results a bigger displacement, 
compared to Figure IV-A-2-3a, at the end of the plate where 
the external tensile load is applied, although this doesn't 
show in Figures IV-A-2-3. Thus, more work is done by the 
applied force, and it can be expected, in a conservative sys- 
tem, that this will result in more strain energy stored in 
the plate. 

On the other hand, by comparing Figures IV-A-2-3a and 
-3c, it can be seen that the effect of assuming to be 

pa C’* t 

bigger than the actual contact angle, 0^ , some hole bound- 

ary points located between points P and P g , are forced un- 
realistically to the surface of the pin. This produces the 
inward hump on the deformed boundary shown in Figure IV-A- 
2-3c. If one imagines that forces on the hole boundary be- 
tween P and P g have done work in forcing these points to the 
pin contour, then, in this case too, the total external work 
and the strain energy stored in the plate would be bigger 
than that in Figure IV-A-2-3 a. 
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This leads to the conclusions that: 

« * 

x) a smooth hole edge deformation shape is pos- 

9.C t 

sible only when 0 =0 and 

c c 

ii) the strain energy stored in the body and the 

potential of (in linear elasticity, two times 

the work done by) applied forces both have 

"t 

a relative minimum when 0 c = © c 


if 0 


c 



Either of these two facts can be used as a criterion to 

determine the correct value for 0 in clearance-fit joint 

c 

analysxs. One such example may be seen in Figure IV-A-2-4. 
This shows the deformed shape of a FEM mesh plotted over the 
original, undeformed mesh for cases in which the assumed 
0 c ’s are equal to (a) 50°, (b) 65° and (c) 70° when a load 
of 2,000 pounds is uniformly applied to a B-type specimen. 
The original mesh is shown as solid lines and the deformed 
mesh as dashed lines. The original, undeformed hole boun- 
daries are shown as dashed lines and the deformed hole boun- 
daries by dot-dashed lines. In these plots, displacements 
are exaggerated by factors of 50, 100 and 100 times the ac- 
tual computed amounts for figure parts (a) , (b) and (c) , 
respectively. As indicated by arrows in the figures, the 
deformed shape of the hole edge boundary of case (a) has a 
kink (0 < 0 ) and case (c) has a kxnk (0 >0 j while 

O C CO 

case (b) shows none. This indicates that the actual contact 
angle is close to 65° xn this example. Indeed, this is con- 
firmed in the applied load versus assumed contact angle 




Figure IV-A-2-4 . Deformed Shape of the Hole Boundary (dot-dash lines) as a Function of Assumed 

Half-Contact Angle (0 ). (B-type [±45°/0°j laminate; applied load 2,000 lbs.) 
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curve of Figure IV-A-2-2, where 0’ is found to be about 63° 
for an applied load level of 2,000 pounds. 

The second approach to a purely analytical determina- 
tion of 0 c is shown in Figure IV-A-2-5, which presents the 
potential of the applied force as a function of assumed con- 
tact angle. In this figure, potential variation in terms of 
assumed contact angle can be seen for four cases of load and 
lamination combinations. The existence of a relative mini- 
mum of potential variation is clearly indicated. Moreover, 
the minima for the cases B4000 and B6000 in Figure IV-A-2-5 
fall at values of 0 c almost identical to those which Figure 
IV-A-2-2 gives for the corresponding loads. 

The above discussion leads to the conclusion that both 

the smoothness of the hole edge deformation and the minimum 

property of potential can be used to determine 0 . 

c 

2) Some Results of Analysis 

Both the semi-empirical and the interactive methods 
described in the previous sections were applied to a single- 
fastener joint model (Figure IV-A-2-1) for two different 
types of graphite-epoxy laminate. These laminates were 
[±30°/0°] or F-type (24 layers) and [±45°/0°] or B-type 
(23 layers) . The elastic properties of the material are 
shown in Table IV-A-2-1, and the strength data is given in 
Table IV-A-2-2. 

In the finite element idealizations, the layers of F- 
type laminates were rearranged as follows: [ (+30° ) _/ (0°) _/ 

O O 

(-30°) g 3, whereas, the layers of B-type laminate were 




Figure IV-A-2-5. Energy Potentials of Applied Loads as a Funtion of Assumed Half-Contact Angle O c ) 



E: Energy potential of applied loads undergoing displacements 

^ ^ 4 

B and F indicate the type of lamination (B: [±45°/0°] , F: 

, ns 

numbers indicate the levels of applied load. 


(lbs ./in. ) 

[±30°/0°] ). 

ns 


The associated 





146 


' table IV-A-2-1 

PROPERTIES OF GRAPHITE/EPOXY TAPE 
(Fiberite Hy-E 1048A1E) 

* 

Properties Symbols Values 

Modulus in Fiber Direction E T 21.306 

* Jj 

(10 6 psi) 

Modulus Transverse Direction E^ 1.4 

(10 6 psi) 


Shear Modulus 
(10 6 psi) 


LT 


0.93 


Poisson’s Ratio 


V 


LT 


0.32 


* 

For a fiber volume fraction of 60%. 


TABLE IV-A-2-2 

STRENGTH PROPERTIES OF GRAPHITE/EPOXY COMPOSITES 


Property 


Explanation 


Value 

(ksi) 


X 


X’ 


y 


y* 


s 


Tensile strength 

in fiber direction 

Compression strength 
in fiber direction 

Tensile strength 

transverse direction 

Compression strength 
transverse to fiber 

Shear strength 
(interlaminar) 


200.0 

>1 4 * 

207.0 

8.0 

18.0 

A* 

19.0 


Hy-E 1048 AiE prepreg tapes (60% fiber volume frac- 
tion) , Fiberite Co. 

These values were suggested by Prof. Diefendorf of the 
Dept, of Materials Engineering, RPI, based on his ex- 
periences with similar fiber-resin systems. The rest 
of the properties were given by Dr. J. Allen of Fiberite 
Co., the manufacturer of the materials, through telephone 
communication with Wpnsub Kim on January 12, 1982. 
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represented as [ (+45° ) _/ ( 0° ) _/ ( -45° ) ] . This was possible 

without loss of accuracy, since laminate plate theory for 
in-plane tension problems recognizes only average values 
through the thickness. 

The isochromatic fringes of Figures IV-A-2-6a, b and c 
showed maximum shear strains which were well-matched with 
the distribution of the corresponding contours generated by 
the FEM analysis (see Figures IV-A-2-7a, b and c, respective- 
ly) . The isochromatic fringe of Figure IV-A-2-6c, for the 
B-type laminate at 6,000 pounds applied load, was obtained 
immediately before failure occurred. Isochromatic pictures 
of the F-type laminate at failure load, 5,760 pounds, however, 
which would have corresponded to Figure IV-A-2-7d, could not 
be taken, since failure took place suddenly (the usual case) 
and before a fringe picture could be taken. 

It can be seen in Figure IV-A-2-6b that the symmetry of 
the fringe pattern of the F-type specimen was disturbed, 
due to small twisting moments which came about because of 
the rotation of the spindle of the Instron testing machine. 
These rotations were not controllable. As a result, ring- 
type isochromatic fringes appeared just above the hole in 
the figure, in which the pin pushes upward. Accordingly, 
the analysis for the F-type laminate was done purely by the 
interactive method. Analyses for the B-type laminate, how- 
ever, were done using both the semi-empirical and interactive 
approaches. Results from both methods for this type of 
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Figure IV-A-2-6. Isochromatic Fringes: (a) B4000, (b) F4000 and (c) B6000 

(one fringe value Y max = 0.0018917 in. /in.) 
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Figure IV-A-2-7. Maximum Shear Strain Distribution Contours from FEM Analysis (a) B4000 & (b) F4000 

(one fringe value y = 0.0018917 m./in.) 
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laminate were almost identical. The results from the inter- 
active method were not included here to avoid duplication. 

It should be noted that the semi-empirical method does 
not require an iterative approach to the failure load condi- 
tion, as regards pin contact angle 0 c , when the analysis has 
been performed at several other load levels, as explained in 
the previous report and indicated in Figure IV-A-2-2. Calcu- 
lations for the failure load level with several different 
assumed contact angles were performed, however, to check the 
smoothness of the deformed contours of the hole edge and the 
level variation of potential with assumed contact angle, a 
necessary procedure for the interactive and minimum potential 
methods, respectively. 

The nomenclature used in Figures IV-A-2-8a through -8d is 
explained in Table IV-A-2-3. These figures show the radial 
and circumferential stresses in each layer of the laminate 
as well as the average values across the thickness of the 
plate along the edge of the hole. Comparing Figures IV-A-2-8a 
and -8b, one observes clearly the nonlinear nature of the 
problem in its various aspects, even though pin-hole clear- 
ance is as small as 0.0015 inch. For example, when the ap- 
plied load level on the B-type laminate increased from 4,000 
pounds (see Figure IV-A-2-8a) to 6,000 pounds (see Figure IV-A- 
2-8b) , the maximum circumferential stresses more than doubled, 
while the peak values of radial stress increased by only rel- 
atively small amounts. Furthermore, the rate at which stress 
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TABLE IV— A— 2— 3 

NOMENCLATURE EXPLANATION FOR FIGURES 8 THROUGH 10 


Symbol 

a . (ksx) : 
net 

CT, (ksx) : 
b 


a 

(ksx) : 

rr 


a ee 

(ksx) : 

a (psi) : 

max e 

a 

(psx) : 

mm 

T 

(psx) : 

max 

°L 

(psx) : 

°T 

(psx) : 

t lt 

(psx) : 


T max (in - /ln/): 
B, P: 

r (xnch) : 
t (xnch) : 


Expl anat ion 

Average net tensxon stress = Applxed load (P) /net cross- 
sectxonal area 

Average pin-bearxng stress = Applied load (P) /pxn-bear- 
ing area (2rt) 

Radxal stress on the hole-edge boundary 

Cxrcumferentxal stress on the hole-edge boundary 

Maxxmum prxncxpal stress 

Minxmum prxncxpal stress 

Maxxmum shear stress 

Stress in materxal (fxber) directxon 

Stress xn transverse fxber dxrectxon 

Shear stress in material coordinates 

Maxxmum shear straxn 

Lamxnatxon type: 

B: [±45°/0°] , F: [±30°/0°] 

ns ns 

Radxus of pxn 
Thxckness of the laminate 


Pxcture Identif icatxon Txtle 


Example: 


B4000 - (-45°) - a 


max 


Plotted quantity 

Fxber orientatxon of the layer 
(ALL) represents whole lamxnate 
or all layers 


Applxed load xn pounds 
Type of lamxnate 
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increases with applied load differs from layer to layer. 
Similar trends occur with F-type laminates (see Figures IV-A- 
2-8c and -8d) . 

Comparing Figures IV-A-2— 8a and -8c shows that the peak 
value of circumferential stresses in each layer, as well as 
the averaged value through the thickness of the whole plate, 
is higher for the F-type laminate than those for the B-type 
laminate at an applied load of 4,000 pounds. This trend, 
however, is reversed for radial stresses. Both the circum- 
ferential and radial stresses occurring in the F-type lami- 
nate are spread more evenly among its layers than those of 
the B-type, except for radial stresses in angle-ply layers. 

On the other hand, comparisons between Figure IV-A-2-8b for 
the B-type laminate under an applied load of 6,000 pounds and 
Figure IV-A-2-8d for the F-type with 5,760 pounds load ap- 
plied, show that the peak circumferential stress in the B- 
type’s zero degree layers is higher than that of the F-type 
laminate; this despite the fact that the peak of the average 
circumferential stress in the F-type remains higher than that 
of the B-type. Radial stresses display basically the same 
trends for higher load levels as at lower level loads com- 
paring the two laminates. The average stresses, however, 
are again exceptions. The maximum value of the average ra- 
dial stress in the F-type plate is about the same as that of 
the B-type laminate; namely, about 103,000 psi in compres- 
sion, while those at lower load levels differ between the two 
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types of laminates. This is particularly interesting, in 
view of the fact that the maximum average radial stresses 
at failure were about the same for both types of laminates, 
and both specimens failed in the same failure mode, i.e., 
bearing. 

Figures IV-A-2-9a through -9d are plots of stresses in 
the direction of load application, i.e., normal stresses on 
the net tension cross-section. Data for both types of lami- 
nate are shown at two levels of applied load. In these fig- 
ures, average tensile stress across the section (calculated 
by dividing the applied load by the net cross-section area) 
is indicated for each case, to allow rapid estimation of 
stress concentration factors, rather than show separate plots 
of stress concentration factor. Surprisingly, tension stress 
concentration factors in both laminates increased with load 
increases, in spite of increasing contact angles. Further- 
more, the growth in stress concentration for B-type laminates 
is bigger than for the F-type. 

Similar plots were made for compression stresses in the 
direction of applied load along the line of symmetry, as shown 
in Figures IV-A-2-10a through -lOd. In these figures, the 
average pin-bearing stress, which was obtained by dividing 
the applied load by the pin-bearing area (diameter of the 
pin multiplied by the plate thickness) , is included for each 
case. It is noteworthy, as may be seen in Figures IV-A-2-10c 
and -lOd, that the maximum compression stress value in the 
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loading direction along the line of symmetry occurred a small 
distance away from the edge of the hole for the F-type spe- 
cimen. This phenomena had also been observed in the photo- ' 
elastic experiments with several types of laminates. 

According to Figure IV-A-2-10a, the compression stress 
for the B-type laminate in zero degree layers at the edge of 
the hole (204,000 psi) has almost reached the compression 
strength limit of the material as given in Table IV-A-2-2 
(207,000 psi), at a load level of only 4,000 pounds. At 
6,000 pounds applied load (see Figure IV-A-2-10b) , zero- 
degree layer compression stress exceeds by far the stress 
allowable of this component. The actual failure, however, 
occurred at 6,000 pounds. Figures IV-A-2-10c and -lOd for 
the F-type laminate, on the other hand, show compression 
stresses at the failure load which check reasonably with the 
material allowable data given in Table IV-A-2-2. 

Similar disagreement between calculated tensile stresses 
at failure load and the material tensile strength data is 
even more severe, as seen in Figures IV-A-2-8 and -9. 

Examination of Figure IV-A-2-8 reveals that the absolute 
maxima of circumferential stress along the edge of the hole 
take place roughly at the points where fibers are tangent to 
the hole boundary. In contrast, the absolute maxima of ra- 
dial stress occur close to the points where fibers are per- 
pendicular to this boundary. 
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The small non-zero value of average radial stress seen 
in Figures IV-A-2-9 on the free edge of the hole, at points 
which are beyond the contact region defined by 0 , requires 
some discussion. Radial stresses should be zero along the 
traction free portion of the hole edge. Therefore, these 
non-zero radial stresses must be regarded as computational 
inaccuracies. The reason for their existence is as follows: 
These stress plots were based on the nodal values along the 
hole boundary. The nodal stresses on this boundary were cal- 
culated from the element centroidal values of each stress 
component in material coordinates by the smoothing technique 
explained in previous reports. The nodal values were then 
transformed into the polar coordinate system which would pro- 
vide radial and circumf erential components of stress for each 
layer. This interpolation (actually closer to extrapolation 
for the hole boundary) causes small errors for each component 
of stress in each layer. These errors accumulate through co- 
ordinate transformations and averaging processes. Thus, the 
non-zero values of the stress resultants at these hole-edge 
nodes, in deviating from traction free boundary conditions, 
may provide an indication of the upper bound for the kind of 
numerical errors likely to be found in the FEM analysis re- 


sults . 
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3 ) Fracture Prediction 

Attempts, which must be considered preliminary, were 
made to examine the ability of two distinct criteria, i.e., 
the 0. Hoff man and the maximum stress criteria, to predict 
ultimate joint strength based on the computed stresses. 

Stresses at each node and at each element centroid were 
transformed into material coordinate systems for each layer, 
and stress components, a , a and x , were put into the 

Jj T L I 

selected failure criteria. For Hoffman's criterion, the 
value of n, as given by the following equation, was calculated 
at each node and element centroid. 
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where, material strength properties are, as also shown in Table 
IV-A-2-2 , 

X = tensile strength in the fiber direction 
X' = compression strength in the fiber direction 
Y = tensile strength transverse to the fibers 
Y' = compression strength transverse to the fibers 
S = inplane shear strength. 

After this, every node and centroid of every element where 

n > i 

was, marked on the FEM mesh model, which was displayed on the 
screen of the graphics computer terminal. These points indicate 
where failures are predicted to occur according to the Hoffman 


criterion. 
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A similar process was followed using the maximum stress 
criterion. In this case, the failure marks were put on the 
points where one of the following conditions were met: 


i) 

CS r 

> 

X 


L 




ii) 


< 

o 

and 1 cr T 1 > X ' 

Ii 




iii) 


> 
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T 




iv) 
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o 

and | cj t |> Y* 

v) 

1 X LT 

l> 
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Some of the results of such failure prediction routines 
are shown in Figures IV-A-2-lla through -llg, in which the 
laminate type, applied load level, and the failure criterion 
used are indicated. Failure plots for the B-type specimen 
at 4,000 pounds applied load, based on the maximum stress 
criterion were not included, because the prediction indicated 
no failure in this case. 

Judging from these plots, both failure criteria were too 
conservative in predicting ultimate strength allowables for 
the joints, but they were correct in predicting the modes of 
failure . 


d. Plans for Upcoming Period 
In the next reporting period, we will concentrate our 
efforts in extending the methods reported here to the anal- 
ysis of multi-hole fastener joints. Consideration will also 
be given to similar analysis of pin-hole joints with inter- 
ference fits. 
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(a) B4000 - (All) 


Figure IV-A-2-11. Failure Prediction: (a) B4000 (All) 

(Small circular marks on nodes and 
elements indicate failure occurred on 
those points.) 
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Figure IV-A-2-11. (Continued) 
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Figure IV-A-2-11. (Continued) 
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3, Compact. Lug Design 

Senior Investigator: D. B. Goetschel 

■» 

a . Introduction 

A critical aspect of many aerospace structural elements 
xs the load transfer that takes place between the connecting 
lugs at the ends of the structure and the portions of the 
structure wherexn the loads are well distributed if not uni- 
form. An example is provxded by the engine drag strut of 
the L-1011 aircraft (see Figure IV-A-3-1) . Such lugs are 
hxghly loaded and have a very complex stress state (see Fig- 
ure IV-A-3-2) . Further, since they must mate with connect- 
ing parts, these lugs are usually designed within rather 
stringent dimensional envelope constraints. As compared to 
lugs which, as in the case of the L-1011 strut, are made 
from 200 ksi steel, meeting the geometric constraints even 
with a composite structural design making maximum use of uni- 
directional graphite-epoxy could prove to be a difficult 
task. This research is xntend.ed to solve the load/volume 
problem in favor of composites rather than either reverting 
to designs using 200 ksi steel lugs attached to a composite 
strut or forcing the dimensional constraints to be relaxed 
in a redesign of the mating parts. 

b. Status 

The design of lugs which are capable of L-1011 drag strut 
loads has been undertaken. Because a flat lay-up appears to 







MAXIMUM TENSILE LOAD - 207.7 KIPS 
MAXIMUM COMPRESSIVE LOAD - 236.8 KIPS 
MAXIMUM LATERAL LOAD - 12.0 KIPS 

DIMENSIONS IN INCHES 

Figure IV-A-3-2. Lug Experiment: Full-Scale L-1011 Drag Strut 
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be the easiest and cheapest to manufacture when compared to 
other proposed designs, we decided, first, to try to use a 
flat, symmetric zero-plus-mmus-angle-ply laminate whose 
lamina each consist of unidirectional graphite epoxy. Uni- 
directional graphite epoxy was selected because of its ready 
availability, reasonable cost and high strength in both ten- 
sion and compression. 

A zero-plus-mmus-angle-ply lay-up has the potential of 
being optimized to resist the potential failure modes of 
bearing shear, compression or tension. The second decision 
was to concentrate initially on the tension design load, 
since it produces the most complex stress state with large 
compressive, shear and tensile stresses. If this design 
loading can be handled, then it should be relatively easy to 
handle the compression design load. 

c. Progress During Report Period 

An initial exploratory testing program was conducted by 
senior undergraduate student Anthony Staniorski to examine 
the behavior of the proposed lug design under tensile loads. 
To accomodate testing machine capacities and conserve on 
fabrication costs, quarter-scale model specimens, as com- 
pared to the L— 1011 strut, were tested. Simple scaling of 
the ultimate design load requires the quarter-scale models to 
withstand 13 kips in static tension. This is well within 
the testing machine capacities. 
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Three specimens were made, each with a different lay-up. 
The material properties of the graphite-epoxy unidirectional ' 
tape used are listed in Table IV-A-3-1. The lay-ups were 
chosen somewhat arbitrarily but with the intent of generat- 
ing various types of failure. The lay-ups and failure loads 
are listed in Table IV-A-3-2. In all cases, the failure 
mechanism was tensile failure in the net section through the 
pin hole. As can be seen, in comparison with the required 
strength of 13 kips, large improvements were needed. 

In parallel with these tests, computer programs were 
written by senior undergraduate students William Dhimitri 
and Paul Adel which would facilitate the design of an opti- 
mum lay-up. As a first approximation, very simple failure 
criteria were used. It was assumed the lugs would fail m 
either pin bearing, pin shear -out or net tensile failure. 

For each condition an average stress state was calculated 
using simple free body diagrams. These average stresses were 
then multiplied by stress concentration factors chosen intui- 
tively. The quadratic interaction criterion was applied to 
these stresses as a final step. These programs predicted 
that all the specimens would fail m either bearing or shear - 
out and that the specimens should have wide margins from 
tensile failure. Since the actual failure mode was, in fact, 
in tension, it was obvious that improved analysis and design 
modifications were called for. 
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TABLE IV-A-3-1 

MATERIAL PROPERTIES OF GRAPHITE/EPOXY USED IN TESTING 


Longitudinal Modulus = E 


Transverse Modulus = E 

Shear Modulus - E 

Poissons Ratio = v 

Ply Thickness = h 


Longitudinal Tensile 
Longitudinal Compressive 
Transverse Tensile 
Transverse Compressive 
Shear 


= 21,300 ksi 
= 1,400 ksi 
= 930 ksi 

= .32 

= .005 inches 

Strength = X = 200 ksi 
Strength = X' = 170 ksi 
Strength ~ Y = 6 ksi 
Strength = Y ' = 36 ksi 
Strength = S = 10 ksi 


TABLE IV-A-3-2 

FAILURE LOADS OF LUG TEST SPECIMENS 


Sample 

Number 

Failure 

Load 

(kips) 

Percent 
of Design 

1 

7.3 

56 

2 

6.1 

47 

3 

9.6 

74 

Design 

Ultimate Load 

= 13 kips 
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Two explanations and corrective actions were postulated 
for the unexpected weaknesses in tension. First, it was ob- 
served that at the sides of the hole, where high normal 
stresses are expected, there is a free edge with associated 
interlaminar stresses. To alleviate these interlaminar 
stresses, it was proposed to bond a thin bushing into the 
pin hole. This should provide some stress continuity from 
lamina to lamina at the otherwise free edge. The bushing 
was made of brass with a wall thickness of .015 inches. Sec- 
ond, it was suspected that the normal stress concentration 
factors at the sides of the hole were much greater than ex- 
pected. Slots were suggested as shown in Figure IV-A-3-3 as 
one means to reduce this stress concentration. By channel- 
ing the load away from the hole, instead of requiring it to 
"flow" immediately behind the hole, the tendency for stresses 
to concentrate next to the hole should be reduced. 

A second round of exploratory tests was performed to 
evaluate these design modifications. Some improvements were 
noticed, but the detailed data will not be reported here 
since variations of results between tests greater than the 
differences attributable to design changes appeared to be 
due to random causes. Multiple samples of each design and 
better consistency in specimen fabrication will be required 
to get more meaningful data. 

The failure mode of the modified samples did, in gen- 
eral, tend to be pin shear -out, however, rather than net 
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tension. Perhaps the most encouraging result was that one 
slotted specimen reached 90% of the ultimate design load, and 
there is no reason to believe that design is very close to 
optimum. 

Finite element analysis was used to examine the ten- 
sile stress at the sides of the pin hole for both slotted 
and unslotted configurations. It was found that slots should 
reduce the tensile stresses next to the hole by thirty to 
forty percent. While this is a dramatic reduction, one would 
not expect to see a corresponding increase in lug strength, 
since other failure modes are likely to govern before the 
full improvement is realized. Other failure modes were not 
examined in the FEM analysis because an unreasonable amount 
of postprocessing of the results had to be done by hand. 

FEM analyses were also used to examine the variations 
that changes in ply lay-up had on net tensile failure. Both 
the percentage and angle of the angle plies were varied. It 
was found that ply lay-up changes did dramatically alter the 
stress concentration next to the pin hole, but these alter- 
ations were roughly proportional to the variations in the 
lay-up 1 ' s tensile strength. The ultimate result was that net 
tensile strength for the lug appears to be only a weak func- 
tion of ply lay-up changes. Here again, other failure modes 
were not examined because of the large amount of hand pro- 
cessing of the data that was required. 
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d. Plans for Upcoming Period 
The optimums that were found in the FEM analyses to 
date will be used as a guide in designing the next sequence 
of tests. The finite element capabilities, however, must be 
improved. Work will begin on automating the postprocessing 
needed to calculate composite material failure criteria so 
that all potential failure modes can be examined. Once the 
composite material failure criterion calculation is auto- 
mated, various criteria will be tested for accuracy. It is 
important to have much more complete results of stress anal- 
ysis to use as a guide in designing new tests. 

All tests are being planned using larger numbers of 
samples to handle the statistical variations in strength 
that occur within a given configuration. In addition, more 
consistent fabrication methods will be developed. 
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IV-B SCALING EFFECTS IN TESTING COMPOSITE STRUCTURES 

Senior Investigator: D. B. Goetschel 

X . Introduction 

It can be extremely valuable in terms of both time and 
money to be able to test small scale models when designing 
large structures. This is especially true for structures 
being designed in composite materials such as graphite epoxy, 
since the material itself is relatively expensive. As a 
first step in investigating scaling effects in design develop- 
ment testing of structures using composite materials, a study 
of scaling of simple tensile specimens was conducted. This 
study was primarily conducted by senior undergraduate stu- 
dent Scott Baxter . 


2 . Status 

Two main effects v/ere identified and hypothesized as 
being the mechanisms resulting in failure stress variations 
due to scaling. First, Weibull statistics are expected to 
be a major factor. These probability distributions are used 
in an attempt to estimate the effect that flaw distribution 
in a material has on the failure load of the specimen. The 
larger the specimen the greater the probability that f lav/s 
will influence the strength of the specimen, since with more 
material there is a greater chance of a critical grouping of 
flaws. Weibull statistics would, therefore, predict that 
larger specimens generally have loWer nominal failure stresses. 
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Brittle materials such as composites are extremely suscepti- 
ble to flaws at the microstructure level, so this effect 
should be relatively important for composites. 

The effect of free edges on interlaminar stresses is 
expected to be another factor affecting scaling in these 
tensile tests. At a free edge there is a large concentra- 
tion of interlaminar shear stresses. Thus, there is a 
"boundary layer" along a free edge that will reach failure 
stress first, as the laminate resists the load. The width 
of the "boundary layer" will be proportional primarily to 
the ply thickness. In a specimen with a smaller cross sec- 
tional width but equal laminate thickness, the "boundary 
layer" will be a greater percentage of the cross sectional 
area. The strength of the smaller specimen would, there- 
fore, be reduced proportionally more than that of the larger 
specimen. 

In any structure where edge-effects on strength are im- 
portant, the full-sized structure will, in general, be de- 
signed with a ply stacking sequence that has a maximum mixing 
of the plies at different angles. Since prepreg composite 
materials generally aren't available in thicknesses less 
than those used for full-sized structures, ply thickness 
can't be scaled when making small-scale structures. 

No other inherent scaling mecahnism has been identified 
in this study on tensile loading. There can always be dif- 
ferences in manufacturing and testing procedures that result 



185 


from specimen size differences. Every effort vri.ll be taken 
to minimize these effects in this test program. Presumably, 
similar care can be exercised when scale model tests are 
performed on real structures. 

3. Progress During Report Period 
For an initial test, dog-bone tensile specimens were 
made from graphite-epoxy with material properties as listed in 
Table IV-A-2-1. Specimens were constructed with three dif- 
ferent ply lay-ups: [0/30/-30] , [0/45/-45] , and [0/60/ 

06 2 0* S 

-60] . These lay-ups were chosen because they seemed typi- 

cal of practical lay-ups used in real structures. The dif- 
ferent sized specimens were scaled in all three dimensions. 

The largest specimen sizes were limited by material costs 
and testing machine capacities. The smallest specimen sizes 
were limited by testing machine sensitivities and grip geom- 
etries. Also, with the above-listed lay-ups, the minimum 
thickness is that corresponding to six prepreg layers. Ini- 
tially, specimens were made in three different sizes with 
a = 1, 2 and 4 (6, 12 and 24 plies). A fourth, larger sized 
specimen was initially planned but was deleted from this be- 
ginning set of tests to save material costs. Larger speci- 
mens can be included in later test programs, once the reli- 
ability of the experimental procedure has been established. 
Three specimens were made for each size of each lay-up for a 
total of 27 specimens. 
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Figure IV-B-1 gives the dimensions for each of the 
specimen sizes an sequence. These dimensions provide for a 
safety factor of two for the shear stress incurred in the 
grips and a relatively gentle, five to one angle ratio an 
going from net section wadth to grip area width. Both of 
these proportions were judged to be sufficiently conserva- 
tive to provide reliable tests results. 

During testing it was discovered that the two smaller 
sizes could not be tested on the large Instron testing ma- 
chine. A smaller Instron machine more suitable for testing 
the smaller sizes was found. It was not effective for test- 
ing the smallest size, however, because the size of the 
teeth on the grips caused too much damage to the specimen 
grip area, which was very thin (.030 inch). Results were 
thus obtained only for the middle two sizes (12 and 24 plies) . 

Calibration of the two Instron testing machines was 
necessary to insure accurate comparisons between the two 
specimen sizes. Calibration was performed with a proving 
ring. The smaller machine was accurately calibrated but the 
larger one was slightly off. To correct for the error in 
the large machine, corrections were applxed to the measured 
loads to obtaxn "true loads". These corrections were a func- 
tion of load and varied from -2% to +1%. Tables IV-B-la 
and IV-B-lb contain the experimental results for the small 
and large specimens, respectxvely . 
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.31 

.62 

1.24 



NOTE : 6 

plies: 

.03 in 
thick 

12 

plies: 

.06 m 
thick 

24 

plies : 

. 12 m 
thick 


.125 nominal 
.250 nominal 
.500 nominal 



Figure IV-B-1 

Dimensions for Three Specimen Sizes 
(all dimensions m inches) 




188 


TABLE IV-B-la 

FAILURE DATA FOR THE 12 -PLY SPECIMENS 


Ply 

Angle 

(degrees) 

Specimen 

Number 

Neck 
Width 
(in. ) 

* 

Area 

( sq . in . ) 

Failure 

Load 

(lb.) 

Failure 

Stress 

(ksi) 

30 

1 

— 

- (damaged while testing) 


2 

0.285 

0.0171 

1987.0 

116.2 


3 

0.288 

0.01728 

2563.0 

148.3 


average 

0.2865 

0.01719 


132.26 

45 

1 

0.254 

0.01524 

1400.0 

91.9 


2 

0.248 

0.01488 

1425.0 

95.8 


3 

0.253 

0.01518 

1600.0 

105.4 


average 

0.2517 

0.01510 


97.68 

60 

1 

0.248 

0.01488 

1300.0 

87.4 


2 

0.253 

0.01518 

1325.0 

87.3 


3 

0.250 

0.01500 

1350.0 

90.0 


average 

0.2503 

0.01502 


88.22 


Ply 

Angle 

(degrees) 

FAILURE 

TABLE IV-B-lb 

DATA FOR THE 2 4 -PLY SPECIMENS 

Specimen 

Number 

Neck 
Width 
(in. ) 

* 

Area 

(sq. in.) 

Failure 

Load 

(lb.) 

Failure 

Stress 

(ksi) 

30 

1 


-(damaged while testing) 


2 

0.53 

0.0636 

7758.0 

122.0 


3 

0.52 

0.0624 

7946.0 

127.3 


average 

0.525 

0.0630 


124.7 

45 

1 

0.47 

0.0564 

5281.0 

93.6 


2 

0.50 

0.0600 

5481.0 

91.5 


3 

0.58 

0.0696 

6233.0 

89.6 


average 

0.517 

0.0620 


91.56 

60 

1 

0.50 

0.0600 

4450.0 

74.2 


2 

0.50 

0.0600 

5123.0 

84.4 


3 

0.53 

0.0636 

5385.0 

84.7 


average 

0.51 

0.0612 


81.41 


* 

Test region cross section 



a. Free Edge Effects 


The magnitude of the interlaminar shear stress near 
the free edge of a laminate can grow to levels which cause 
large plastic deformations well before loads which can cause 
failures due to normal stresses. An inelastic solution is 
not readily available at this time, however, so an elastic 
solution was used to obtain an estimate of the effects of 
these free edge stresses. 

rn * 

Tang and Levy give an elastic solution for T , 

zx 

the shear stress oriented parallel to the free edge: 


T 


zx 


E° 

x 


Q 


Q 12 Q 26 


16 


Q22 


(tt/2) 


sin h(frg/2) 
cos h{7r£ b /2) 


sin (ttti/2 ) (1) 


where: 



and 

Q. . are the off-axis elastic moduli. 

This equation can be used to estimate the size of the bound- 
ary layer (i.e., the width of the region near the free edge 
of the laminate) that doesn't contribute to resisting axial 
forces. One measure of this might be the width for which 
T zx exceeds 5 ksi. It is assumed that this is the shear 
stress at which plastic strains become significant. The 
second columns of Tables IV-B-2a and -2b give these values. 


* 

. Numbers in brackets in this section refer to the references 
which are listed on page 197 . 



TABLE IV-B-2a 


WIDTH OF BOUNDARY LAYER FOR 12-PLY SPECIMENS 


Ply 

Angle 

(degrees) 

Width for 
T 

zx 

Less Than 
5 ksi 
(in.) 

Estimated 
Boundary 
Layer 
Width 
(in. ) 

30 

.00375 

.01125 

45 

.00239 

.00717 

60 

.00078 

.00234 


TABLE IV-B-2b 

WIDTH OF BOUNDARY LAYER FOR 24 

-PLY SPECIMENS 

Ply 

Angle 

(degrees) 

Width for 
T 

zx 

Less Than 
5 ksi) 
(in. ) 

Estimated 
Boundary 
Layer 
Width 
(in. ) 

30 

.0037 

.0111 

45 

.0024 

.0072 

60 

.00067 

. 00201 
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Since the large shear stresses predicted by the elastic 
analysis would never be achieved because of large plastic 
deformations, the actual boundary layer would be larger than 
these values. Lacking a more reliable non-linear analysis, 
physical intuition must be relied on to estimate the size 
of this increased boundary layer. For the time being, a mul- 
tiplying factor of three has been chosen to approximate the 
size of the actual boundary layer. This estimate is con- 
tained in column three of Tables IV-B-2a and -2b. A re- 
duced area (A^) can now be calculated as 


A r = t(w - 2B Z ) (2) 

where t is the laminate thickness, w is the average of the 
measured test section widths and is the boundary layer 
width. The nominal area is 

A = tw (3) 

These free edge effects predict a ratio of failure 
stress due to scaling, of the small to large specimens equal 
to 


(A r /A) s 
CA r /A) ~ 


(4) 


Table IV-B-3 shows the predicted scaling effects due to 


free edge effects in terms of a ratio of the failure stress- 
es. As can be seen, smaller specimens show a failure stress 
reduction of between one percent and four percent. 
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TABLE IV-B-3 

PREDICTED FREE EDGE SCALING EFFECTS 
Small Specimens Large Specxmens 


Ply 

(degrees) 

A r 

(sq. in.) 

A 

r 

A 

A 

r 

(sq. in.) 

A 

r 

A 

a 

s 

30 

.01584 

.9215 

.06034 

.9577 

.9622 

45 

. 01424 

.9430 

.06031 

.9721 

.9701 

60 

.01474 

.9813 

.06072 

.9921 

.9891 


TABLE IV-B-4 

PREDICTED WEIBULL SCALING EFFECTS 


Ply 

Angle 

M* 

Ms 

a 

s 

(degrees) 

a £ 

30 

7.62 

1.114 

45 

8.47 

1.120 


60 


8.24 


1.119 
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b. Weibull Statistical Effects 

For brittle materials such as graphite-epoxy composites, 
the size effect on strength through flaw distribution is 
much larger than it is for more ductile materials. Since 
this size effect is due to variability of strength of the 
constituent material , the larger the scatter in these prop- 
erties, the greater the effect. A flaw in a brittle material 
causes a stress concentration and the material is unable to 
flow plastically to relieve the stresses at these concentra- 
tions as a more ductile material does. Thus, it can be pre- 
dicted that a smaller specimen will fail at a higher average 
stress since the probability of encountering a large imper- 
fection grows as the volume of the specimen increases. 

. 12 3 4 51 

Basically, Weibull Theory ' ' ' J says that a speci- 
men of volume, V, containing an even distribution of flaws 
throughout its volume has a probability of survival , S , when 
subjected to a stress, a, equal to 


S 


exp[- f v 



m 

dxdydz] 


( 5 )' 


where cr Q is the characteristic strength of the material and 
- m is the Weibull Modulus, a measure of the variability of 
the material. When a large number of specimens are tested, 
this equation reduces to 



] 


m 


S = exp[- KV 


( 6 ) 
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In this study, specimens of the same shape but different 
volumes are compared. Thus, the equation becomes 


r y 

a i 

= V 

M 

a 

2 

a 

O 


and rearranging gives 


(7) 



V_ 


V, 


1/m 


( 8 ) 


To calculate the Weibull Modulus, m. Equation (6) can 
be manipulated to become 
1 

log log — = m log a - (m log a - log KV - log log e) (9) 
s o 

. 1 

Thus, if log— is plotted versus o on log-log paper, m is 
s o 

the slope of the straight line. To find S it is seen that 
S = 1 - F, where F is the probability of failure. 

One equation used in Weibull Theory to calculate F is 


F = 


- j " 0.3 


( 10 ) 


n + 0.4 

where n is the total number of samples and j is the rank of 
the specimen . 

These tests gave a value for m of 18.78 with a correla- 
tion of .98. In previous tests, m for graphite-epoxy has 
been found to be approximately 20, which is in very close 
agreement . 


Table IV-B-4 shows the expected scaling effects due to 
Weibull flaw distributions as calculated using Equation (8) 
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and accounting for the reduced areas calculated from free 
edge effects. As can be seen, smaller specimens show a 
stress increase of about twelve percent. 

c . Correlation of Theory and Experiment 

Table IV-B-5 shows the predicted and measured scaling 
effects. The predicted values are based on the combined 
effect of free edge interlaminar stresses and Weibull sta- 
tistics as noted in the preceding sections. These results 
indicate that smaller scale specimens would be expected to 
show a failure stress increase of between five percent and 
ten percent for a scaling factor of two. The predicted and 
measured values are in reasonably good agreement. 

4 . Plans for Upcoming Period 
The postulated scaling mechanisms do appear to model 
the measured tensile loading behavior closely enough to al- 
low scaled models to be used reliably when designing large 
structures subjected to pure tensile loading. In the coming 
period, an additional series of more accurately controlled 
tensile tests over a wider range of parameters will be per- 
formed to confirm these preliminary results. Attempts will 
be made to establish a non-linear analysis to predict the 
free edge interlaminar stresses more reliably. 

Investigations of other types of structural behavior 
will be planned as regards scaling effects. It would be 
particularly important to concentrate on those structural 
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TABLE IV— B- 5 


COMPARISON 

OF PREDICTED 

TO MEASURED 

VALUES 

Ply 

Angle 

(degrees) 

Measured 

Failure 

Stress 

Ratios 

‘V 0 *’ 

Calculated 

Failure 

Stress 

Ratios 

{a s /a t) 

% 

Error 

30 

1.061 

1.072 

1.0 

45 

1.067 

1.087 

1.9 

60 

1.0835 

1.107 

2.2 
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elements that show larger effects due to scaling. A member 
of Rensselaer's Industrial Technical Advisory Committee in- 
dicated that very large effects had been found when making 
scale models of joints. We plan, therefore, to first inves- 
tigate bonded double lap joints. 

The failure data for different ply lay-ups were compared 
against several failure criteria. Strain in the fiber direc- 
tion gave much better correlation than the other popular cri- 
teria, like quadratic stress interaction. This seems reason- 
able for these types of lay-ups with a significant number of 
plies in at least three widely spaced directions, since matrix 
failure should have little effect on overall structural fail- 
ure in a single static test. It is planned to investigate 
failure criteria further. 


5. References 


1. Tang, S. and A. Levy, "A Boundary Layer Theory - Part II: 
Extension of Laminated Finite Strip", Journal of Compos- 
ite Materials, Vol. 9, January 1975, 42-52. 

2. Hitchon, J. W. and D. C. Phillips, "The Effect of Speci- 
men Size on the Strength of CFRP", Composites, April 
1978, 119-124. 

3. Hitchon, J. W. , W. H. McCausland and D. C. Phillips, "The 
Variability of Tensile Strength of Unidirectional Carbon 
Fiber Reinforced Plastics", AERE-R 8217, November 1975. 

4. Braidon, P. M. , "An Introduction to Weibull Statistics", 
AERE'-R 7165, December 1975. 

5. Weibull, W. , "A Statistical Distribution Function of Wide 
Applicability", Journal of Applied Mechanics, Vol. 18, 
1951, 293-297. 



198 




199 


park him k m? e Preceding page blank 


IV-C LOCAL BUCKLING OF COMPOSITE TUBES 

Senior Investigator: D. B. Goetschel 


1. Introduction 

Because of their high unidirectional strength, contin- 
uous fiber composite materials are frequently used in thin- 
walled structures. When these types of structures are sub- 
jected to large compressive loads, such as is encountered on 
the L-1011 wing-mounted engine drag strut, for example, de- 
signing the member to be adequate for local buckling is very 
important . 


2 . Status 

In the course of the preliminary design of a structure 
such as the L-1011 drag strut, it often will appear that the 
central portion of the strut can be relatively thin-walled. 
Local buckling, then, becomes a prime consideration. The 
best cross-sectional shape usually isn't obvious; I-beams 
or round, elliptical or rectangular tubes are all possibil- 
ities. The local buckling characteristics associated with 
these shapes could be the deciding factor m choosing among 
them. We concluded, based on results of a preliminary 
search, that the available equations are inadequate to the 
task of making that choice. To further investigate this, we 
decided to first concentrate on establishing the proper for- 
mulas for round tubes, in a form suitable for design purposes. 
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3. Progress During Report Period 

This research is the work of graduate student Shu-Ping 
Chen. The earlier literature survey was extended to make 
an exhaustive collection of equations for predicting the 
local buckling of laminated, anisotropic round tubes. The 
papers perused are listed on pages 201 through 204. 

A large number of theoretical studies were found in 
which various simplifying assumptions were made. Some in- 
vestigators made use of empirical factors derived through 
correlations with experimental results to account for ini- 
tial imperfections. 

A test program was set up in order to make an initial 
screening of these various theories. Ten specimens were made 
from prepreg graphite-epoxy. They each had a two-inch inside 
diameter, a ten-inch length and were either three or four 
layers (.0165 or .022 inches) thick. These dimensions were 
chosen to eliminate Euler buckling, inelastic material res- 
ponse and prebuckling deformation as complicating factors. 

The lay-ups used all had zero degrees and plus-minus 30, 45 
or 60 degree plies. 

The buckling loads obtained in tests of these specimens 
were only 40 to 70 percent of those predicted by the various 
theories. Such large discrepencies can hardly be accepted 
in a design method. Since the theories being evaluated in- 
clude a reduction factor to allow for. initial imperfections, 
their prediction should have provided a reasonable lower 
bound to the test results . 
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4 . Plans for Upcoming Period 
It seems clear that there must have been something 
wrong with either the tests or the predictions to have such 
large discrepancies. Perhaps the equations used for the pre- 
dictions have been empirically calibrated in a different 
region of the design space. The reduction in buckling load 
due to initial imperfections would be expected to be a func- 
tion of the thickness to diameter ratio. Another possibility 
is that the initial imperfections in the specimens used in 
previous studies may have been unreasonably small compared 
to what would be expected in production parts. Alternative- 
ly, the initial imperfections in the current tests may have 
been unreasonably large. These possibilities as well as 
the possibility that other faults existed in the experimental 
procedure used here will have to be explored further. In 
addition, we hope to begin developing design equations for 
the other candidate cross-sections. 
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IV-D ACOUSTIC EMISSION OF LOCAL BUCKLING IN GRAPHITE TUBES 

Senior Investigator: H. A. Scar ton 

1. Introduction 

The use of acoustic emission (AE) testing for detailed 
analysis of failure scenarios in advance of ultimate fail- 
ure has proven successful in delineating their mechanism of 
[1-51 * 

failure . Because of this capability, the AE technique 

is being used in the study of local buckling in slender pre- 
preg graphite-epoxy circular cross section tubes, reported 
in Part IV-C of this report. 

2, Status 

This acoustic emission study, initiated during the cur- 
rent reporting period, was primarily performed by Mr. Chi-Min 

r 6 1 

Chang 1 J , in support of the local buckling research of Mr . 

[71 

Shu-Ping Chen . The major emphasis in this effort is to 

assist in assessing the degree of local buckling occurring 

in the specimens. Examples of local and Euler buckling modes 

* * 

appear in Figure IV-D-1 

3 . Progress During Report Period 
a: Specimens and Test Set-Up 

[81 

Ten specimens made of prepreg epoxy (1048 AE prepreg L J ) 
were used for acoustic emission testing. The first specimen 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 217 . 

** 

For the sake of clarity, figures in this part are placed at 
the end of its text, pages 219 through 243. 
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was constructed using four layers of 0.559 mm (0.022 in.) 
thickness. The remaining nine specimens were in three layers 
of 0.419 mm (0.0165 in.) thickness. All ply lay-ups used 
axially oriented layers, i.e., where fibers are aligned at 
0°, along with balanced plies, with one of the following or- 
ientation pairs: ±30°, ±45° or ±60°. Table IV-D-1 shows 

the identification numbers used for the specimens. 

In the initial testing of the first two tubes, failures 
occurred at the ends of the specimens. Each of the remaining 
specimens, therefore, was reinforced with glass epoxy mate- 
rial on both of the ends, as shown in Figure IV-D-2a. Fig- 
ure IV-D-2b is a photograph of the nine three-layer specimens 
tested in this study. The one four-layer specimen, not 
shown, was tested without acoustic emission instrumentation. 

Two AE sensors of the R15 type (resonant at 150 kHz) , 

r 9 1 

Physical Acoustic Corporation , were taped directly onto 
the specimens, using an appropriately curved "shoe" and high 
vacuum grease, 50A4084, from Automation Industries, Inc.^^ 
(see Figure IV-D-3) . The AE signal was bandpass -filtered 
(100 ~ 300 kHz) with a preamplifier (PAC model 1220A) . The 
total system amplification was maintained at 98 dB and a 
constant threshold voltage of 2V (after amplification) was 
used. The test dead time was set at 2 ms. Schematic dia- 
grams of the experimental set-up are given in Figures IV-D-4 


and -5. 
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TABLE IV-D-1 

STACKING SEQUENCE OF THE NINE SPECIMENS 


Specimen 
Tube Number 

1 

2 

3A 

3B 

3C 

4A 

4B 

5 


Stacking Sequence 
Outer-To-Inner Surface 

(-30V0°/30°) t 

(30 o /0°/“30°) t 

(45°/0°/“45 o ) t 

(-45V0°/45 o ) t 

(-45 o /0°/45°) t 

(45°/0V-45°) t 

(45 o /0V“45°) t 

(-60°/0°/60 o ) t 

(60°/0°/-60 o ) t 


6 
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To assure proper alignment of specimens with the load 
axis of the MTS Tension/Torsion Servo hydraulic System^ 11 "' 
(20,000 lbs. max. tensile or compressive load, 10,000 in. -lbs. 
max. torque) during the test, an alignment attachment was 
designed to be inserted in the top and bottom of the speci- 
men. 

112 ] 

Strain gages were installed on Specimen No. 3C, as 
shown in Figure lV-D-2b. A strain gage, installed on an un- 
loaded specimen. No. 4A, was used for temperature compensa- 
tion. Instrumentation used in the test was a Data Logger 
model 2760 from the FX System Corporation 1 J . 

The MTS machine was tested for background v noise 

under loaded conditions. In this test, the PAC 3400 Sys- 
rg] 

tern was set for a total system amplification gain of 63 

dB, with a constant threshold voltage of 2 V. Noise, in the 

range 38-45 dB, was measured; its probable source is the hy- 

* 

draulic system. This may be compared to a Pentel test re- 
sult of 95-97 dB. 

To reduce the MTS machine noise reaching the sensors, 

AE isolators 1, 2 and 3 were installed at the test. Figure 
IV-D-4 shows this schematically. AE isolators 1 and 2 con- 
sisted of Vibrex Pad (1/4" thickness) produced by Vibra- 
sciences Incorporated^^ , and AE isolator 3 was an E. A. R. 
ISODAMP C-1002-12 (1/8" thickness), produced by E . A. R. 

* 

Lead break occurred using five clicks on a 0.5 mm Pentel 
pencil held at 45° adjacent to the sensor. 
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r 14 151 

Corporation ' . By this means, the machine noise was 

reduced at the sensors to the 30-40 dB range. 

b. Initial Test Results 

The computer controlled tests were conducted with mono- 

* 

tonically increasing stroke control conditions. A typical 
load stroke curve for tube number 3A is shown in Figure IV- 
D-6. It appears that the glass ends underwent load readjust- 
ment at points A and B (Note the increasing slope of the 
curves after events A and B . ) . Local buckling occurred at 
point C. Figure IV-D-7 shows strain versus time curves for 
tube 3C. (Note that the other example uses tube 3A.) Symme- 
try is in evidence during most of the test, verifying pure 
compressive loading conditions. For reasons of economy, 
tube 3C was the only specimen instrumented for strain meas- 
urement. The more interesting AE data occurred in tube 3A, 
which has a different time history. 

The AE sensors indicated AE events all through the load- 
ing sequence. Discussion of these results is divided, here, 
into the following three sections. 

i) Overall test time, 

ii) The local buckling period, 

iii) Other selected points and time periods. 

i) Overall test time . From AE displays "Parametric (Load) 
Rate (i.e.. Load Change per Division of Time) versus Time", 
as shown in Figure IV-D-8 for tube number 3A, three peaks of 
acoustic emission activity were found (point A at 40 sec.. 
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* 

point B at 51 sec. and point C at 115.5 sec.) . Similarly, 
event rate (events per tick) , energy rate and count rate 
versus time, Figures IV-D-9 through 11, respectively, show 
much activity at these times. The first two peaks (A and B) 
were associated with two end 'effects which occurred during 
the delamination between the graphite/epoxy and the glass/ 
epoxy tape reinforcement. The delamination (shown in Figure 
IV-D-2a) had been inspected after tube number 3A was tested. 

Analysis of the amplitude distributions during these 
times (see Figures IV-D-12 through -18) reveals that many 
events occurred in the range 40-60 dB. Two peaks can be 
easily found, one at 45 dB, the other at 55 dB. Failure 
mechanisms cannot be assigned with confidence on the basis 
of the overall time period data record, since a number of 
different phenomenon occurred. It is necessary to compare 

/ 

amplitude distributions obtained with tube number 3A: (1) 

before the end-effect events occurred and (2) before local 
buckling occurred. Note that an increase in events at ampli- 
tudes around 40-60 dB coincided with delamination between 

$“ 

Note that the PAC 3400 plotting system labels the horizontal 

time axis in 120 divisions or tick marks, with the numbers 

1 

being written vertically, such as 2. This represents 120 

0 

ticks x 1.5 sec. /tick = 180 seconds. Figure IV-D-8 shows 
only parametric rate data. Since an absolute calibration 
was not deemed necessary, only the temporal sequence of 
events is needed to provide insight into the actual time 
history. Rate data analysis, rather than cumulative data, 
is used, since this technique increases the sensitivity of 
the instrument to particular failure scenarios. 
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the two different materials at the ends. This is discussed 
further in Section iii, page 212. 

ii) Local buckling period . The PAC 3400 has the capa- 
bility of "labelling" data recorded within certain time per- 
iods. Data can then be recovered for analysis purposes 
from that aggregated without respect to time, by a technique 

called "time-filtering". This was done, in the Post Test 
[9] 

Mode , for each tube so as to eliminate events which oc- 
curred at times associated with postbuckling behavior. For 
this reason, the events from between time ticks 85 and 120 
(127.5 to 180 seconds) have been suppressed in the data dis- 
cussed in this section. 

It will be seen from the time-filtered display of "Event 
Rate versus Time" (Figure IV-D-19) or "Energy Rate versus 
Time" (Figure IV-D-20) that the event and energy rates are 
not at their highest levels at the final load drop point C, 
for tube number 3A. The end effects for that specimen had 
occurred before point C, and their energy rates are almost 
the same as those associated with point C. However, before 
point C, there is an event scenario associated with point D 
which has a still larger AE energy rate. In fact, when the 
precurser load drop point (point D m Figure IV-D-8) occurs, 
the material is seen to release the largest AE energy rate 
(see Figure IV-D-20) . Point D might be explained as the on- 
set of local buckling and the following load drop point C as 
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the point at which local buckling has been completed and its 

[7] 

mode shape fully developed 

Analysis of mode shapes indicates that the number of 
pulses of rise time versus time may be twice the experiment- 
ally measured circumferential mode number N. For example. 

Figure IV-D-21 shows specimen 3A after testing; a mode num- 

[71 

ber of three is in evidence. Figure IV-D-22 shows the 
identical post-test buckled shape under the action of an 
axial preload to demonstrate its deformed shape. Study of 
this shape for tube number 3 A reveals that the circumferen- 
tial mode number N is, in fact, 3. Analysis of the number 
of pulses of rise time rate versus time, as seen m Figures 
IV-D-23 and -24, between the time tick divisions 60 and 84 
(90-126 seconds) , is found to be approximately 6 or 2N (two 
times the circumferential mode number) . Note that Figure 
IV-D-24 shows only the 90-126 second time intervals. The 
same conclusion was observed by studying data from the other 
tubes. 

Plots of duration rate (Figure IV-D-25) , amplitude rate 
(Figure IV-D-26) and count rate (Figure IV-D-27) versus time 
for specimen 3A also show the specific characteristics of 
local buckling . 

iii) Other selected points and time periods . As has 
been mentioned previously, subsegments of AE data at any 
time can be filtered by the "post" and "exam" analysis of 
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[9] 

software available on the PAC 3400 . For tube number 3A, 

the segments of interest for initial end effects occur at 
points A, at time tick 30 (45 seconds) and B , at time tick 
35 (37.5 seconds); the local buckling period is a segment in 
the time tick range from 60-84 (90-126 seconds) . Figures 
IV-D-28 and -29 show the display of Energy Rate versus Time 
and a Time-Filtered Amplitude Distribution for point A alone; 
Figures IV-D-30 and -31 for B alone; and Figures IV-D-32 and 
-33 for C alone. Figures IV-D-34 and -35 show the AE display 
for the local buckling segment between 60-74 ticks (90-111 
seconds) . 

Recall that the AE amplitude range of delamination for 

composite materials such as glass/epoxy and graphite-epoxy 

* 

can be in the range of 45-60 dB , while fiber breakage occurs 
above 8 0 dB'^ . The same conclusions can be obtained from 
points B and C, as seen in Figures IV-D-30 through -33. 

Since Figure IV-D-33 shows an amplitude distribution m the 
range of 45-60 dB, the initiation of local buckling could be 
associated with debonding. 

Figures IV-D-36 and -37 portray events which occurred 
at a maximum constant stroke. Analysis of these plots veri- 
fied the visual observation that after local buckling occur- 
red, intense fiber breakage (amplitudes above 80 dB) took 
place. Figure IV-D-37 depicts events which occurred early 

* 

Attenuation with distance tests using Pentel-lead breaks 
were performed on the tube, and they showed less than a 5 dB 
drop in amplitude over the length of the specimen. This 
lends credence to arguments based on amplitude characteris- 
tics. 
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in the constant-hold stroke phase and shows only fiber break- 
age. Figure IV-D-38 shows the amplitude distribution when 
the stroke was held to maximum time; it would appear from 
this that more delamination then occurred. 

Table IV-D-2 sumarizes results from the remainder of 
the specimens studied. Analogous conclusions were drawn by 
either analysis similar to that done for tube 3A or by in- 
terrupting the test prior to one of the above effects for 
visual inspection and appropriate AE analysis. 

c. Preliminary Conclusions 

Experience with AE analysis of specimens subject to de- 
bonding suggests that the major amplitude distribution during 
incipient local buckling is m the range associated with de- 
bonding (around 45-60 dB) . Postbuckling behavior, on the 
other hand, seems associated with AE amplitude distributions 
typical of fiber breakage. It also seems possible that AE 
parameters such as energy rate, rise time rate, duration rate 
and count rate may have a numerical relation with the circum- 
ferential mode indices which describe the local buckling 
shapes of composite tubes. Taping the ends of specimens with 
the same material as they are made of promises to eliminate 
the end effect for the future tests. 

4, Plans for Upcoming Period 

A major study in unaxial tension-tension, tension-com- 
pression, biaxial tension-torsion low cycle fatigue and 
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TABLE IV-D-2 

SUMMARY OP THE LOCAL BUCKLING TEST 


Buckling 

Load 



Test Item 

(lbs . 

1 . 

Tube No. 1 
(-30°/0 o / 30°) t 1 11 111 1 

1935 

2. 

Tube No. 2 
(30°/0V-30 °) t 

1550 

3. 

Tube No. 3A 

(45 o /0°/-45°) t 

2095 

4. 

Tube No. 3B 
(-45°/0°/ 45°) T 

2605 

5. 

Tube No. 3C 
(-45°/0°/45 °) t 

2100 

6. 

Tube No . 4A 
(45 o /0°/-45 o ) t 

1690 

7. 

Tube No. 4B 
(45V0°/-45 °) t 

1950 

3. 

Tube No . 5 
(-60°/0°/6 0° ) T 

2365 

9. 

Tube No. 6 
(60V0°/-60 o ) t 

2800 


Acoustic 





Emission 

(events) 

Inspection 
After Test 

Note . 
Number 

— — 

m 

— 

liii 

(1) 


n 

= 

4 iv 


— 

m 

= 

1 

(1) 


n 

- 

4 


1330 

m 


1 

(2,3) 


n 

= 

3 


843 

m 

= 

1 

(4) 


n 

- 

4 


— 

m 

= 

1 

(5) 


n 


3 


— 

m 

= 

1 

(6) 


n 

= 

3 


1664 

m 

= 

1 

(7) 


n 

= 

3 


1246 

m 

= 

1 

(8) 


n 

= 

3. 


— 

m 

- 

1 

(9) 


n 

= 

3 



1 See following page for notes 

i;L Stacking sequence, outer-to-inner surface 

111 m = Number of axial half-waves for a displacement mode in 
a simply- supported shell 

1V n = Number of circumferential waves as seen in Figure 
4 . 2 ^ 
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TABLE IV-D-2 (Continued) 


Note Information ; 

(1) Tubes nos- 1 and 2 fractured at both ends. Therefore, 
each of the tubes was reinforced on both ends with fila- 
ment (glass/epoxy) winding before testing, and Tubes 1 
and 2 were not monitored with AE instrumentation. The 
test sequence is Nos. 3C, 6, 4A, 4B, 3A, 3B and 5. 

(2) As shown in Figure IV-D-6, load drop points A and B are 
due to end effects caused by delamination between speci- 
men and the end reinforcing material (glass/epoxy) . 

(3) Tube no. 3A has good AE results, for local buckling and 
end effects. 

(4) Tube no. 3B has good AE results, for local buckling with- 
out end effects (the load-stroke diagram for this tube 

is shown in Figure IV-D-6.). 

(5) Tube no. 3C failed rapidly under load control (500#f/min) ; 
after testing of this specimen, stroke control of the MTS 
machine was advised. 

(6) Tube no. 4A lost AE signals due to misoperation of AE 
instruments . 

(7) AE recording was stopped when local buckling noise from 
the test of tube no. 4B was heard. 

(8) Tube no. 5 has good AE results for local buckling and 
end effects. 

(9) AE recording was stopped just before the first end effect. 
No AE data for local buckling was obtained. 
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monotonic loading of cylindrical composite tubes is underway. 
We expect to be able to report results from these studies 
during the next project period. 
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Figure IV-D-1. Buckling Mode Shapes 


220 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



Figure IV-D-2a. Specimen Tube No. 3A with Reinforced Glass/Epoxy Ends 



Figure IV-D-2b. Tubular Test Specimens with Glass/Epoxy End Reinforcements 
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Figure IV-D-3 . AE Sensors, AE Shoes, Alignment Attachment for Tube No. 3A 
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Figure IV-D-4. Schematic Diagram of Test Set-Up 
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Figure IV-D-5 . Schematic of the Set-Up for the Acoustic Emission System 









OUAUlY 


Where 

AjB = Load drop points due to end effect 

C = Local buckling load 

A = 1025 lbs (corresponding to AE 
417 events) 

B = 1862 lbs (corresponding to AE 
613 events) 


S = Stop MTS system and inspect 
the tube for a short period 




Strain (10 in. /in 
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Figure IV-D-7 . Curves of Strain vs. Time of Tube No. 3C 


Load Change/Tick 
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4 ft, ft-tap is ON LOCAL BUCKLING OF 45/9M5 TUBE 3A 23:55*37 04/28/82 



Tiss per Tick tX 180 as. 3 ■ 8915 

First Arrival = ft Event I 92116 

Mi-. Replay of Paraseiric vs. Ties i flax = 838128 felts s 8.828 
Stanford Display*' Data replay existed 

Rate Data 


Figure IV-D-8. Parametric (Load) Change/Tick vs. Time for Tube No. 3A 
(Note the pause in stroke ramp at the 36th to the 50th 
tick mark for visual inspection of the specimen) 


4 ft, fi-Duap is ON LOCAL BUCKLING OF 454H5 TUBE 3 h 64:39:47 04/26/82 



Tiers per Tick E X 188 us, 3 s 0815 

First Arrival = ft Event § 82116 

ffiDE: Replay of Events us, Tits? fox = mm Volts “ 8,898 

Stanford Display Data replay completed 

Rate Data 


Figure IV-D-9. Event Rate vs. Time of the Whole Test for Tube No. 3A 
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4 Ch, Mm is 08’ L0C8L RICKLIHfi OF 45/B/-45 TUBE 3ft @&.42'5? 04/26/82 
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Rate Data 



Figure IV-D-10. Energy Rate vs. Time of the Whole Test for Tube No. 3A 


4 tt, Wkap is OH LOCH BUMS OF 45/&H5 TUBE 3A 86=53: 15 . 84/26/82 



Figure IV-D-11. Count Rate vs. Time of the Whole Test for Tube No. 3 A 
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4 CK, <H5us? is OH LOCAL BUCKLING OF 45/fr'-45.TlfflE 3fi 63-31 46 04/33/82 



168 db. Full Seals £ Cat A= m PE Cnt 8= 8891 
First Arrival = B Loc = 088844 Event I 68413 
HCGE: Events vs. ftspliluds Ha? = $1889 Mts= 8.144 


Figure IV-D-12. Amplitude Distribution for Tube No. 3A 
Before End Effect Point A 


4 Ch, iHkirp is OH LOCAL BUCKLING OF C4M5 T1EE St 03-32-1J 64/39/32 



Figure IV-D-13. Amplitude Distribution for Tube No. 3A 
After End Effect Point A 
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4 Ch, is OH LOCAL 80CKLIHG OF 4S/8H5 TUBE 3ft 83.3348 84/38/82 



8 1 2 3 4-5 6 7 8 91 

Scale 3888808808 

S ^ * Amplitude (dB re lyV) ® 


103 db. Full Scale 
First Arrival = ft 
M&-- Events vs, Amplitude 


AE Cnt fc 8813 AS Cnt B= 8033 
Loc = M2 Event I 88SB5 
Hax = 801888 Mis* 8.196 


Figure IV-D-14. Amplitude Distribution for Tube No. 3A 
Before End Effect Point B 


LOCAL BUCKL1HS OF 45/8M5 TUBE 3ft tt&ti mm 


,? 1 \ ! ? « 1 8 5 1 


Scale 8 0 0 0 6 0 0 

, Amplitude (dB re lyv) 


183 db. Full Scale 
First Arrival = B 
BODE; Events vs. ftaplitude 


AE Cnt Pr 0803 ftTCtitls 


Event I 03617 
Holts = 6.202 


Figure IV-D-15. Amplitude Distribution for Tube No. 3A 
After End Effect Point B 
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4 Oil Mop is SI LOCftt BUCKLING OF 45/8M5 TLSE 3ft 63=10-37 94/29/82 



31234557891 
Seal? 0886008800 
s v 4 8 

Amplitude {dB re.ljjV) 


193 db. Full Scale ftECntft=8977 AE Qtt B='§982 

First Arrival s B Loc = 00888$ ftent H 81310 

RODE; Events vs. Aaplitude M - 801689 tbits = 8.345 

Figure IV-D-16. Amplitude Distribution for Tube No. 3A 
Before Local Buckling - Point C 


4 Ch, Mup is ON LOOfiL HOIK OF 45/0M5 TUBE 1 86.88:45 94/29/82 



Scale 8908088800 
-U v 0 

Amplitude (dB re lpV) 

m db. Full Scale AE Cr* ft= 9889 A£ Cut B=' 080 j' 

First Arrival = B Loc = 8S0852 Event it 01349 

t®DE: Events vs. fteplitude Hax = 591693 tbits =8,353 

Figure IV-D-17 . Amplitude Distribution for Tube No. 3A 
After Local Buckling - Point C 
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4 Ch, ft-CuEP is OH LOCAL BUM OF 45/8M5 TUBE 3A 04; 12=33 M/26/82 
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Figure IV-D-18. amplitude Distribution for Tube No. 3A 
for the Whole Test 
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Standard Display ' 

Rate Data 

Figure IV-D-19. Event Rate vs. Time for Times Associated with 
Local Buckling of Tube No. 3A 
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Figure IV-D-20. Energy Rate vs. Time of Local Buckling 
for Tube No. 3A 
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Figure IV-D 21. Posttest Photograph of Buckled Specimen 3A 
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Figure IV-D-22. Additional Posttest Photograph of Buckled Specimen 

. [7] 

3A Under Posttest Axial Preload Condition 
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4 Ch, (Htap is OH LOCAL BUCKLING Of 45/0M5 TUBE 3d 03:34:21 94/29/82 



Ti Mper Tide [ X 196 is.] r'MB 
First Arrival = B 
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12 
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Figure IV-D-23 . 


Rise Time Rate vs. Time of Local Buckling 
for Tube No. 3 A 


4 Ch, fHtop is OK LOCAL BUCKLING OF 4MHS TUBE 3A 06:13:54 04/30/82 



Figure IV-D-24. Rise Time Rate vs. Time of Local Buckling 

for Tube No. 3A (An amplified or rescaled time) 
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4 Ch, (Hluq> is m tOCftL SMS OF 4MM1 TUBE 3ft 33:37:38- WML 
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Figure IV-D-25. Duration Rate vs. Time of Local Buckling 
for Tube No. 3A 
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Figure IV-D-27 . Count Rate vs. Time of Local Buckling 
for Tube No. 3A 
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Figure IV-D-28. Energy Rate vs. Time of End Effect 
Only at Point A for Tube No. 3A 
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4 Ch, ft-Ouep is ftl LOCAL BUCKLIKG OF 45W-45 TUBE 3A 15:37:14 84/25^2 
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Figure IV-D-29. Amplitude Distribution of End Effect 
Only at Point A for Tube No . 3A 
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Figure Iv-D-30. Energy Rate vs. Time of End Effect 
Only at Point B for Tube No. 3A 
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Figure IV-D-31. Amplitude Distribution of End Effect 
Only at Point B for Tube No. 3A 
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Event § 0D016 
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Figure IV-D-32. Energy Rate vs. Time of Local Buckling 
Only at Point C for Tube 3A 
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4 ft, Muap 3? M LOCAL BUCKLING OF 45^-45 TUSE 3A 28-22--34 84/29/82 
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Figure IV-D-33 . Amplitude Distribution of Local Buckling 
Only at End Point C for Tube No. 3A 


4 Cn, fHhfip is ON LOCAL BUCKLIffi OF 45/0/-4S TUBE 1 06-M.35 04/1/82 
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Figure IV-D-34. Energy Rate vs. Time for Only the Local 
Buckling Period for Tube No. 3A 
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4 Ch, flrftip is OH LOCAL BUCKLING OF 45'6M5 ME 1 35.63=22 64/30/82 



Figure IV-D-35. Amplitude Distribution for Only the Local 
Buckling Period for Tube No. 3A 
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Figure IV-D-36. Energy Rate vs. Time in the Postlocal 
Buckling Period for Tube No. 3A 
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Figure IV-D-37. Amplitude Distribution m Early Portion of Postlocal 
Buckling Phase for Tube No. 3A 
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Figure IV-D-38. Amplitude Distribution for Total Period of 
Postlocal Buckling Phase for Tube No. 3A 




245 


Preceding page blank 


PART V 

PROCESSING SCIENCE AND TECHNOLOGY 

V-A INITIAL SAILPLANE PROJECT: RP-1 

V-B SECOND SAILPLANE PROJECT: RP-2 
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INTRODUCTION 


The sailplane project involving design, fabrication and 
testing of advanced, full-scale, f light -worthy , all-composite 
aircraft was conceived as a highly productive way to conduct 
research on processing science and technology for low-cost 
aerospace structures. A series of innovative designs and 
processes has resulted from that research to date. It in- 
cludes: 

a) A simplified fabrication process without molds 
for accurate, smooth surface airfoils in wings 
and tail surfaces up to roughly twenty feet 
semispan, 

b) Single sided, composite-foam "sandwich" struc- 
tures which can be successfully joined to 
sub- structures to form ultralight-weight sur- 
faces capable of carrying useful distributed 
loads , 

c) Composite hinge configurations that are integ- 
rally bonded into airfoil contours to minimize 
peel-strength problems, 

d) Proper mixtures of hardeners and experimental, 
especially low viscosity resins to achieve 
both minimum-depth bond and surfaces for 
ultralight - weight and acceptable strength, 
and 

e) A technique for fabricating smooth- surfaced 
graphite/epoxy tubes of relatively small dia- 
meter, without the use of release agents or 
external molds, through the use of Teflon 
rods as mandrels. 
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The impetus toward innovation in processing science is 
being preserved in this project by undertaking continually 
more challenging structural functions in sailplane design . 
Further, research on problems within the disciplined context 
of actual manned flight testing is contributing to funda- 
mental knowledge as regards monitoring and maintaining struc- 
tural integrity of functioning, low-cost aerospace structures. 
These aspects are described in the following sections dealing 
with the first and second sailplane projects. 
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V-A INITIAL SAILPLANE PROJECT: THE RP-1 

Senior Investigators: F. Bundy 

R. J. Diefendorf 
H. Hager up 

Flight tests through the summer of 1981 produced per- 
formance data for the 1981 configuration of the RP-1 showing 
best glide ratio equal to 11.3 and a minimum sink rate of 
4.2 fps. This performance was sufficiently below design 
values that a drag clean-up program, involving light-weight 
secondary structure, was initiated. A two and one-half foot 
extension was added to the fuselage to improve airflow in 
the rear, between the two tail booms. The original design 
concept included a cockpit which was open below the pilot to 
permit operation in the foot-launch mode. This area has now 
been closed to insure that the induced drag of the wing cor- 
responds to the full span, rather than half span. Fairings 
have also been provided around the supporting structure for 
the skid and around the wheel. We expect a twenty percent 
improvement in glide ratio to result from these modifications 
to the RP-1. 

Before the RP-1 aircraft flies again, we are committed 
to performing the yearly proof-test to ensure that no struc- 
tural degradation has occurred. This is anticipated to in- 
clude both a repeat of past static tests and an additional 
test intended to investigate the failure condition on the 
RP-2 aircraft (see the following section, V-B) as it may 
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pertain to the RP-1. In addition, the first indication of 
possible strength degradation of typical structural elements 
used on the RP-1, and hidden from view, will be possible af- 
ter the next reporting period. Five identical elements of 
this kind were tested in the Fall of '81, and five more will 
be tested prior to September 20, 1982. 

Resumption of flight tests is scheduled for July, 1982. 
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V-B SECOND SAILPLANE PROJECT; THE RP-2 

Senior Investigators: F. Bundy 

R. J. Diefendorf 
H . Hagerup 
C . Muser 

The completion of the structural framework of the fuse- 
lage and the finished empennage and wings of the RP-2 all- 
composite, high-performance, ultralight sailplane early last 

t 

fall was followed by structural proof -tests of the wings and 
fuselage-wing connection. The fuselage was put upside down 
on an elevated platform. Once the wings were assembled on 
the fuselage, sandbags were used to load the structure. De- 
flections were measured at several stations along the wings. 
Preliminary loadings revealed -the need for another transverse 
compression member in the fuselage framework at the forward 
part of the wing root rib. After this was installed, the 
structure was tested for the highest vertical bending to be 
encountered within the flight envelope. At that point, the 
maximum positive design load factor was 5.9. The sailplane 
successfully withstood this test without incident. 

One week later the plane was tested for vertical bending 
and torsion combined. The magnitude of the load was the same 
as in the previous test, corresponding to the maximum design 
load point on the flight envelope, but the loading line was 
closer to the trailing edge to induce torsion. In this com- 
bined bending and torsion test, the center box of the wing 
connection failed catastrophically at 98 percent of design 
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load, causing serious damage to the wing connection itself, 
and also serious secondary damage to the outboard one-third 
of the wings, which were inadequately supported when the 
wing carry-through structure failed. 

Immediately after this sobering test, a small task force 
was formed to investigate the possible sources for the unex- 
pected failure of the structure. The failure was obviously 
not due to straightforward causes . The carry-through struc- 
ture and associated fuselage design was devised so that the ' 
former would carry only bending in the vertical plane . The 
bending test was completed successfully, yet that same struc- 
ture failed during a new test condition in which torsion was 
added. The task force consisted of three students; Keith 
Applewhite, Andrew Bertolazzi and Thomas de Mint. They were 
supervised and guided by Postdoctoral Associate Christoph 
Muser. The results of this investigatory group’s efforts 
were summarized in a report submitted along with notes of 
all the calculations necessary to determine the cause of the 
structural failure. 

The conclusion reached by the task force was that the 
cause of failure was the combined action of three separate 
loadings: vertical bending, torsion and horizontal bending. 

Horizontal bending moments were induced by a nose-down atti- 
tude of the sailplane on its platform. The severe implica- 
tions of these combined loadings was not recognized earlier, 
because the load paths in the wing carry-through structure 
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were obscured by deflections whose magnitude had been under- 
estimated. In other words, a load applied on the wings — 
despite the intention of the designers to the contrary — 
could take different paths through the fuselage depending on 
deformations and on how "tight" various fittings had been 
made. 

Work in the fabrication laboratory pertaining to the 
RP-2 during the remainder of the present reporting period 
has been mainly concerned with (1) development of structural 
repair techniques, (2) wing connection failure analysis and 
redesign and (3) reconstruction of the damaged parts of the 
RP-2 . 

After evaluating all the damage done to the aircraft, 
it was decided to use all salvagable parts and to just re- 
place the broken sections, rather than to rebuild a completely 
new set of wings. This had two advantages. First, costs 
could be kept at a minimum, and second, and perhaps most im- 
portant, ther.e would be the new challenge of repairing the 
primary structure of a light-weight composite wing. 

The major areas which needed close attention could be 
singled out. The first one was the approach to splicing the 
capstrips of the I— beams used to carry vertical bending mo- 
ments, on the one hand, and splicing the wing skins which 
carry horizontal bending and torsion, on the other hand. The 
second area was the wing-fuselage connection and carry-through 
structure which had to be redesigned to eliminate the design 
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faults revealed in the test-to-failure, within the limits 
imposed by the desire to salvage parts . 

The capstrips of the I-beam are made of unidirectional 
graphite prepreg and are tapered towards the wing tip. At the 
splice locations, the capstrips are eight layers thick (0.04 
in.) , and the taper of the salvaged part of the capstrip re- 
stricted the allowable length of the splices to 40 cm (see 
Figure V-B-l) . After studying the literature (see References 
on page 263) , it was decided that the inboard and the newly 
built outboard section of the capstrips should have blunt 
ends which butt up against each other. On the top and bottom 
sides of the capstrips positioned in this way, reinforcements 
of the same material and with the same thickness as at the cen- 
terline of the capstrips should be used. The total length, 
the step-lengths and the manufacturing methods were arrived 
at in tests guided by approximate calculations. About twenty 
splice specimens were subjected to simple tension tests. 

Some specimens had reinforcements only 12 cm long, and some 
had nonuniform step-lengths. Some had step-ends which were 
arrowhead shaped in planform, others were straight and normal 
to the capstrip Cj~. . The splice designs selected for the ex- 
tensive tests required by use in a flight vehicle had rein- 
forcements 32 cm long, with straight ends , and a constant step- 
length of 4 cm. These specimens withstood a tensile load 
five times greater than the sparcap load present at the splice 
location under a design load factor of 5.9 g. It is important 



Legend - I - Salvaged part of capstrip 

2 - Reinforcements 

3 - Rebuilt part of capstrip 

4 - Room temperature cured glue line 

(also untapered length of spar cap) 


AM material 5 unidirectional graphite prepreg. 
(Drawing not in proportion) 
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Figure V-B-l. Schematic of Spar Capstrip Splice 
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to recall that the governing design criteria for the wings 
is stiffness and not strength. In this regard, the spliced 
spar will actually be stiff er than the original. 

An I-beam section of about 3 m length was built with 
spliced upper and lower capstrips. That beam was subjected 
to three times the bending moment experienced by the wings 
at 5.9 g, and no failure occurred. At this point, we are 
confident enough to plan to use this splicing technique — 
which joins pressure vessel-cured parts together with a 
room-temperature-cured bond — for the repair of the wings. 

The failed carry-through structure design had several 
shortcomings. In principle, the two I-beams of the two wings 
were connected by means of a tight-fit sleeve. That is, one 
I-beam slid into a box-end built into the other I-beam (see 
Figure V-B-2) . These two beams were prevented fr om d isen- 
gaging, as a result of spanwise motion, by a safety pin. 

Short vertical shear pins or "studs", mounted in the root 
ribs, transferred both lift forces and torsion into the tor- 
que box of the fuselage. Horizontal bending moments caused 
by the components of lift and drag m the plane of the wing 
chord, were to be taken through these same pins, into the 
fuselage sides, and resisted in compression by the fuselage 
torque box at the aft set and by a compression strut between 
the front two torsion pins. Whether this horizontal moment 
was caused by forward or aft acting loads, the safety pin 
holding the wings together would always provide the tensile 




SIDE VIEW 
(CROSS SECTION) 


Starboard wing 


PLAN VIEW 


Port wing 


U Legend : I - Torsion pin 

2 - Inserted I-beam 
3 - Wing connection box 
4 - Safety pin 
5 - Fuselage torsion box 
6 - Compression strut 
7 - Fuselage "tusks" 

Figure V-B-2. Previous Design of Wing Carry-Through Structure 
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reaction. However, the elasticity of the fuselage and the 
safety pin made it possible for part of that horizontal mo- 
ment to be transferred through the wing vertical bending 
carry-through structure's connection box. This additional 
stress concentration was most dangerous at point A of the 
open end of the box (see Figure V-B-3) . 

The redesign of the carry-through structure, conducted 
m this report period, took into account the fact that the 
capstrip sections to be used for the wing connection were of 
a given length and width and had only unidirectional fibers. 
There were similar salvage restrictions associated with the 
side of the fuselage and the positions of the torsion pins. 

The new wing connection design (Figure V-B-4) uses two steel 
bolts to transmit the vertical bending moment, which is the 
principal loading. These bolts connect the webs of the box 
of the one I-beam to the web of the other I-beam. This wing 
connection, however, is also designed to withstand the ten- 
sion loads experienced from the horizontal moment. The same 
loads, but in compression, have to be taken by the compres- 
sion strut between the front torsion pins. A careful redesign 
of that particular load path has been made. There is no 
longer any direct contact between the box at the end of one 
I-beam and the other I-beam. It is intended that the clear- 
ance at all points other than the two-bolt attachment will 
be sufficient that no unwanted loads will be carried by the 
spar connections, even in the presence of elastic deflections. 



Legend ! I - Wing connection box 

2 - Graphite cloth wrap 

3 - Capstrip of box 

4 - Web of box 

5 - Capstrip of I-beam 

6 - Web of I-beam 

M V " Vertical bending moment 
M H -Horizontal bending moment 
A - Point of highest stress concentration 


Figure V-B-3. 
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TOP VIEW 



SIDE VIEW 



w i:^==«fss=cEa-- 

4 - Steel bushings 

5 - Steel washers 

6 - Steel pins 

7 - Root ribs -G r /E (fabric) 

8 - Capstrip of w i ng connection box-6 r /E (prepreg. onidirect.onal tape) 


Design of Wing Carry-Through Structure 


Figure V-B-4. Modified 
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Flat tensile specimens were tested to determine the 
number of layers required by the spar webs to transmit the 
loads transmitted to them by the bolts. These tests demon- 
strated that it is possible to transmit loads three times 
greater than the loads experience at 5.9 g and that failures, 
when they do occur, are in bearing. Steel bushings are used 
to prevent damage during assembly and disassembly. Steel 
washers are mounted on either side of the webs (not used for 
the tensile specimens) as an additional protective measure. 

Presently, the two principal activities are repair of 
the wings and fabrication of the test spar for the new wing 
connection. The status of reconstruction is that the dam- 
aged outboard portions of the wings have been removed; the 
necessary replacement capstrips, shearwebs, stringers, ribs 
and panels have been fabricated; and assembly has begun. 

Local damage in one wing root, where the mainspar capstrips 
separated from the web for a distance of one and one-half 
inches into the left wing, has been excised and repaired, 
and a new root rib is being installed. Redesign of the 
center box of the wing connection also necessitated minor 
modification in the fuselage framework, which is in progress. 
Work is continuing, meanwhile, on the aerodynamic fairing of 
the fuselage framework and wing intersection, using RPI's 
graphics computer-aided -design capability. 

We plan to have the assembled RP-2 aircraft retested in 
October, 1982. After that, the final assembly of the fuselage 
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can be carried out and the cockpit enclosure built. 

In summary, a failure during structural proof -test has 
cost the RP-2 project six to eight months of time. An im- 
portant and positive outcome, however, is that after success- 
ful flight operations with the RP-1 for two summers without 
significant structural damage, we have been required through 
this temporary setback on the RP-2 to concentrate on the im- 
portant area of advanced composites structural repairs . The 
first flight of the RP-2 is anticipated now for the Summer 
of '83. 
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VI COMPUTER SOFTWARE DEVELOPMENTS 

Senior Investigator: M. S. Shephard 

I. Introduction 

The objective of this portion of the project is to pro- 
vide advanced and specialized computer "tools" for the anal- 
ysis and design of composite structural elements and the study 
of structures-related phenomena in composite materials. The 
major thrust to date has been and will continue to be in the 
finite element area, with effort directed at implementing 
finite element analysis capabilities and developing inter- 
active graphic preprocessing and postprocessing capabilities. 

2. Status 

The results achieved in this area of research during re- 
cent reporting periods have been used for the numerical anal- 
ysis of moisture effects and the numerical investigation of 
the micromechanics of composite fracture discussed in Parts 
III— D and -E of this report. In addition, the POFES system, 
which is available to the composites program, continues to be 
developed through student class projects and projects supported 
by the Graphics Center's Industrial Associates Group. 

3. Progress During Report Period 

Class projects with potential usefulness for the compos- 
ites program include both specific analysis studies and 

* | ' 

POFES (People Oriented Finite Element Software) is RPI ' s 
finite element software system. 
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improvements to our general in-house analysis capabilities. 
Class composite analysis studies included: 

1. "A Finite Element Analysis of a Graphite/Epoxy 
Aircraft Engine Drag Strut Lug", R. E. Walters. 

2. "Composite Ring Stress Analysis", R. K. Phelan. 

3. "The Natural Frequency and Mode Shape Investi- 
gation of an Octagonal, Orthotropic Sandwich 
Plate with the Finite Element Model Analysis", 

R. A. Valicenti. 

4. "Finite Element Analysis of L-1011 Engine Drag 
Strut", F. Ogunlau. 

Class Improvements to our general in-house analysis 
capabilities include: 

1. an eigenvalue and eigen vector solution algo- 
rithm for natural frequency and mode shapes . 

This algorithm is based on the subspace iter- 
m* 

at ion method 

2 . a new triangular plate bending element based 

r 2 31 

on discrete Kirhoff assumptions ' . This 

element should prove very efficient for the 
analysis of composite plates and shells. 

Although these class projects are limited in scope, they 
do serve to give the students involved additional insight 
into the numerical analysis of composites, they provide ad- 
ditional input into other portions of the project and/or they 
help improve our in-house facilities for the analysis of com- 
posite materials and structures. 


Numbers in brackets in this section refer to the references 
which are listed on page 277 , 
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Projects of interest to the composites program supported 
out of the graphics center or elsewhere are discussed m the 
following subsections. 

a, Three-Dimensional Preprocessor, Based on Extrusion Capa- 
bilities 

The extrusion preprocessor allows for the definition 
and meshing of three-dimensional objects by first defining 
components with planar faces that are then extruded in the 
third dimension and patched to other, similarly-defined com- 
ponents. At this time only linear extrusion capabilities 
and some initial mesh generation capabilities have been im- 
plemented. The process of defining and meshing an object 
begins by establishing a three-dimensional planning grid 
(Figure VI-1) that delimits the object space. This grid is 
also used to help maintain numerical accuracy during curve 
definition. Next, the user defines planar faces and extrudes 
them to form components (Figure VI-2) . Various extruded com- 
ponents are then interactively patched together (Figure VI-3) 
and passed on to the mesh generator where three-dimensional 
mesh patches are interactively created (Figure VI-4) . 

b. Shell Surface Program Redesign to Include Improved Model- 
ing Capabilities 

The shell preprocessor is being redesigned to allow for 
the improved representation of a wider class of surfaces and 
to improve the user interface. The basic curve type has been 
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Figure VI-4. Three-Dimensional Component Mesh 
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changed to rational parametric cubics, which allow for error- 

r 4 ] 

free representation of all conic sections . The use of 
rational parametric cubic surfaces is also planned. The 
user interface is being improved by the addition of new 
control/display features and the formalization of a menu 
hierarchy. 

c. Three-Dimensional Postprocessor Development 

An initial version of a three-dimensional postprocessor 
for deflected shapes has been developed. The program is 
designed for both static and dynamic analysis results, with 
the dynamic analysis results including both mode shapes and 
time-histories. Basic interactive operations include single 
and multiple view options, scale change, rotation, trans- 
lation, windowing and sectioning. Figure VI-5 shows a 
multiple view display of a shell surface that has been sec- 
tioned by user-defined cutting planes. 

d. Automatic Mesh Generator 

A new approach to the automatic generation of two- 
dimensional finite element meshes has been developed. The 
method is referred to as the modified quadtree approach 
[ 5 , 6 ] and appears to be extendable to three-dimensional 
problems through a modified octree method. The process 
begins with the interactive definition of the problem geom- 



Figure VI-5. Multiple View Option in Three 
Dimensional Postprocessor 



Figure VI -6. Boundary Curves for Plate with Hole 
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etry in terms of boundary curves. Figure VI-6 shows the 
curves used to define a plate with a hole in it. Curve 
tolerances are then defined interactively to produce 

'•A 

desired mesh gradations. From this point the mesh is" . ih* 
generated automatically. Automatic mesh generation pro- 
ceeds by defining a modified quadtree (Figure VI-7) in an 
integer data space, then breaking the modified quadtree 
into a valid finite element mesh (Figure VI-8) . The final 
steps involve, first, pulling the boundary nodes from 
their locations in the integer data space and moving 
them precisely to the geometric boundaries and, second, 
smoothing the internal node locations to improve the 
element aspect ratios. Figure VI-9 shows the final mesh 
for the plate with a hole in it, after pulling and smooth- 
ing. Figure VI-10 shows a mesh generated for a piston cross 
section, while Figure VI-11 shows a mesh generated for a 
typical bracket. 

4 . Plans for Upcoming Period 
During the next reporting period, programming on the 
extrusion preprocessor and shell surface program will continue. 
As mentioned in the section in this report dealing with the 
numerical investigation of the micromechanics of composite 
fracture. Part III-E, the modified quadtree mesh-generating 
capability will be considered for use m the analysis method 
for following the process of crack propagation. 



3 



Figure VI-9. Final Triangular Mesh 
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Figure VI -10. Piston Cross Section 
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INTRODUCTION 


Technical meetings, on- and off-campus, provide impor- 
tant opportunities for interchange of technical information. 
Because of the large number of composites meetings, a cen- 
tral catalog with all upcoming meetings is being maintained 
and distributed periodically. In this way we help to assure 
that a Rensselaer faculty/staff member can participate in im- 
portant meetings. The calendar for this reporting period is 
shown in Table VII-1. Meetings attended by RPI composites 
program faculty/staff /students during the reporting period 
are shown in Table VII-2. Some meetings particularly rele- 
vant to composites, held on-campus with special speakers, 
are listed in Table VII-3 . A list of composite-related vis- 
its to relevant organizations by RPI faculty/staff /students , 
with the purpose of each visit outlined, is presented in 


Table VII-4. 
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TABLE VI I -1 

CALENDAR OF COMPOSITES -RELATED MEETINGS 
(September 30, 1981 through April 30, 1982) 


1981 

10/5-7 

10/7 

10/9 

10 / 10-11 

10/12-15 

10/13-15 

10/19-22 

10/20 

10/27 

10/28-30 

11/12-13 

11/15-20 

11/16-19 

11/23-24 

11/3Q-12/3 


International Symposium on Aero elasticity 7 , 
Nuremberg, Fed. Rep. of Germany. 

Symposium III - Design and Use of Kevlar in Air- 
craft, Los Angeles, CA. "Sponsored by DuPont." 

Symposium III - Design and Use of Kevlar in Air- 
craft, Philadelphia, PA. "Sponsored by DuPont." 

Annual Materials Science Seminar, Louisville, KY. 
"Sponsored by American Society of Metals." 

3rd World Congress on Finite Element Methods, 

Los Angeles, CA. 

13th National Technology Conference, Mt. Pocono, 
PA. "Sponsored by SAMPE." 

International Symposium on Optimum Structural 
Design, Tucson, AZ . 

Symposium on Producibility and Quality Assurance 
of Composite Materials, St. Louis, MO. "Spon- 
sored by ASTM. " 

ASCE Annual Convention and Exposition, St. Louis, 
MO. 

7th Annual Mechanics of Composites Review, Day- 
ton, OH. "Sponsored by Air Force Wright Aero- 
nautical Labs . " 

NASA Workshop on "High Temperature Nonlinear Con- 
stitutive Equations", Lewis Research Center, 
Cleveland, OH. 

ASME Winter Annual Meeting, Washington, DC. 
"Sponsored by ASME," 

Composites Symposium on Fabrication, Design and 
Testing, Wilmington, DE. "Sponsored by SPE and 
SESA." 

Annual Review of the Ceramic Materials Research 
Program on Brittle Material Design, Seattle, WA. 
"Sponsored by NASA and the University of Washing- 
ton . " 

13th Annual Department of Defense Manufacturing 
Technology Conference, San Diego, CA. "Sponsored 
by Uniter States Department of Defense." 
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1981 
12/7-9 

1982 
1/6 

1/11-13 

1/11-15 

1/18-22 

3/2-5 

3/9-10 

3/22-26 

3/24-26 

3/29 

3/29-4/1 

4/7-9 

4/28 


V/STOL Conference, Palo Alto, CA. "Sponsored 
by AIAA." 


American Chemical Society and Society of Plastics 
Engineers, Albany, NY 

20th Aerospace Sciences Meeting, Orlando, FL. 
"Sponsored by AIAA." 

37th Annual Conference and Product Showcase of 
the SPI Reinforced Plastics/Composites Institute, 
Washington, DC. "Sponsored by the Society of 
the Plastics Industry, Inc." 

Gordon Conference on Composite Materials, Ven- 
tura, CA. 

Design Criteria for the Future of Flight Con- 
trols, Dayton, OH. "Sponsored by Flight Dynam- 
ics Lab. " 

Symposium on Long-Term Behavior of Composites , 
Williamsburg, VA. "Sponsored by ASTM. " 

Conference on Computers /Graphics in the Building 
Process, Washington, DC. "Sponsored by National 
Academy of Science and World Computer Graphics 
Association. " 

International Conference on Forward Swept Wing 
Aircraft, Bristol, U. K. "Sponsored by the 
University of Bristol." 

International Conference on Deformation, Yield 
and Fracture of Polymers, Churchill College, 
Cambridge , England . 

American Chemical Society Meeting, Las Vegas, NV. 

11th Southeastern Mechanics Conference, Hunts- 
ville, AL. 

ASTM Task Group E 24.04.09 Meeting on Crack 
Growth in Adhesive Joints, Philadelphia, PA. 



282 


TABLE VI I -2 

COMPOSITES -RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
(September 30, 1981 through April 30, 1982) 


1981 


10/9 Symposium III - Design and Use of Kevlar in Air- 

craft (Prof. Diefendorf and Mr. Paedelt) , Phila- 
delphia, PA. 

10/12-15 3rd World Congress on Finite Element Methods (Prof 
Shephard) , Los Angeles, CA. 

Professor Shephard presented the paper, "The 
Finite Element Modeling Process - Will it be 
Automated " . 


10/27 ASCE Annual Convention and Exposition, (Prof. 

Shephard) , St. Louis, MO. 

Professor Shephard presented the paper, "Adap- 
tive Element Analysis Accounting for Material 
Interfaces" . 


10/27-29 7th Annual Mechanics of Composites Review (Prof. 
Loewy) , Dayton, OH. 

Professor Loewy, Keynote Speaker, presented 
"Composites: Coming of Age". 

11/12-13 NASA Workshop on "High Temperature Nonlinear Con- 
stitutive Equations" (Prof. Krempl) , Cleveland, OH 

11/15-20 ASME Winter Meeting (Profs. Krempl and Scarton) , 
Washington, DC. 

Professor Scarton presented the paper, "Acoustic 
Emission in Composite Materials". 


1982 

1/6 ACS and Society of Plastics Engineers Meeting 

(Prof. Diefendorf) , Albany, NY. 

Professor Diefendorf presented the paper, "Aero- 
space Applications of Composite Materials". 

1/18-22 Gordon Conference on Composite Materials (Prof. 
Diefendorf) , Ventura, CA. 

Professor Diefendorf presented the paper, 
"Residual Stresses in Carbon Fibers". 

3/28-4/1 American Chemical Society Meeting (Prof. Diefen- 
dorf) , Las Vegas, NV. 

Professor Diefendorf presented the paper, "In- 
terfaces and Liquid Crystalline Polymers". 
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3/29 

4/7-9 

4/28 


International Conference on Deformation, Yield and 
Fracture of Polymers (Prof. Sternstein) , Churchill 
College, Cambridge, England. 

Professor Sternstein presented the paper, 

"Matrix Dominated Mechanical Properties of 
-High Performance Composites". 

11th Southeastern Mechanics Conference, Prof. 
Krempl) , Huntsville, AH. 

ASTM Task Group E 24.04.09 Meeting on Crack Growth 
in Adhesive Joints (Prof. Shephard) , Philadelphia, 
PA. 
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TABLE VII -3 

COMPOSITES -RELATED MEETINGS/TALKS HELD AT RPI 
(September 30, 1381 through April 30, 1982) 


Topic 

Composites Applications 
to Rotary Wing Aircraft 

Kinking of Fibrous Com- 
posites 


The Optical Method of 
Caustics as Applied to 
The Fracture of Duc- 
tile Materials 


Date Speaker (s) 

10/6/81 Erv7in Durchlaub 

Boeing Vertol, Philadel- 
phia 

10/20/81 Bernard Budiansky 

Gordon McKay Professor of 
Structural Mechanics 
Harvard University 

2/16/82 Ares J. Rosakis 

Division of Engineering 
Brown University 


Design and Automated 2/23/82 

Manufacturing of Com- 
posite Airframe Struc- 
tures 


James O. Warren and 
David L. Smith 
McDonnell Douglas Corp. 


Pulse Propagation in 3/2/82 

Composites Viewed as 
Interpenetrating Solid 
Continua 


Evidence for Crack Tip 3/4/82 
Controlled Stable 
Crack Growth 


Continuum Damage Me- 3/9/82 

chanics (A Systematic 
Life Prediction Method) 

RPI Faculty Workshop on 3/15/8: 
Instructional Computer 
Graphics 


Matthew F. McCarthy 
Department of Engineering 
Science and Mechanics 
Virginia Polytechnic In- 
stitute and State Univer- 
sity 

George Green 
General Electric R & D 
Center, on leave from 
Central Electricity Gener- 
ating Board 

F. A. Leckie 

Professor, University of 
Illinois 

M. S. Shephard 
Assistant Professor of 
Civil Engineering 
Rensselaer Polytechnic 
Institute 


Finite Element Modeling 3/24/82 
Design 


M; S. Shephard 
Assistant Professor of 
Civil Engineering 
Rensselaer Polytechnic 
Institute 
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Fatigue Damage Meehan- 4/2/82 
isms in Carbon Fiber 
Reinforced Plastics 

Fibrous Skins 4/6/82 

An Analysis of the In- 4/26/82 
fluences of Time and 
Loading Rate on Crack 


A Finite Element Time 4/27/82 
Integration Method for 
the Theory of Visco- 
plasticity Based on In- 
finitesimal Total 
Strain and Over stress 


Karl Schulte 
Research Associate 
Institut f. Werkstoff- 
Fors chung, W. Germany 

S. A. Wainwright 
Professor, Duke University 

Mark ’M. Little 
Doctoral thesis disertation 
Department of Mechanical 
Engineering , Aeronautical 
Engineering and Mechanics 

Robert M. Zirin 
Doctoral thesis dissertation 
Department of Mechanical 
Engineering, Aeronautical 
Engineering and Mechanics 
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TABLE VI I -4 

COMPOSITES -RELATED VISITS TO RELEVANT ORGANIZATIONS 
by RPI Faculty /Staff /Students 
(September 30, 1981 through April 30, 1982) 


Visited 

Date 

By 

Purpose 


1981 



Carnegie Insti- 
tute of Technol- 
ogy, Pittsburgh, 
PA 

10/20 

Prof. M. S. 
Shephard 

Presented at Nippon 
Univac Kaisha Semi- 
nar: "Academic Re- 
search in Cad/Cam" 

Norton Company, 
Worchester , MA 

10/22 

Prof. R. J. 
Dief endorf 

Presented talk: 
"Ceramic Composites" 

General Electric 
Company, Sara- 
toga , NY 

10/28 

Prof. R. J. 
Dief endorf 

Presented talk: 
"Composite Materi- 
als" 

Du Pont Pioneer- 
ing Research Lab., 
Wilmington, DE 

11/3 

Prof. R. J. 
Dief endorf 

Presented talk: 

"The Structure of 
Carbon Fibers and 
the Consequences on 
Performance" 

Louisville, KY 
(S. P. E. Short 
Course) 

11/4 

Prof. S. S. 
Sternstein 

Gave invited lecture: 
"Matrix Dominated 
Mechanical Properties 
of Composites" 

Exxon Enterprises, 
Fountain Inn, SC 

11/4-5 

Prof. R. J. 
Dief endorf 

Presented talk: 

"The Structure of 
Carbon Fibers and the 
Consequences on Per- 
formance 

CIBA-Geigy Re- 
search Center, 
Tarry town, NY 

11/19 

Prof . S . S . 
Sternstein 

Presented talk: 
" Composites" 

Naval Research 
Laboratory, Wash- 
ington, DC 

12/1 

Prof. S. S. 
Sternstein 

Presented talk: 
"Viscoelastic Charac- 
terization of Compos- 
ites" 


1982 



Physical Acous- 
tics Corporation, 
Montreal, Canada 

1/5 

Prof. H. 
Scarton 

Presented talk: 
"Acoustic Emission in 
Composite Materials" 
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UCLA, Los Angeles, 
CA 

Aerospace Corpor- 
ation, El Segundo, 
CA 


Lehigh University, 
Bethlehem, PA 


General Motors Re- 
search Laboratory, 
Detroit, MI 

l 

Kent State, 

Kent, OH 
(Liquid Crystal 
Institute) 

Cornell Universi- 
ty, Ithaca, MY 
(Program for 
Computer Graphics) 

Du Pont Corpora- 
tion, Wilmington, 
DE 

Iowa State Uni- 
versity, Ames, IA 


1982 


1/12-13 

Prof. R. J. 
Diefendorf 

1/14 

Prof. R. J. 
Diefendorf 

3/1 

Prof. M. S. 
Shephard 

3/8 

Prof. M. S. 
Shephard 

3/17 

Prof. R. J. 
Diefendorf 

3/19 

Prof. M. S. 
Shephard 

4/18 

Prof. R. G. 
Loewy 

4/21-22 

Prof. M. S. 
Shephard 


Presented talk: 
"Composite Materials" 

Presented talk: 
"Structure of Carbon 
Fibers and Conse- 
quences on Perform- 
ance 

Presented at seminar: 
"Computer Graphics 
in Education and Re- 
search" 

Presented talk: 
"Progress on Auto- 
mated 3-D Mesh Gener- 
ation" 

Presented at seminar : 
"Discotic Liquid 
Crystals in Pitch" 

Presented talk: 
"Automatic Finite 
Element Mesh Gener- 
ation" 

To Discuss Composites 
Program with Dr. A. 
Dhingra 

To attend : "Work- 

shop on Computer- 
Aided Design in Edu- 
cation" 
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Preceding page blank 


PERSONNEL 


Co-Principal Investigators 

Ansell, George S., Ph.D. 
Loewy, Robert G., Ph.D. 
Wiberley, Stephen E., Ph.D. 

Senior Investigators 

Brunelle, E. J., Jr., SC.D. 
(Aeroelastic and structural 
design and analysis, applied 
mechanics of composite 
structures) * 

Bundy , F. , Ph.D. 

{Physical chemistry and 
structures testing) * 

Diefendorf 1 , R. J., Ph.D. 
(Fabrication, resin matrix, 
f,iber behavior, interfaces)* 

Feeser 1 , L. J., Ph.D. 
(Computer applications and 
graphics, computer aided 
design, optimization) * 

Goetschel, D. B. , Ph.D. 
(Structural analysis design 
and testing ) * 

Hagerup , H . J . , Ph.D. 
(Aerodynamics, configura- 
tion, pilot accomodation, 
flight testing) * 

Krempl, E., Dr.Ing. 

(Fatigue studies, failure 
criteria). * 

Scar ton, H., Ph.D. 

(Acoustic emission NDE)* 


Dean, School of Engineering 
Institute Professor 
Professor of Chemistry 


Associate Professor of 
Aeronautical Engineering 


Research Professor of Mate- 
rials Engineering 

Professor of Materials 
Engineering 

Professor of Civil Engineer- 
ing 


Assistant Professor of 
Mechanical Engineering 

Associate Professor of 
Aeronautical Engineering 


Professor of Mechanics and 
Director of Cyclic Strain 
Laboratory 

Associate Professor of Me- 
chanical Engineering and Me- 
chanics 


Fields of Speciality 

^Member of Budget Committee together with Co-Principal Inves 
tigators 
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Senior Investigators 

Shephard, M. S., Ph.D. 
(Computer graphics , finite 
element methods) * 

Sternstein 1 , S. S., Ph.D. 
(Failure analysis, matrix 
behavior, moisture effects) * 


Associate Director, Center for 
Interactive Computer Graphics 
and Assistant Professor of 
Civil Engineering 

William Weightman Walker 
Professor of Polymer- Engineer- 
ing 


Research Staff 

Manager & Master Technician, Composites Laboratory 
Paedelt, Volker 


Research Associate 

Muser, Christoph, Dr. Eng. 

Graduate Assistants 

Bobal , Gail, M.E. 

Chang, Chi-Mm, M.S. 

Chen , Chikuang , B . S . 

Chen , Kuang - j ung , B . S . 
Chen, Shu-ping, B.S. 
Helmer, James, B.S. 

Lumban Tobing, Frida, M.S. 
Niu , Tyan-Min , M.S. 

Undergraduate Assistants - 

Adel , Paul 
Applewhite, Keith 
Baxter, Scott 
Bergman , Richard 
Bertolazzi, Andrew 
Bristol, Brian 


Ogunlan, Fola, B.S. 
Singh, Sachchida, B.S. 
Uzoh, Cyprian, B.S. 
Valicenti , Raymond , B.S. 
Wmckler, Steven, M.E. 
Yehia, Nabil, M.S. 
Yurgartis, Steven, M.S. 

Seniors 

Chin , Hong 
Cox , Mary 
DeTaranto, Francis 
Dhimitri, William 
Fairchild, Kendall 
Fisher, Steven 


Fields of Speciality 

1 Member of Budget Committee together with Co-Principal Inves 
tigators 
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Newmark, Glenn 
Niederer, Melvin 
Northrop, Steven 

Undergraduate Assistants 

Baucom, Allan 
Benson, Robert 
Burdett, John 
Cioffoletti, Anthony 
DeMint, Thomas 
Ehrgott, Darhl 
Ficarra, Fobert 
Krrschenbaum, Susan 
Koch, Stephen 
Kuntsmann, Debra 
Lafreniere, Suzanne 
LeBlanc , James 

Undergraduate Assistants 


Stamorski, Anthony 
Vaney , Phil ippe 

Juniors 

Lundquist, Cheryl 
Malusa , Stephen 
Melendez, Hermenegildo 
Powell, David 
Quinlivan, Patrick 
Sirkxn, Stephen 
Spinelli, Gregory 
Taf finder, Douglas 
Toomey, John 
West, Jeffrey 
White, Eric 

- Sophomores 


Tobiason, Scott 
Goguen, Gary^ 
Hsi, Ten-Li 


Martin, Randall 
Schroeder, Hans 



294 


! 

/ 



( 



AUTHOR INDEX 


! Preceding page blank 1 295 


Page (s) 

Brunelle, E. J., Jr 71 

Bundy, F . 249, 251 

Diefendorf, R. J 9, 249, 251 

Goetschel, D. B 173 , 183, 199 

Hagerup, H. J . 249, 251 

Krempl, E 79 

Loewy, R. G 135 

Muser, C 129, 251 

Scarton, H. A 205 

Shephard, M. S 101, 119, 267 

Sternstein, S, S. 89 


' ,l > > I i I ri' ’WO 


, i f „* s:<ji hj,wi> 


I ’it'it 



